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Abstract. We review recent studies of the evolution of collective excitations in 
atomic nuclei far from the valley of /3-stability. Collective degrees of freedom govern 
essential aspects of nuclear structure, and for several decades the study of collective 
modes such as rotations and vibrations has played a vital role in our understanding 



' of complex properties of nuclei. The multipole response of unstable nuclei and the 



possible occurrence of new exotic modes of excitation in weakly-bound nuclear systems, 
present a rapidly growing held of research, but only few experimental studies of these 
phenomena have been reported so far. Valuable data on the evolution of the low- 
energy dipole response in unstable neutron-rich nuclei have been gathered in recent 
experiments, but the available information is not sufficient to determine the nature of 
observed excitations. Even in stable nuclei various modes of giant collective oscillations 
had been predicted by theory years before they were observed, and for that reason it 
is very important to perform detailed theoretical studies of the evolution of collective 
modes of excitation in nuclei far from stability. We therefore discuss the modern 
theoretical tools that have been developed in recent years for the description of 
collective excitations in weakly-bound nuclei. The review focuses on the applications 
of these models to studies of the evolution of low-energy dipole modes from stable 
nuclei to systems near the particle emission threshold, to analyses of various isoscalar 
modes, those for which data are already available, as well as those that could be 
observed in future experiments, to a description of charge-exchange modes and their 
evolution in neutron-rich nuclei, and to studies of the role of exotic low-energy modes 
in astrophysical processes. 
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1. Introduction 

Studies of nuclear vibrational modes over several past decades have provided valuable 
information on the structure of the nucleus and the forces of cohesion that are responsible 
for the nuclear binding. Collective degrees of freedom characterize many aspects of 
nuclear structure. The response of a nucleus to external forces often exhibits a degree of 
simplicity associated with collective modes: rotations and vibrations. Even the simplest 
excitations, giant vibrations or giant resonances characterized by a coherent oscillation 
of all the nucleons, probe not only global nuclear properties such as the size, the shape, 
the distributions of protons and neutrons, the compressibility of nuclear matter, but also 
the details of the in-medium modification of the nucleon-nucleon interaction, and the 
interplay between different degrees of freedom in dissipative processes that determine the 
damping mechanism. The excitation energies and decay properties of giant resonances 
have been measured mostly for nuclei along the valley of /^-stability, and the extension 
of these studies to regions of unstable nuclei is still in its infancy. 

The multipole response of nuclei far from the /3-stability line and the possible 
occurrence of exotic modes of excitation presents a rapidly growing field of research. 
Characteristic ground-state properties (weak binding of the outermost nucleons, 
coupling between bound states and the particle continuum, nuclei with very diffuse 
neutron densities, formation of neutron skin and halo structures) will also have a 
pronounced effect on the multipole response of unstable nuclei. For instance, the dipole 
(El) response of neutron-rich nuclei is characterized by the fragmentation of the strength 
distribution and its spreading into the low-energy region, and by the mixing of isoscalar 
and isovector components. While in light nuclei the onset of dipole strength in the 
low-energy region is caused by non-resonant independent single-particle excitations of 
the loosely bound neutrons, several theoretical analyses have predicted the existence 
of the pygmy dipole resonance (PDR) in medium- mass and heavy nuclei, i.e., the 
resonant oscillation of the weakly-bound neutron skin against the isospin saturated 
proton-neutron core. The interpretation of the dynamics of the observed low-energy 
El strength in nuclei with a pronounced neutron excess is currently very much under 
discussion. 

Of course, not only pygmy states, but also other possible exotic modes are being 
studied: isoscalar dipole, toroidal, giant pairing vibrations, low-energy monopole and 
quadrupole, and spin-isospin excitations in unstable nuclei. The isoscalar giant dipole 
resonance (ISGDR) corresponds to a second order high-energy compression mode 
and therefore provides information on the nuclear matter compression modulus, but 
the existence of a low-energy component has recently been experimentally confirmed. 
Several theoretical studies have predicted that the low-energy isoscalar dipole vibration 
is not sensitive to the nuclear compressibility and that, in fact, it could correspond 
to the toroidal dipole resonance. The toroidal dipole mode is a transverse zero-sound 
wave and its experimental observation would invalidate the hydrodynamical picture of 
the nuclear medium, since there is no restoring force for such modes in an ideal fluid. 
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Pairing vibrations are induced in the nucleus by the addition or removal of a pair of 
neutrons, and can be associated with the fluctuation of the pairing field. High-energy 
collective pairing modes - Giant Pairing Vibrations (GPV), have been predicted and 
analyzed theoretically, but have never been observed in studies of reactions induced by 
beams of stable isotopes. Entirely new types of collective excitations might arise in 
nuclei near the particle emission threshold: di-neutron vibrations close to the neutron- 
drip line, and proton pygmy resonances in proton-rich nuclei. Several new theoretical 
approaches have recently been developed, which provide a fully microscopic description 
of low-energy collective excitations in weakly bound nuclei. This review presents an 
opportunity to compare the results and predictions of various models, and to discuss the 
development of modern theoretical tools based on the interacting shell-model, the time- 
dependent non-relativistic and relativistic self-consistent mean-field framework, and the 
extensions of the latter models beyond the mean-field approximation. 

Theoretical predictions of exotic modes have also prompted the design of 
experiments with radioactive beams, and a number of studies of low-energy multipole 
response in unstable nuclei have been reported in recent years. Low-lying El strength 
has been observed in neutron-rich oxygen isotopes, exhausting about 10% of the classical 
dipole sum rule below 15 MeV excitation energy. In heavier systems data have 
recently been reported on the concentration of electric dipole strength below the neutron 
separation energy in N = 82 semi-magic nuclei. The experimental information which 
is presently available, however, is not sufficient to determine the dominant structure of 
the observed states. The Sn isotopes present another very interesting example of the 
evolution of the low-lying dipole strength with neutron number. Very recently the dipole 
strength distribution above the one-neutron separation energy has been measured in the 
unstable 130 Sn and the doubly-magic 132 Sn. In addition to the giant dipole resonance 
(GDR), evidence has been reported for a PDR structure at excitation energy around 
10 MeV both in 130 Sn and 132 Sn, exhausting a few percent of the El energy-weighted 
sum rule. Obviously this is a rapidly expanding field and many new experiments are 
being planned and designed at existing or future radioactive-beam facilities, which will 
allow the study of the evolution of collective modes in nuclei far from stability, and the 
discovery of new exotic modes of excitation. 

Besides being intrinsically interesting as new structure phenomena, exotic modes 
of excitation might play an important role in nuclear astrophysics. For example, the 
occurrence of the PDR could have a pronounced effect on neutron capture rates in 
the r-process nucleosynthesis, and consequently on the calculated elemental abundance 
distribution. Even though its strength is small compared to the total dipole strength, 
the PDR significantly enhances the radiative neutron capture cross section on neutron- 
rich nuclei, as shown in recent large-scale QRPA calculations of the El strength for 
the whole nuclear chart. The latest theoretical and computational advances in nuclear 
structure modeling have had a strong impact on nuclear astrophysics. More and more 
often calculations of stellar nucleosynthesis, nuclear aspects of supernova collapse and 
explosion, and neutrino-induced reactions, are based on microscopic global predictions 
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for the nuclear ingredients, rather than on phenomenological approaches. The ability 
to model the Gamow- Teller response, for instance, is essential for reliable predictions 
of /3-decay rates in neutron-rich nuclei along the r-process path. The calculation 
of GT strength, however, can also be used to constrain the spin-isospin channel of 
energy density functionals. When approaching the neutron (proton) drip lines, an 
increasing fraction of the GT_ (GT + ) strength is found within the /3-decay window, 
and a consistent study of this phenomenon has yet to be carried out in the framework 
of microscopic self-consistent models. 

The low-energy El strength could also play a role in the photodisintegration of 
Ultra-High Energy Cosmic Rays (UHECR). Under the assumption that UHECR are 
extra-Galactic nuclei accelerated to energies up to 10 21 eV, their interaction with 
the 2.7 K Cosmic Microwave Background (CMB) leads to photoabsorption reactions, 
followed by nucleon emission. Recent calculations have shown that the photo- 
disintegration path proceeds through regions of unstable nuclei, and that the nucleon 
emission rate is very sensitive to the low-energy dipole strength. 

The ultimate exotic modes far from stability could be collective excitations in nuclei 
far beyond the drip-line. Such systems can be expected to exist in the inner crust of 
neutron stars, where nuclear clusters are immersed in a dilute gas of neutrons and 
electrons. Model calculations have predicted the existence of super-giant resonances 
(SGR) at very low energies, typically around 3 MeV, and exhausting more than 70% 
of the EWSR. The SGR can have a pronounced effect on the specific heat of the crust, 
and could therefore affect the cooling time of the neutron star. 

The most accurate description of nuclear vibrations is provided by the time- 
dependent mean-field theory, and thus we begin this review with an outline of the 
theoretical tools which are based on the self-consistent theory of small-amplitude 
vibrations (Sec. H}, and its extension beyond the mean-field approximation (Sec. [3]). An 
extensive review of recent studies of the evolution of low-energy dipole vibrations, and 
the possible occurrence of pygmy modes in nuclei far from stability, is presented in Sec. HJ 
Various isoscalar modes, those already observed in experiments, as well as those that so 
far have only been predicted in theoretical studies, are reviewed in Sec. [5j A discussion 
of charge-exchange modes and their evolution in neutron-rich nuclei is included in Sec. EJ 
The possible role of exotic low-energy modes in astrophysical processes is described in 
Sec. [7] and, finally, Sec. [8] contains the concluding remarks and ends with an outlook for 
future studies. 

2. Self-Consistent Theory of Small Amplitude Vibrations 

Modern nuclear structure theory has evolved from macroscopic and microscopic studies 
of phenomena in stable nuclei towards regions of exotic, short-lived nuclei far from 
the valley of stability, and nuclear astrophysics applications. The principal challenge 
is to build a consistent microscopic theoretical framework that will provide a unified 
description of bulk properties, nuclear excitations and reactions. 
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The ab-initio approach, which starts from accurate two-nucleon and three-nucleon 
interactions, adjusted to nucleon-nucleon scattering data and spectroscopic data on 
few-nucleon systems, respectively, provides the basis for a quantitative description of 
ground-state properties, excited states and transitions in relatively light nuclei with 
A < 16. Improved shell-model techniques, which employ accurately adjusted effective 
interactions and sophisticated truncation schemes, are used in large-scale calculations 
of structure phenomena in medium-mass nuclei, including properties which are relevant 
for astrophysical applications. The structure of heavy nuclei with a large number of 
active nucleons, however, is best described in the framework of self-consistent mean-field 
models. A vast body of data, not only in medium-heavy and heavy stable nuclei, but also 
in regions of exotic nuclei far from the line of /3-stability, has been successfully analyzed 
with mean-field models based on the Skyrme and Gogny non-relativistic interactions, 
and on relativistic meson-exchange effective Lagrangians. The self-consistent mean-field 
approach to nuclear structure represents an approximate implementation of Kohn-Sham 
density functional theory, which enables a microscopic description of the nuclear many- 
body problem in terms of a universal energy density functional. When compared with 
ab-initio and shell-model approaches, important advantages of the mean-field framework 
include the use of global effective nuclear interactions, the ability to describe arbitrarily 
heavy systems including superheavy nuclei, and the resulting intuitive picture of intrinsic 
nuclear shapes. 

The unique structure properties which characterize highly unstable nuclei as, for 
instance, the weak binding of the outermost nucleons and the coupling between bound 
states and the particle continuum, the modification of the effective nuclear potential and 
the formation of nuclei with very diffuse neutron densities, the occurrence of neutron 
skin and halo structures, will also affect the multipole response of these systems, and 
new modes of excitation could arise in nuclei at the limits of stability. Therefore a 
quantitative description of properties of ground and excited states in weakly-bound 
nuclei, and especially studies of exotic modes far from stability, necessitate using the 
time-dependent self-consistent mean-field framework. In this section we present an 
outline of the theoretical tools that have been employed in most studies included in 
this article. We start with the non-relativistic Hartree-Fock-Bogoliubov theory, extend 
this approach to the relativistic mean-field framework, and derive the (continuum) non- 
relativistic and relativistic quasiparticle random phase approximations in the small- 
amplitude limit of the self-consistent time- dependent mean-field theory. For a more 
detailed introduction we refer the reader to several excellent monographs [TJ [2j [3j HJ EJ [6] , 
and recent review articles E] . 

2.1. The Hartree-Fock-Bogoliubov Method with Effective Nuclear Forces 

In addition to the self-consistent mean-field single-nucleon potential, the inclusion of 
pairing correlations is essential for a quantitative description of structure phenomena in 
open-shell spherical and deformed nuclei. In weakly-bound systems far from stability, 
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in particular, the Fermi surface for one type of nucleons is found close to the particle 
continuum. The single-nucleon separation energies become comparable to the pairing 
gaps, and this results in the lowest particle-hole (ph) and particle-particle (pp) modes 
being embedded in the continuum. A unified and self-consistent treatment of both the 
mean-field and pairing correlations becomes necessary, and the coupling between bound 
and continuum states has to be taken into account explicitly. 

The Hartree-Fock-Bogoliubov (HFB) theory [1] |2] provides a unified description 
of ph- and pp-correlations in nuclei and, when the self-consistent HFB equations are 
formulated in coordinate space, allows for a treatment of continuum effects in the 
presence of pairing. In the HFB framework two average potentials are taken into 
account: the self-consistent Hartree-Fock field T which encloses all the ph correlations, 
and the pairing field A which sums up the effects of the pp interaction. The ground 
state of a given nucleus is described by a generalized Slater determinant |$) of single- 
quasiparticle self-consistent solutions of the HFB equations, and represents the vacuum 
with respect to independent quasiparticles. The quasiparticle operators are defined by 
the unitary Bogoliubov transformation of the single-nucleon creation and annihilation 
operators: 

=^2u lk ct + V lkCl , (1) 
i 

where Z7jfc, Vj*. are single-quasiparticle wave functions that satisfy the HFB equation. 
The index I denotes an arbitrary basis, for instance the harmonic oscillator states. In 
the coordinate space representation I = (r,a, r), with the spin-index a and the isospin 
index r. The HFB wave functions determine the hermitian single-particle density matrix 

p lv = m4 Cl \<s>) = (V*V T ) U >, (2) 

and the antisymmetric pairing tensor 

k w = ($|c,,c,|$) = (V*U T ) U >. (3) 
These two densities can be combined into the generalized density matrix 

* = U' !-,«)• W 
For a nuclear Hamiltonian of the form 

H = ei °l ° l 4 Vlm'l'm c t c m' c m c l' > (5) 

I ll'mm' 

where v denotes a general nucleon-nucleon interaction, the expectation value ($|//"|$) 
can be expressed as a function of the hermitian density matrix p, and the antisymmetric 
pairing tensor k. The minimization of this energy functional with respect to p and k 
leads to the Hartree-Fock-Bogoliubov equations 

h - X 



-A* 





(6) 
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The single-nucleon Hamiltonian reads h = e + T, and the two self-consistent potentials 
T and A are defined by 

f«' = Virn>l'mf>mm> , (7) 

mm 1 

and 

A«/ = viii mm ik mm i . (8) 

m<m' 

The chemical potential A is determined by the particle number subsidiary condition, in 
such a way that the expectation value of the particle number operator in the ground 
state equals the given number of nucleons. The column vectors denote the quasiparticle 
wave functions, and Ej. are the corresponding quasiparticle energies. 

In the framework of Kohn-Sham density functional theory (DFT) [91 [101 [EE], 
of which the self-consistent HFB represents a particular implementation, the nuclear 
many-body problem is defined in terms of a universal energy density functional. Self- 
consistent mean-field models approximate the exact energy functional, which includes all 
higher-order correlations, with powers and gradients of ground-state nucleon densities. 
Although it models the effective interaction between nucleons, a general density 
functional is not necessarily related to any given microscopic nucleon-nucleon potential, 
i.e. it is rather the density functional that defines the effective nuclear interaction. 
This means that in the DFT formulation of the HFB framework one does not start 
with a Hamiltonian defined by a two-body interaction as in Eq. (jHJ), but rather from the 
energy functional E[JZ] = E[p, k] that depends on the densities p and k. The generalized 
Hamiltonian 7i is then obtained as a functional derivative of the energy with respect to 
the generalized density: 

H = — = ( h A ) (9) 
n 5K [ -A* -h* J ' 1 } 

where the single particle Hamiltonian h results from the variation of the energy 
functional with respect to the hermitian density matrix p 

h = (10) 
5p 

and the pairing field is obtained from the variation of the energy functional with respect 
to the pairing tensor 

A = £ (ii) 

OK 

In a compact form the stationary HFB equation is given in terms of the generalized 
density: 

[H,n] = Q. (12) 

In principle, a universal energy density functional can be built as an expansion 
in terms of local densities and currents, including all terms allowed by the underlying 
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symmetries, and without direct reference to any specific nucleon-nucleon interaction. By 
employing global effective nuclear interactions, with a small set of parameters adjusted 
to reproduce empirical properties of symmetric and asymmetric nuclear matter, and bulk 
properties of few stable spherical nuclei, the current generation of self-consistent mean- 
field models has achieved a high level of accuracy in the description of ground states and 
properties of excited states in arbitrarily heavy nuclei, including rare isotopes with a 
large neutron to proton asymmetry. In the non-relativistic framework, in particular, two 
classes of effective nucleon-nucleon interactions have become standard in self-consistent 
Hartree-Fock and Hartree-Fock-Bogoliubov calculations. The first is the finite-range 
Gogny force [El [13]: 

2 

^Gogny(ri 2 ) = J2 e ' (ri2/H)2 ( W J + B A - HjP T ~ M 3 P a P T ) 

'ri + r 2 ~ 



+ t 3 (l + x P a )5(ri2) P a 

+ iW ls {bi + £2) • k f x 5{r 12 ) k (13) 

where P a = ~(1 + &i ■ a 2 ) is the spin-exchange operator, P T = |(1 + fx ■ f 2 ) the isospin 
exchange operator, r 12 = r x — r 2 , and k = — |(Vi — V 2 ). The interaction includes the 
sum of two Gaussians with space, spin and isospin exchange, the term which includes 
an explicit density dependence, and the spin-orbit term. Wj, Bj, Hj, Mj, pj, t 3 , xo, a 
and Wis are adjustable parameters of the interaction. 

The second class of effective interactions is based on the zero-range, momentum- 
dependent Skyrme force: 

^sk(r 12 ) = t (1 + x P a ) 5(r 12 ) 

+ l tl ( 1+Xl P a )(# 2 5(r 12 ) + 5(v 12 ) k 2 

+ t 2 (l+x 2 P a )ti-6(r 12 ) k 

+ \ h {l+x,P (T )5{v l2 )p a(i:r Yl 



+ iW (d- 1 + 6- 2 )-ti x 5(r 12 )k. (14) 

Standard Skyrme interactions include ten adjustable parameters which determine the 
central term, the velocity dependent terms, the density dependent term, and the spin- 
orbit term [TJ [TH EES, EES]. The Skyrme energy density functional can be derived from 
the Hartree-Fock expectation value of the zero-range momentum dependent two-body 
force Eq. (fl4l) . or it can be parameterized directly without reference to an effective two- 
body force [7]. In the latter case the universal functional contains systematically all 
possible bilinear terms in the local densities and currents (plus derivative terms) which 
are invariant with respect to parity, time-reversal, rotational, translational and isospin 
t r ansf or mat ions . 

The pairing field in Eq. ([8]) is determined by the effective interaction in the pairing 
channel. In applications of the HFB model with the Gogny force, the same effective 
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interaction is used both in the ph and pp channels, with the exception of the density- 
dependent zero-range term which, with the choice of the parameter Xq — 1, does not 
contribute to the pairing channel. There is no physical reason, however, to use identical 
terms of the energy density functional in the the mean-field and pairing channels. 
Skyrme-type forces, for instance, generally exhibit unrealistic pairing properties, and 
thus an additional effective pairing interaction has to be specified in Skyrme-HFB 
calculations of open-shell nuclei. A standard choice for the pairing interaction is a 
zero-range local force, often including an explicit density-dependence: 



fpair(ri 2 ) = — (1 - Pa) 



2 



1 — T] 



Pc 



<*(ri-r 2 ). (15) 



Depending on the value of the parameters rj, p c and /3, pairing is more active in the 
volume of the nucleus, or on its surface. The strength V is adjusted to reproduce the 
odd-even staggering of binding energies in selected isotopic chains. Usually this results in 
slightly different values of the pairing strengths for protons and neutrons, thus breaking 
the isospin invariance of the pairing energy functional. 

For nuclear systems with time-reversal invariance, the HFB method can 
be considerably simplified by employing the Bardeen-Cooper-Schrieffer (BCS) 
approximation. In the BCS approximation the pairing potential d, defined by the 
relation: 



A= T , (16) 




where A is the pairing field of Eq. ([SJ), is diagonal in the basis of the eigenstates of the 
mean-field Hamiltonian h 

dnfh ^nmdnrii hip n £ n ip n , (1 V 

and rh (n) denotes the time-conjugate partner of the single-particle state m (n). The 
resulting two-component HFB wave functions read U n = u n ip n and V n = v n (p n , and the 
occupation amplitudes (u n ,v n ) are determined by the gap equation 

(e n - \)(u 2 n -vl) + 2d nn u n v n = 0, (18) 

and the normalization condition: 

u 2 n + v 2 n = l. (19) 

In the case of well-bound nuclei close to the stability line, pairing correlations are 
often treated in the BCS approximation, with the strength of the pairing force adjusted 
to the experimental odd-even mass differences. This approach, however, presents only 
a poor approximation for weakly-bound nuclei far from the valley of /^-stability. In 
particular, for nuclei at the limits of particle stability (drip-line nuclei), the Fermi level 
lies close to the continuum, and the coupling between bound and continuum states has to 
be taken into account explicitly. The BCS model does not provide a correct description 
of the scattering of nucleonic pairs from bound states to the positive energy continuum, 
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and as a result levels high in the continuum become partially occupied. Including the 
system in a box of finite size, i.e. solving the HFB equations in coordinate space, leads 
to unreliable predictions for quantities that crucially depend on the size of the box, e.g. 
nuclear radii. The reason is that in the BCS approximation the pairing field does not 
vanish asymptotically. Thus for weakly-bound nuclei the full HFB theory, including the 
continuum, has to be employed [T71 [18]. 

The asymptotic behavior of the HFB wave function is determined by the physical 
condition that at large distances from a nucleus the mean field T(r) and the pairing 
field A(r) vanish. For a bound system (negative chemical potential A < 0), two 
distinct regions characterize the quasiparticle spectrum. This is illustrated in Figure 
[TJ Between and —A the quasiparticle spectrum is discrete and both the upper and 
lower components of the radial HFB wave function decay exponentially at r — > oo. 
In the continuum region above —A one finds two types of HFB quasiparticle resonant 




HF HFB 

Figure 1. Illustrative representation of single-particle (left) and single-quasiparticle 
(right) spectra. A is the chemical potential. The arrows depict the correspondence 
between the single-particle states and the single-quasiparticle states in the discrete 
and continuous regions of the spectra. 

states. First, those which correspond to the single-particle resonances of the mean field. 
The low-lying single-particle resonances are particularly important in the treatment 
of pairing correlations of weakly-bound nuclei, because the pairing interaction scatters 
pairs of nucleons into positive-energy states above the particle threshold. A second type 
of resonant states is specific to the HFB method, and originates from bound single- 
particle states which, in the absence of pairing correlations, are found at energies e < 2A. 
With the inclusion of the pairing field these bound states couple with the continuum 
single-particle states and, therefore, acquire finite widths. 

The HFB equations are usually solved by imposing box boundary conditions, 
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i.e. the HFB wave functions are assumed to vanish beyond some distance, usually 
chosen to be a few times the nuclear radius. The energy continuum is thus replaced 
by a spectrum of discrete energies, whose density increases with the size of the box. 
Continuum HFB calculations can be performed either in the complex energy plane by 
employing Green's function techniques (TH], or on the real energy axis p2l HE]. It is 
possible to treat the continuum exactly in HFB calculations on the real energy axis, by 
imposing the correct boundary conditions for the HFB wave functions 17% and [20J. 
Far from the nuclear potential, the Uk and Vk which belong to the single quasiparticle 
continuum (E^ > —A), take the asymptotic form of a scattering state and a function 
that exponentially decreases at infinity, respectively. These asymptotic forms must, of 
course, be matched with the corresponding radial functions in the inner region of the 
potential. By comparing the exact treatment of the continuum with different levels 
of approximation (discretization in a box, HF-BCS including the resonant part of the 
continuum), it has been shown that the different ways in which the coupling to the 
continuum is treated in HFB strongly affect the resulting pairing correlations in nuclei 
close to the particle drip-lines. 

2.2. The Relativistic Hartree-Bogoliubov model (RHB) 

Self-consistent mean-field models based on the relativistic mean-field theory (RMF) [2U 
[22], have been very successfully employed in analyses of a variety of nuclear structure 
phenomena, not only in nuclei along the valley of /5-stability, but also in exotic nuclei 
with extreme isospin values and close to the particle drip-lines. RMF-based models have 
reached a level of sophistication and accuracy comparable to the non-relativistic Hartree- 
Fock-Bogoliubov framework based on Skyrme and Gogny effective interactions [8]. 

Most applications have used the finite-range meson-exchange representation of the 
RMF theory, in which the nucleus is described as a system of Dirac nucleons coupled to 
effective mesons and the electromagnetic field. A medium dependence of the effective 
interaction can be introduced either by including non-linear meson self-interaction terms 
in the Lagrangian, or by assuming an explicit density dependence for the meson- 
nucleon couplings [8J. An alternative RMF representation is formulated in terms of 
point-coupling (PC) (contact) nucleon-nucleon interactions [23J, without the inclusion 
of meson fields. The medium dependence of the interaction can be taken into account 
by terms of higher-order in the nucleon fields, or it can be encoded into the effective 
couplings, i.e. in the strength parameters of the interaction in the isoscalar and isovector 
channels. On the phenomenological level, when applied in studies of properties of finite 
nuclei, the meson-exchange and nucleon point-coupling representations produce results 
of comparable quality. The applications reviewed in this article are based on the finite- 
range meson-exchange picture of effective nuclear interactions, and in the following we 
outline its basic features. 

The isoscalar scalar cr-meson, the isoscalar vector cu-meson, and the isovector vector 
p-meson form the minimal set of meson fields which, together with the electromagnetic 
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field, are necessary for a quantitative description of bulk nuclear properties. The model 
is defined by the Lagrangian density 

£ = Cn + C-m + Cint- (20) 

Cn is the Lagrangian of the free nucleon 

C N = ^{i^d^m),p , (21) 

m is the bare nucleon mass and i/j denotes the Dirac spinor. C m is the Lagrangian of 
the free a, u and p meson fields, and the electromagnetic field 

C m = \d„o-d»a - l -mla 2 - \% v n^ + l -mlu^ 

- Ir^BT + \mlp,j? - \F» V F»\ (22) 

with the corresponding masses m a , m u , m p . The field tensors fl^, R pv , F^ u read 
Qfiv = d^<jj v — d v ui p 

R^u = djiPv - dypp (23) 
F^v = d fJ- A l/ — dyA^ , 

where arrows denote isovector quantities. The minimal set of interaction terms is 
contained in 

C int = -^r CT cn/> - - - iiT^A^ , (24) 

where the vertices read 

T a = g a , r^ = ^ 7 M , ^ p =g p rY, T» = e^^Y, (25) 

with the coupling parameters g a , g^, g p and e. The phenomenological a-meson 
approximates a large attractive scalar field that results from complicated microscopic 
processes, such as uncorrelated and correlated two-pion exchange. The tu-meson 
describes the short-range repulsion between the nucleons, and the p-meson carries the 
isospin quantum number. The latter is required by the large empirical asymmetry 
potential in finite nuclear systems. Because of parity conservation there is no direct 
contribution from the pion field. The self-consistent RMF approach represents a 
particular realization of the relativistic Kohn-Sham density functional theory [H)l [TT] , in 
which one attempts to effectively include in the nuclear energy density functional effects 
which go beyond the Hartree approximation (Fock terms, short-range correlations, 
vacuum-polarization effects). The many-body correlations of the energy density 
functional can be represented by a medium dependence of the corresponding effective 
nuclear interaction. An effective density dependence can be included, for instance, 
through meson self-interaction terms. Over the years a number of non-linear meson- 
exchange interactions have been adjusted to the nuclear matter equation of state 
and bulk properties of a set of spherical closed-shell nuclei, and applied in the 
description of nuclear properties along the /3-stability line. One of the most successful 
phenomenological interactions of this type is the non-linear parameter set NL3 
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which has been employed in many studies of ground-state properties and collective 
excitations both in stable nuclei and in exotic systems far from the line of /9-stability. 

Another class of medium- dependent effective interactions is characterized by an 
explicit baryon-density dependence of the meson-nucleon vertices. Such an approach 
retains the basic structure of the relativistic mean-field framework, but can be more 
directly related to the underlying microscopic description of nuclear interactions. 
The functional form of the meson-nucleon vertices can be deduced from in-medium 
Dirac-Brueckner interactions, obtained from realistic free-space NN interactions, or a 
phenomenological approach can be adopted, with the density dependence for the a, 
10 and p meson-nucleon couplings adjusted to properties of nuclear matter and a set 
of spherical nuclei. The former represents an ab-initio description of nuclear matter 
and finite nuclei, and the corresponding density- dependent relativistic mean-field model 
has also been applied to asymmetric nuclear matter and exotic nuclei [25]. In the 
latter approach very accurate phenomenological density-dependent relativistic effective 
interactions have recently been adjusted [26], E3, [28], and employed in analyses of both 
bulk nuclear properties and collective excitations. A number of recent studies have 
shown that, in comparison with non-linear meson self- interactions, relativistic models 
with an explicit density dependence of the meson-nucleon couplings provide an improved 
description of asymmetric nuclear matter, neutron matter and nuclei far from stability. 

An extensive review of the relativistic extension of the HFB theory - the relativistic 
Hartree-Bogoliubov (RHB) model - including numerous applications, can be found in 
Ref. [8]. In this framework the ground state of an open-shell nucleus is determined by 
the single-quasiparticle solutions of the RHB equations (cf. also Eq. (jSj)): 



density of Eq. (f20l . A denotes the pairing field, and (Uk,Vk) are quasiparticle spinors. 
Pairing effects in nuclei are restricted to a narrow window of a few MeV around the 
Fermi level, and their scale is well separated from the scale of binding energies which 
are in the range of several hundred to thousand MeV. There is no empirical evidence for 
any relativistic effect in the nuclear pairing field and, therefore, pairing can be treated 
as a non-relativistic phenomenon. In most applications of the RHB model for spherical 
and deformed nuclei the Gogny force (cf. Eq. ffTBl ) has been employed in the pp channel. 

The RHB equations are solved self-consistently, with mean-field potentials in 
the particle-hole channel determined from solutions of the stationary Klein-Gordon 
equations 





(27) 




(29) 
(30) 
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for the cr-meson, the tu-meson, the p-meson and the photon field, respectively. Because 
of charge conservation, only the 3-rd component of the isovector p-meson contributes. 
For the ground-state solution of an even-even nucleus there are no contributions from 
currents (time-reversal invariance) and the spatial components u>, p 3 , A of the vector 
fields vanish. The source terms in equations (12711 to ( [301) are sums of bilinear products 
of baryon amplitudes 

Ps(r) ^#)A(r) (31) 

fc>0 
k>0 

Ps(r) = XX(r)7*V*(r) (33) 

fc>0 

p em (r) = ^^(r)^y*(r), (34) 

k>0 

where the summation is performed only over occupied orbitals of positive energy. This 
is the no-sea approximation, in which the Dirac sea of negative-energy states does not 
contribute to the densities and currents in an explicit way. The self-consistent solution of 
the Dirac-Hartree-Bogoliubov integro-differential equations and Klein-Gordon equations 
for the meson fields determines the ground state of a nucleus. 



2.3. Continuum QRPA 

Small amplitude collective excitations of arbitrarily heavy nuclei can be accurately 
described by the random phase approximation (RPA) or, in the case of open-shell nuclei, 
by the quasiparticle random phase approximation (QRPA) [TJ [2J [U El E] • As it has been 
already emphasized in the introduction to this section, a quantitative description of 
ground-states and excitations in weakly-bound nuclei characterized by the closeness of 
the Fermi surface to the particle continuum, must take into account the effects of the 
coupling between bound states and the particle continuum. Here we derive the QRPA 
response based on the continuum HFB description of the nuclear ground state in the 
coordinate space representation, which is naturally suitable for the treatment of the 
coupling to continuum states [291 [30]. The QRPA represents the small amplitude limit 
of the general time-dependent (TD) HFB theory, and we therefore start from the HFB 
equation which describes the response of the generalized density matrix to an external 
periodic perturbation pQ: 

ih^ = [H(t)+F(t),K(t)}. (35) 

TZ and 7i are the time- dependent generalized density and HFB Hamiltonian, 
respectively, and T is the external periodic perturbation given by 

T = Fe- iujt + h.c. (36) 
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In the presence of pairing correlations, the fluctuating field F in Eq. (1361) is a generalized 
one-body operator which includes both particle-hole and two-particle transfer operators 

F = E 4 lc fe + + Ffaci) , (37) 

ij ij 

and the single-particle creation and annihilation operators, respectively. We 

assume that the external field induces small oscillations around the stationary solution 
of the HFB equation (fl2|) . which we denote here as 1Z° (TC° being the corresponding 
Hamiltonian) . Accordingly, 

nit) = n° + n'e~ iujt + h.c. , (38) 
nit) = n° + n'e- iut + h. c ., (39) 

and the TDHFB equation ( 135]) becomes 

huW = [n', n°] + [H°, n'} + [f, n°] . (40) 

The variation of the generalized density has the form 

* ( ; ■ ) • w 

where = (0|c|q|') is the variation of the particle density, = (0|c,Cj|') and 
R'ij = (0\cjcl\') are the corresponding fluctuations of the pairing tensor, and |') denotes 
the change in the ground-state wavefunction |0) caused by the external field. Instead 
of the variation of just one quantity (p') as in RPA, in QRPA we must specify the 
variations of three independent quantities, namely p', k' and R' . If we use the notation 

( p'\ 

P'= «' , (42) 
\k' J 

for the density variations, the corresponding variation of the HFB Hamiltonian reads 

W = Vp' , (43) 

where V is the matrix of the residual interaction, expressed in terms of the second 
derivatives of the HFB energy functional E, defined in Sec. 12.11 

V^iyxt, rVV) = - , -, a, (3=1, 2, 3. (44) 

dpp(r'a'r')dp a (raT) 

Here the notation at means that whenever a is 2 or 3, then a is 3 or 2. In this three 
dimensional space, the first dimension represents the particle-hole (ph) subspace, the 
second one the particle-particle (pp), and the third one the hole- hole (hh) subspace. 

The QRPA Green's function G relates the perturbing external field to the density 
variations 

p' = GF , (45) 
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where F is the three dimensional column vector 

I FU \ 

F = F 12 . 

\ F21 J 

G is the solution of the Bethe-Salpeter equation 

G = (1 - GoV)- 1 Go = Go + G VG , 
and the unperturbed Green's function Go reads 



G a/3 (rar, rVV; 



U?}{rar)U^\v'a'r') m 2 (raT)U*f 2 (r'a'r') 



Huj-(Ei + Ej) +ir/ hw + (Ei + Ej) + it] 



(46) 



(47) 



(4? 



where Ei are the quasiparticle energies and Uij are 3 by 2 matrices expressed in terms 
of the two components U and V of the HFB wave functions 

/ U.iva^iva) U*{ra)V*{va) \ 
U i3 {jo)= Vi{ra)V^a) V* (ra)V* (ra) , (49) 
\ -K-M^-(m) -U*(ra)U*(ra) J 

with the notation f(ra)=—2af(r — a). The symbol in Eq. fT451) indicates that 

the summation is taken both over discrete and continuum quasiparticle states, i.e., the 
unperturbed Green's function Go takes into account the resonant states. This can be 
done either in the complex plane by integrating on a contour [30] , or on the real axis by 
integrating the resonant states with the corresponding widths [29] . 

Go in Eq. (|48l) is constructed from the solutions of the HFB equations, i.e. the 
quasiparticle energies and the corresponding wave functions U and V. For a given 
residual interaction in Eq. (1441) . G is then calculated from the Bethe-Salpeter equation 
(JIT]) . G can be used for calculating the strength functions associated with various 
external perturbations. For instance, in the case of transitions from the ground-state to 
excited states within the same nucleus, only the (ph,ph) component of G is acting. If 
the interaction does not depend on spin variables, the strength function is then given 
by 

L / 'i5 ,11 *(r)G 11 (r,r';a;)F 11 (r0drdr / (50) 



S(u) 



-Im 



7T 



where G 11 is the (ph,ph) component of the QRPA Green's function. Examples of such 
calculations can be found in Refs. [29l [30] . 

It should be noted that generally the notation includes the neutron-proton 
formalism, i.e. each supermatrix (Go, G, V) consists of nine blocks which correspond 
to the ph, pp, hh channels, and each of these blocks is further divided into four 
sub-blocks corresponding to the neutron-neutron, neutron-proton, proton-neutron and 
proton-proton channels, respectively. 

An important point that must be emphasized is the issue of self-consistency. 
Namely, since the characteristic ground-state properties of weakly-bound nuclei strongly 
influence the multipole response of these systems, it is particularly important that for 
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the QRPA residual interactions, both in ph and pp channels, the same interactions are 
used which also determine the ground-state solutions of the HFB equations. 

It is, however, difficult to achieve full self-consistency in continuum RPA or QRPA 
calculations with Skyrme forces. While the zero-range terms of the interaction do not 
present any particular problem, the velocity-dependent terms introduce serious technical 
difficulties that are often avoided by approximating the residual interaction in the 
(ph,ph) subspace by its Landau-Migdal limit, in which the momenta of the interacting 
particle and hole are equal to the Fermi momentum, and the transferred momentum is 
zero. Taking the Landau-Migdal form for the ph interaction simplifies the numerical 
task considerably, however at the cost of losing the full self-consistency. This deficiency 
of the QRPA approach is cured by renormalizing the residual interaction. It should 
be noted, however, that an exact treatment of continuum effects within the general 
framework of the fully self-consistent QRPA is presently not available. 

2.4. Discrete QRPA 

In addition to the linear response formalism based on Green's functions, the QRPA 
can also be formulated in a discrete basis. We therefore start from Eq. (|40p . and 
choose a discrete basis of two-quasiparticle configurations \kl) = c4aj |0), where ct* 
is a quasiparticle creation operator and |0) is the quasiparticle vacuum. The energy of 
the configuration \kl) is E k + Ei. The quasiparticle space consists of discrete states, even 
when it describes configurations in the continuum. This discretization is achieved either 
by enclosing the nuclear system in a finite box, or by expanding the HFB wave functions 
in a discrete basis, e.g., in terms of eigenfunctions of a harmonic oscillator potential. In a 
QRPA description of transitions to low-lying excited states in open-shell weakly-bound 
nuclei, in particular, the two-quasiparticle configuration space must include states with 
both nucleons in discrete bound levels, states with one nucleon in a bound level and one 
nucleon in the continuum, and also states with both nucleons in the continuum. 

The QRPA equation is obtained by taking the matrix elements of Eq. (1401) between 
(kl\ and |0), and then also between (0| and \kl). Inserting a completeness relation in 
terms of a sum over all configurations labeled by k'l', the QRPA matrix equation finally 
reads 

Aki,k'i' Bki,k'v \ ( Xk'v \ _ E ( Xki 
- B ki,k>i> - A *ki,k'i> ) \ Y k'i' J \ Yki 

A and B contain matrix elements of the HFB or HF-BCS Hamiltonian. In the simplest 
case in which pairing is treated in the BCS approximation, the explicit expressions for 
the QRPA matrices read 

Aki,k'v = (E k + Ei)5 k k'0~ii' + 

+ V kWV {u k UiU k 'U V + V k ViV k tV V ) + 
+ W kWV {u k ViUyV V + V k UiV k >Ui>), 

B k i tk >i> = - V k i k >i>(u k uiv k >vi> + v k viu k >u v ) + 
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+ W k l h r V {u k VlU k ,V V + V k UiV k <Ui<), (52) 

where V and W denote the matrix elements of the interactions in the pp and ph channels, 
respectively, and u and v are the corresponding BCS occupation factors of single-nucleon 
states. 

The solution of the eigenvalue problem of Eq. ( IBTl) determines the energies E n of 
the excited vibrational states |n), and the corresponding wave functions expressed in 
terms of the forward and backward amplitudes Xj^' and 

H = E( x H )a ^ t+ ^rw^)io). (53) 

kl 

If instead of the simpler HF-BCS model, the full HFB framework is used as a 
basis on which the QRPA is formulated, the matrix equation (I5ip remains valid, but 
the matrices A and B display a more complicated structure. In place of the simple 
u and v occupation factors which depend on a single index, one finds the Bogoliubov 
matrices U and V. One way to circumvent the technical difficulties of working in 
the quasiparticle basis is to formulate the QRPA in the canonical single-nucleon basis 
which diagonalizes the density matrix. The corresponding matrix equations require 
only the evaluation of matrix elements of the residual ph and pairing pp interactions 
in this basis, multiplied by certain combinations of occupation factors. The fully self- 
consistent QRPA, formulated in the HFB canonical single-particle basis, was introduced 
in Ref. [31] • Self-consistency requires the QRPA residual interaction to be derived from 
the same force or energy functional that determines the HFB solution. This is crucial for 
the decoupling of spurious states, associated with symmetry breaking by the mean field 
solution, from the spectrum of physical excitations. The details of the implementation 
of the fully self-consistent HFB+QRPA framework in the canonical basis, and accurate 
tests using Skyrme energy density functionals and density-dependent pairing functionals, 
have recently been reported in Ref. |32j . 

The principal source of arbitrariness in the matrix representation of QRPA is the 
truncation of the basis. In realistic applications it is, therefore, necessary to verify 
the stability of the results with respect to variations of parameters that determine 
the discretization and basis truncation (the size of the box in coordinate space or the 
harmonic oscillator parameter, and the upper limit for unperturbed energies of two- 
quasiparticle configurations). 

Vibrational states can be excited by acting on the nucleus with an external operator 
F(r). The corresponding strength function is defined as 

S(u) = ^2\(n\F\0)\ 2 6(tuv - E n ). (54) 

n 

If the explicit wave function of the state \n) (Eq. (|53|) ) is inserted in this expression, the 
strength of each peak can be evaluated in terms of the single-particle matrix elements 
of the operator F(r), the coefficients of the BCS or Bogoliubov transformations, and 
the forward X and backward Y amplitudes. 
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2.5. Relativistic QRPA 

Relativistic RPA calculations have been performed since the early 1980s, but it is 
only more recently that non-linear meson self-interaction terms or density- dependent 
meson-nucleon couplings have been included in the RRPA framework [33j [3H [35] . As 
in the case of ground-state properties, the inclusion of a medium dependence in the 
residual interaction is necessary for a quantitative description of collective excited states. 
Another essential feature of the RRPA is the fully consistent treatment of the Dirac sea 
of negative energy states. Within the no-sea approximation, in addition to the usual 
particle-hole pairs, the RRPA configuration space must also include pair-configurations 
built from positive-energy states occupied in the ground-state solution, and empty 
negative-energy states in the Dirac sea [36]. Collective excitations in open-shell 
nuclei can be analyzed with the relativistic quasiparticle random-phase approximation 
(RQRPA), which in Ref. [37] has been formulated in the canonical single-nucleon basis 
of the relativistic Hartree-Bogoliubov (RHB) model. An alternative derivation of the 
RQRPA in the response function formalism, but with pairing correlations treated only 
in the BCS approximations, has recently been formulated in Ref. |38j . 

The RQRPA represents the small amplitude limit of the time-dependent relativistic 
Hartree-Bogoliubov (RHB) framework. The RQRPA matrix equations in the 
quasiparticle basis are, however, rather complicated and require the evaluation of the 
matrix elements of the Dirac Hamiltonian in the basis of the Hartree-Bogoliubov spinors 
Uk{r) and Vfe(r). A considerably simpler representation is provided by the canonical 
single-nucleon basis. Namely, any RHB wave function can be expressed either in 
the quasiparticle basis as a product of independent quasiparticle states, or in the 
canonical basis as a highly correlated BCS-state. The canonical basis is specified by the 
requirement that it diagonalizes the single-nucleon density matrix. The transformation 
to the canonical basis determines the energies and occupation probabilities of single- 
nucleon states that correspond to the self-consistent solution for the ground state of a 
nucleus. Since it diagonalizes the density matrix, the canonical basis is always localized. 
It describes both the bound states and the positive-energy single-particle continuum. 

Taking into account the rotational invariance of the nuclear system, the matrix 
equations of the RQRPA read [37J: 



For each RQRPA energy u u , X u and Y v denote the corresponding forward and 
backward two-quasiparticle amplitudes, respectively. The coupled RQRPA matrices 
in the canonical basis read 




(55) 



A J 



■kk'W ~ 




7,7 



(56) 
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+ <kmi>(-lY l - il ' +J v£ kfl ,. (57) 

H 11 denotes the one-quasiparticle terms 

H li = ( u kUi - v k vi)h kl - (u k vi + v k u t )A kl , (58) 

i.e. the canonical RHB basis does not diagonalize either the Dirac single-nucleon mean- 
field Hamiltonian ho, or the pairing field A. The occupation amplitudes v k of the 
canonical states are eigenvalues of the density matrix. V and V are the particle- 
hole and particle-particle residual interactions, respectively. Their matrix elements are 
multiplied by the pairing factors ^ and £, defined by the occupation amplitudes of 
the canonical states. The relativistic particle-hole interaction V is defined by the same 
effective Lagrangian density as the mean-field Dirac single-nucleon Hamiltonian ho- V 
includes the exchange of the isoscalar scalar <r-meson, the isoscalar vector u;-meson, the 
isovector vector p-meson, and the electromagnetic interaction. The two-body matrix 
elements include also contributions from the spatial components of the vector fields. 
The pairing factors read 

vL'Vlx' for a ' u ° >/A if J is even 

for oj>, p, A; if J is odd . . 

^"■' AA ' ^ < V^Vxy for a, W °, p°, A ; if J is odd ! ' * 

for u>, p, A; if J is even 

with the ^-coefficients defined by 

Vkk> = ± v k u k > , (60) 

and 

^fcfc' = U k u k> =F v k v k > . (61) 

a, uj°, p°, and A denote the time-like components, and u), p, A the spatial components 
of the meson and photon fields, respectively. 

The RQRPA configuration space must also include the Dirac sea of negative energy 
states, i.e. pair-configurations formed from the fully or partially occupied states of 
positive energy and the empty negative-energy states from the Dirac sea. The inclusion 
of configurations built from occupied positive-energy states and empty negative-energy 
states is essential for current conservation and the decoupling of spurious states, as 
well as for a quantitative comparison with the experimental excitation energies of giant 
resonances. 

The RQRPA model is fully self-consistent: the same interactions, in the particle- 
hole and particle-particle channels, are used both in the RHB equations that determine 
the canonical quasiparticle basis, and in the RQRPA equations. The parameters of the 
effective interactions are completely determined by the RHB calculations of ground- 
state properties, and no additional adjustment is needed in the RQRPA calculations. 
This is an essential feature of the RHB+RQRPA approach and it ensures that RQRPA 
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amplitudes do not contain spurious components associated with the mixing of the 
nucleon number in the RHB ground state, or with the center-of-mass translational 
motion. 

2.6. Multipole Transition Strength and Transition Densities 

We conclude this section with a brief overview of the operators associated with nuclear 
collective excitations, and provide a summary of basic definitions and useful relations. 

In the simplest case the number of protons and neutrons does not change when 
an external perturbation acts on the nucleus. The one-body operator F(r) associated 
with the external field can then induce, even in a superfluid system, only ph excitations. 
Later on we will also consider more general cases, including pairing vibrations and 
charge-exchange excitations. 

The case which is best defined is the electromagnetic excitation of a nucleus with 
real photons. A detailed discussion of this process can be found in many textbooks, 
for instance in Refs. [U El [6]. The multipole decomposition of the photon field plane 
wave leads to terms which contain Bessel functions jj(kr), where k is the photon 
momentum. In the cases of interest here the photon energy is of the order of MeV, 
and r can be the size of the nuclear radius, so that kr m 10 -2 and the Bessel function 
is approximated by ^j+iy.\ • I 11 ^ ms a PP rox i ma tion the dominant part of the matrix 
element which corresponds to the electric multipole J, is proportional to 

J p(r)r J Y JM (f) dr (62) 

where p is the charge density. The transition amplitude is therefore the matrix element 
of the operator 

z 

Fj M (r) = Y,erfY JM (r i ) (63) 

i=l 

between the initial and the final state (the sum is over protons). For states with 
good angular momentum, the response to the electric multipole transition operator 
is described by the reduced transition probability 

B{EJ,J t -> J f ) = - J—K/HFjIIOI 2 . (64) 

This relation is, of course, also valid for magnetic multipole transition operators. 

In many cases studies of excited states with electromagnetic probes are not possible, 
either because forbidden by selection rules, like for example the isoscalar monopole 
resonance which has zero angular momentum and therefore cannot be excited by 
photons, or because strongly hindered by nuclear excitations. Transitions induced by 
the strong interaction are typically studied by means of hadron inelastic scattering. In 
analogy with the electromagnetic case, the strong field is described by a plane wave 
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and, in the limit of small momentum transfer, the nuclear multipole transition operator 
reads 

A 

F JM (r) = J2rfYjM(h) (65) 

i=l 

when neutrons and protons are excited in phase (isoscalar excitation), or 

A 

F JM (r) = J2 r i Y JM(h)r z (i), (66) 

i=\ 

when neutrons and protons are excited with opposite phases (isovector excitations), and 
where r z (i) is the third component of the isospin operator. 

There is no compelling theoretical justification for the choice of the "effective" 
nuclear operators of Eqs. fl65l) and fl66|) . In those cases when the transferred momentum 
is not small, excitations of "overtones" induced by the component of the operator 
proportional to r/ +2 are also observed. Nevertheless, it has been shown that properties 
of nuclear giant resonances can be consistently described with the excitation operators 
Eqs. fl65l) and fl66l) [5, 6J. Concerning the isospin degree of freedom, some nuclear 
reactions are selective and either isoscalar or isovector states are excited. This is the 
case, for instance, in (a, a') inelastic scattering which predominantly excites isoscalar 
states. In other cases, like (p,p') inelastic scattering, the reaction itself is not isospin 
selective but, if the target is a light nucleus with good isospin quantum number, then 
isoscalar and isovector excitations can be separated. In those nuclei where isospin is no 
longer a good quantum number, e.g. light neutron-rich nuclei far from stability, nuclear 
excitations will correspond to a mixture of isoscalar and isovector modes. 

When excitations include the spin degree of freedom, the corresponding isoscalar 
and isovector multipole operators read 

A 

F JM (r) = ^Trf [Y L {h) ® &{i)] JM , (67) 

i=l 

and 

A 

F JM (r) = J2 r il Y ^i)®t(i)}jMT z (i) , (68) 

i=l 

respectively, with J — L, L ± 1. 

The strength function S(u) associated with the transition operator F is defined by 
Eq.( l54l) . It is often useful to consider moments of the strength function: 

m k (F) = I (hio) k S(co)d(hiu) = J2 E t\(n\F\0)\ 2 5(hu - E n ). (69) 

The first moment, or the energy-weighted sum rule (EWSR), is very important because 
its value equals the ground-state expectation value of the double commutator [F, [H, F}}, 
and therefore it can be evaluated without actually calculating the strength function. 
This is the well known Thouless theorem, and its non-trivial extension to the HFB- 
QRPA case has been proven in Ref. |29j . The ratio mi /mo is the quantity often compared 
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with the experimental excitation energy of the corresponding resonance although, of 

course, this is only correct if there are no multiple peaks within the energy interval over 

which the integration in Eq. ( 1691) is performed. One often finds the notation Eq 

for m,i/m (y'mi/m-i). 

When F represents the isovector electric dipole operator, and the residual 

interaction does not include velocity-dependent or exchange terms, a very simple result 

is obtained for the EWSR: 

2ir 2 e 2 hNZ N Z 
mi = — ~ 60—- MeV ■ mb. 70 

mc A A v ' 

This is the well known Thomas-Reiche-Kuhn (TRK) sum rule, which only includes 
the numbers of neutrons and protons, and is therefore completely model independent. 
For Skyrme and Gogny forces, which exhibit velocity-dependent and exchange terms, 
respectively, the value of m\ in the above expression is multiplied by (1 + «), where k 
is typically ~ 0.2-0.3. In the general case of a spin-independent isoscalar operator of 
multipole L, one can derive the following result: 

mi = £^^- A(r(2i " 2)) ' (71) 
where (r^ 2i ~ 2 ^) denotes the ground-state expectation value. 

Essential information on the dynamics of a nuclear collective mode is contained in 
the transition density. For the state \v) this quantity is defined as the matrix element 
of the density operator: 

^(r) = (i/|^tf(r-p < )|0>, (72) 

i 

where |0) denotes the ground state. The proton (neutron) transition density includes 
summation over protons (neutrons) in Eq. ( 1721) . and the isoscalar (T = 0) and isovector 
(T = 1) transition densities are defined by: 

5pY = 5pY + (-lf5ff/. (73) 

Assuming spherical symmetry, the transition density reads 

5p^(r)=5p n J ^(r)Y J * M (r), (74) 

where the radial factor can be expressed in terms of the forward (X) and backward 
(Y) QRPA amplitudes, and the single-nucleon radial wave functions multiplied by the 
corresponding occupation factors: 

5p^(r) = {Xl l + Y^)(u k v l + v k u l ){k\\Y J \\l)R k (r)R l (r). (75) 

3. Beyond the Mean-Field Approximation 

3.1. Extensions of the (Q)RPA 

The mean field theories which have been discussed in the previous Section, cannot 
provide a complete description of the physical phenomena associated with collective 
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nuclear excitations. The stationary models based on effective interactions, e.g. the 
HF or HFB, certainly include a large amount of correlations in their phenomenological 
parameters, and predict nuclear binding energies with remarkable accuracy. However, 
these models fail to reproduce the empirical features of single-particle levels around 
the Fermi surface [391 1^0] . which is essential for a complete description of collective 
excitations. 

Time- dependent theories, and their small amplitude limit such as the RPA and 
QRPA, provide an accurate description of the two principal integral characteristics 
of giant resonances: the total strength m and the energy- weighted sum rule mi. 
Accordingly, the calculated centroid energies mi /mo will not be modified by including 
additional correlations. On the other hand, the second moment of the strength 
distribution, which in the case of a single peak is associated with its full width at half- 
maximum (FWHM), or other details of the response function, cannot be described at 
the level of a mean-field theory. (Q)RPA calculations do not reproduce the experimental 
values of the total width of vibrational states. 

At the (Q)RPA level, nuclear collective motion is represented as a coherent 
superposition of Ip — lh (or two-quasiparticle) states. The energy and angular 
momentum of these vibrations can be released to other degrees of freedom, because 
vibrational states are embedded in a dense background of excited states. When the 
energy of a vibrational state lies above the particle emission threshold, the state can 
decay by neutron or proton emission. This damping mechanism is associated with the 
escape width T\ which can be taken into account within the framework of continuum- 
(Q)RPA (see Sec. 12. 3j) . The spreading width arises because the energy and angular 
momentum of coherent vibrations can be transferred to more complicated nuclear states, 
of 2p — 2h (and eventually 3p — 3h ... np-nh) character. In order to describe a 
theoretical framework must include the coupling to these complex configurations. 

A vast amount of data on widths of giant resonances in stable nuclei has been 
accumulated. Escape widths are large in light nuclei, but become less important in 
medium-heavy and heavy systems, where typical values are of the order of 1 MeV or 
less. On the other hand, spreading widths of several MeV are typical for resonances in 
heavier nuclei. For an extensive review on this subject the reader is referred to Ref. [39J. 
Much less is experimentally known about transition strengths of giant resonances in 
unstable nuclei, however, it has been pointed out that in neutron-rich nuclei far from 
stability spreading widths could be enhanced with respect to stable nuclei [U]. 

In the following we will briefly review several models which go beyond the mean-field 
approximation in the description of collective excitation phenomena. Simple lp — lh, 
or two-quasiparticle states are coupled by the residual interaction to more complicated 
configurations. A straightforward approach to this problem would be the diagonalization 
of the effective Hamiltonian in the configuration space which includes at least the 2p—2h, 
or four-quasiparticle states. This avenue, however, is usually not feasible because the 
number of 2p — 2h configurations can become very large (~ 10 2 or 10 3 per MeV). 
The second RPA (SRPA) is based on non-interacting 2p — 2h configurations and, in 
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Figure 2. Diagrams which correspond to the coupling of the p — h components of a 
giant resonance with phonon states. 

addition to interaction terms acting in the lp — lh space, only interactions between 
2p — 2h and lp — lh configurations are explicitly taken into account. The SRPA 
retains basic properties of the RPA, e.g. the conservation of the EWSR [12]. Practical 
implementations of this framework have been restricted to relatively light nuclei, and 
only with radical truncations of the configuration space (for a detailed discussion, cf. [13] 
and references therein). In addition, the current SRPA models are not self-consistent, 
e.g. in Ref. [13] the Woods-Saxon potential is used to determine the single-particle 
energies and wave functions, whereas the RPA residual interaction and terms that mix 
complex configurations are derived from a G-matrix. A more complete SRPA description 
of collective excitations necessitates a fully self-consistent implementation of the effective 
interaction, both in the ground-state and in the construction of the SRPA matrix, and 
also the inclusion of contributions from complex configurations that have been omitted 
in the current versions of the SRPA. 

Another theoretical framework which takes into account the spreading width is 
based on the concept of particle-vibration coupling. Since the nucleus is a highly 
correlated system, a nucleon which propagates in the nuclear medium can excite the 
whole nucleus. The low-lying nuclear excitations correspond predominantly to surface 
vibrations, and theories which are based on particle- vibration coupling take into account 
the fact that the nucleons are confined inside the nucleus because of the surface, but at 
the same time they strongly couple to the dynamical fluctuations of the surface. 

The particle-vibration coupling concept lies at the basis of Nuclear Field Theory 
(NFT) [44]. At the lowest order of NFT the nucleons and the collective vibrations are 
taken as independent degrees of freedom, and nuclear dynamics is determined from the 
hierarchy of their couplings. These couplings are described by a (non-relativistic) field 
theory. Diagrams which correct for the violation of the Pauli principle and appear in 
lowest order, because vibrations are microscopically built from p — h pairs, are also 
included. In the NFT framework the natural extension of the RPA is a model in which 
the lp — lh configurations are coupled with states composed of lp — lh pair plus a 
phonon (instead of 2p — 2h states). This coupling is expressed by the sum of the four 
diagrams depicted in Fig. [21 

The more recent Extended Theory of Finite Fermion Systems (ETFFS) [4"5| 146] 
takes into account lp — lh, complex lp — l/i(g>phonon configurations, the single-particle 
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Figure 3. Diagrams included in the model of Ref. 45J. The upper row contains 
the same diagrams as Fig. [2j whereas the lower row displays the "ground-state 
correlations" diagrams. 

continuum and ground-state correlations. The starting point is the exact Bethe-Salpeter 
equation for the p — h Green's function, which is approximated in such a way that it 
basically corresponds to the RPA solution, plus the contribution of the coupling of 
lp — lh configurations with phonon configurations. These couplings are depicted by 
the diagrams in Fig. [3j The ETFFS approach includes diagrams associated with the 
presence of lp — lh plus phonon components in the ground-state (bottom row of Fig. [3]), 
i.e. "ground-state correlations". 

Another theoretical framework based on the idea that the interaction between 
nucleons and phonons determines nuclear dynamics, is the Quasiparticle-Phonon Model 
(QPM) [4"7] 148] . Within the QPM the excited states of even-even nuclei are treated as 
phonons built from two-quasiparticle pairs, i.e. they correspond to solutions of QRPA 
equations. The ground state is considered as a phonon vacuum and the QRPA yields 
several collective phonons of each multipolarity, as well as many non- collective (or almost 
pure two-quasiparticle) solutions; for simplicity all these solutions are referred to as 
phonons. The main advantage of the QPM approach is that it accounts for the coupling 
between simple (one-phonon) and complex (two- or three-phonon) configurations, even 
in cases which are numerically very demanding. The wave function of the excited 
states represents a combination of one-, two- and three-phonon configurations, with 
the one-phonon configurations corresponding to the set of all QRPA solutions (of given 
multipolarity J n ). The two- and three-phonon configurations are composed of the 
phonons of different multipolarities coupled to the given J 71 . The model Hamiltonian 
is diagonalized in this basis, and the result are the eigenvalues and the microscopic 
structure of each excited state (i.e., the amplitudes of one-, two- and three-phonon 
configurations in the wave function). The NFT and QPM produce similar results. In 
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Figure 4. Basic diagram which gives rise to the spreading width of one-phonon states 
in the QPM. 

the QPM the basic process which gives rise to the spreading width of giant resonances, 
is the coupling of one- and two-phonon configurations. The leading diagram is depicted 
in Fig. HI It has been shown that the largest contributions come from configurations 
in which one of the two phonons in the intermediate state is non-collective. In this 
case the equivalence with the NFT diagrams in Fig. [3] can be demonstrated [39] . The 
model Hamiltonian contains terms which correspond to the mean-field for protons and 
neutrons (with some simplified form of pairing in the case of open-shell nuclei), plus 
a residual effective nucleon-nucleon interaction. In most applications the mean-field 
is a phenomenological Woods-Saxon potential U, and the residual nucleon-nucleon 
interaction has a separable form. A separable interaction is characterized either by 
a form factor r L , or by the Bohr-Mottelson form factor and is used both in the 
(Q)RPA and as the particle-vibration coupling interaction: 

c[c k >, (76) 

k,k';n,LM v ~" 1 " 

where k and k' (n, A and /i) label single-nucleon states (phonons). Rq is the nuclear 
radius, and j3 n L is proportional to the square root of the reduced transition probability 
Eq. (El]), i.e. it is a measure of the collectivity of the phonon state. 

The time-dependent density-matrix (TDDM) model [19] is an extension of the time- 
dependent Hartree-Fock theory beyond the mean-field level. The model describes the 
time evolution of one-body and two-body density matrices and, therefore, includes the 
effects of both a mean-field potential and two-body correlations. TDDM has originally 
been formulated to describe large-amplitude collective motion, but it has also been 
applied in studies of small-amplitude oscillations, in particular low-energy excitations 
in unstable oxygen isotopes [3D] . 

Concluding, we emphasize that the basic advantage of mean-field models lies in the 
fact that they can be formulated in a fully self-consistent way, and easily applied to 
nuclei all over the periodic table. They reproduce global properties of nuclear collective 
excitations, but fail to describe specific phenomena in which particular orbitals around 
the Fermi surface play a special role, or when the inclusion of damping mechanism 
becomes essential. The models that we have described above extend the physical 
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picture beyond the mean-field approximation. They can be considered microscopic 
in the sense that they are based on rigorous formal schemes but, on the other hand, in 
most applications the particle-vibration coupling vertex is phenomenological. 

3.2. Illustrative calculations 

There are, of course, many applications of the NFT theory [51] . We would like to 
mention, in particular, a recent study of the effective pairing interaction which arises 
from the exchange of vibrations between nucleons close to the Fermi surface [52]. Here 
we review another NFT-based approach to nuclear collective excitations, which uses 
Skyrme effective interactions |53j. Starting from the self-consistent Skyrme HF-RPA, 
this model includes the couplings associated with the diagrams shown in Fig. El Both the 
phonons in the intermediate states, and the particle-vibration couplings are consistently 
calculated using the same Skyrme force. 

In the original formulation of the model, which does not include pairing correlations, 
the coupling with the one-particle continuum is considered. In this way both the escape 
and the spreading width of the vibrations are taken into account. In some cases [M] it 
has been possible to describe the particle decay of giant resonances (i.e. partial escape 
widths to definite hole channels in the residual A — 1 system). More recently this 
model has been extended to include pairing, but without considering the coupling to 
the continuum. As an illustrative example, in Fig. [5] we display the calculated IVGDR 
strength distribution in 120 Sn, obtained with the Skyrme force SIII. Both the QRPA 
discrete peaks (vertical bars), and the result obtained with the inclusion of phonon 
coupling - (Q)RPA-PC, are shown. We notice that the dipole strength is significantly 
redistributed by the phonon coupling, but the position of the centroid energy E does 
not change from the QRPA value. The (Q)RPA-PC result is in very good agreement 
with the experimental photoabsorbtion cross section for 120 Sn, which is also included 
in the figure. The most important QRPA-PC additional contributions originate from 
the coupling to the low-energy density vibrations with J n equal to 2 + , 3~, 4 + and 5~. 
When these phonons are collective, corrections associated with the violation of the Pauli 
principle are less important. In addition, the coupling to giant resonances can also be 
taken into account, but this gives minor contributions. The reason can be understood 
from the diagrams shown in Fig. [2j the coupling to low-energy phonons is associated 
with smaller energy denominators. 

In the QPM [18] the model Hamiltonian is composed of terms corresponding to 
the mean-fields for protons and neutrons, a monopole pairing, and a residual nucleon- 
nucleon interaction. The mean field is a phenomenological Woods-Saxon potential U. 
The residual nucleon-nucleon interaction is separable, and can be characterized either 
by a form factor r A , or by a Bohr-Mottelson form factor. The strength parameters of 
the residual interaction are adjusted in each particular nucleus to reproduce the energy 
position and the B(EX) value of the low- lying 2f and 3[" levels. The spurious 1~ state is 
excluded from the excitation spectra by adjusting the isoscalar strength of the residual 
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Figure 5. Photoabsorbtion cross section for 120 Sn, calculated with the QRPA (vertical 
bars) and QRPA-PC (solid curve). The theoretical results are shown in comparison 
with experimental values. 
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Figure 6. Fragmentation of the low- lying electric dipole strength in 138 Ba. 
Calculations are performed in the one-phonon approximation (top panel), and taking 
into account the coupling to two-phonon configurations (middle panel), or to two- and 
thrcc-phonon configurations (bottom panel). 



interaction for A 71 " = 1 , so that this state has zero energy. 

Here we present an example of the fragmentation of one-phonon configurations in 
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QPM calculations [52] • In Fig. E]the low-energy portion of the electric dipole transition 
strength of 138 Ba is plotted for three QPM calculations. In the upper panel the B(E1) 
strength distribution corresponds to the one-phonon approximation, and only five QRPA 
states carry visible El strength. When the coupling to two-phonon configurations 
is included, the El strength becomes strongly fragmented, and is distributed over a 
hundred of states (middle panel). The fragmentation becomes even stronger when three- 
phonon configurations are included in the wave function (bottom panel). The i?(£'l)'s 
of the strongest 1~ states in the bottom part of Fig. [6], are in quantitative agreement 
with the experimental values [56j. In the calculation presented in the bottom panel 
all two- and three-phonon configurations with an energy below 8.5 MeV are included. 
These configurations include phonons with multipolarity from 1* to 9^, and their total 
number is about 1200. At higher energies the density of complex configurations increases 
rapidly, and for a feasible calculation truncation become necessary. 

4. Low-Energy Electric Dipole Strength 

4-1. Low-Energy Response of Light Nuclei 

Evidence of unusually strong dipole response at low-energy in light nuclei was first 
reported in fragmentation experiments of halo nuclei on heavy targets with a large 
number of protons. In the reaction 208 Pb( 11 Li, 9 Li) at 800 MeV/nucleon [57J a large 
Coulomb excitation cross section of 0.9 b was extracted (see also [58]). This cross 
section is associated with a large peak in the B(E1) distribution that appears around 
the two-neutron separation energy, which is of the order of 300 keV. Such a small value 
of the separation energy could be correlated with the enhancement of the low-lying 
dipole strength. Namely, the wave functions of the two weakly bound neutrons which 
form the halo structure are extended far beyond the 9 Li core. As a result, the excitation 
of the halo neutrons to the continuum is intensified, and the low-lying dipole response 
is decoupled from the IV GDR. Initial calculations could reproduce this phenomenon 
only qualitatively. In Ref. |59j for instance, a simple RPA calculation was performed on 
top of a Woods-Saxon mean-field potential, whose parameters were adjusted in such a 
way that the pi/2 neutron state was located at —0.2 MeV. Such an approximation is not 
satisfactory, of course, because recent experiments have shown that, in addition to the 
(P1/2) 2 orbital, the two halo neutrons outside the 9 Li core distribute their spectroscopic 
amplitudes between the (si/2) 2 and (g? 5 / 2 ) 2 configurations [60] . 

The enhancement of the low-lying multipole strength is expected to be a general 
phenomenon in nuclear systems characterized by small values of particle (e.g. neutrons) 
separation energies. For a simple spherical potential well, the lp—lh transition strength 
can be evaluated analytically, and it has been shown that weakly-bound particles 
produce a non-resonant concentration of strength just at the threshold energy [61J. 
The origin of this strength is in the optimal matching of the continuum wave functions 
with the tail of the wave functions of the outermost weakly-bound orbitals (orbitals are 
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said to be weakly bound when the corresponding value of k = ^2m\e\/h, associated 
to the single-particle energy e, is much smaller than the inverse of the nuclear radius). 
Such a simple model could explain the experimentally observed "threshold effect" in the 
response of one-neutron halo nuclei, like n Be |62j. The occurrence of threshold strength 
is also predicted by more realistic continuum RPA calculations. For instance, RPA 
transition densities and currents in n Be have been studied in Ref. [63]. In Fig. [7] we 
plot the strength distributions and transition densities for the lowest isoscalar monopole, 
isovector dipole and isoscalar quadrupole modes. The transition densities differ from 
the prediction of the macroscopic model, not only because of the long tail associated 
with the excitations of the weakly-bound neutrons, but also due to the node of the 
2si/ 2 radial function. Macroscopic models, by definition, do not include single-particle 
shell structures. For n Li, 11 Be, and other light nuclei that we mention below, the term 
"threshold effect" probably applies much better to the phenomenon under study than 
"soft mode". The latter expression is generic, and has also been used for low-lying 
states in stable nuclei, which absorb only a tiny fraction of the total strength for a 
given multipole. Even though the occurrence of pronounced threshold strength can 
be expected also for other multipoles when the separation energies become small, in 
Ref. [M] it has been pointed out that the contributions of L ^ 1 to the cross section are 
negligible in experiments with Coulomb excitations on high-Z targets. 

An interesting question is whether the response of a two-neutron halo system 
like 11 Li, is very different from that of a one-neutron halo, e.g. 11 Be. Two-neutron 
correlations are, of course, essential for the binding of n Li, considering that 10 Li is not 
a bound system. In Ref. [65j the ground state and the dipole response of 11 Li were 
studied by solving a Hamiltonian that includes a Woods-Saxon potential for the two 
neutrons outside the structureless 9 Li core, plus a density- dependent zero-range pairing 
force. With a careful treatment of the continuum [66], it was shown that in this model 
the dipole response of the correlated system differs from the free response by 15%-20%. 
In the description of the dipole response the recoil of the 9 Li was not taken into account. 
This term has been shown to play a role, however, in the quantitative description of 
the ground state of n Li [57] . More generally, one expects that not only the recoil, but 
also the polarization of the 9 Li core is important for a quantitative description of the 
dynamical response [68J. 

Electromagnetic dissociation measurements have been performed for a series of 
neutron halo nuclei: 6 He [69], 8 He [70], n Li pH [72], [73] , 12 Be [74], 14 Be [75], 19 C [76], and 
for the proton halo nucleus 8 B [77J. Appreciable El strength is generally found already 
at low excitation energy, far below the domain of the GDR. For n Li, in particular, the 
El strength observed below 4 MeV excitation energy corresponds to ~ 8% of the TRK 
sum rule, and can be decomposed into at least two broad structures with peak energies at 
w 1.0 MeV and at ~ 2.4 MeV. In a (p,p') experiment, the corresponding peaks have been 
observed at 1.3 and 2.9 MeV [78]. A long-standing issue in n Li is the description of the 
low-energy El excitations in terms of its halo structure and, particularly, the question 
of whether the two halo neutrons are subject to strong correlations and eventually form 
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Figure 7. Left panel: low-energy isoscalar monopole, isovector dipole and 
isoscalar quadrupole strength distributions in 11 Be, calculated with the self-consistent 
continuum RPA |63j . For the peak states the radial transition densities are shown in 
the panel on the right, and compared with collective transition densities. 

a dineutron cluster. While several early experiments [7T1 172] did not find evidence 
for strong correlations between the halo neutrons, a new and significantly improved 
measurement of the low- lying B(E1) distribution [73] has revealed a strong low-energy 
El excitation peaked at E x = 0.6 MeV with B(E1) = 4.5(6) Weisskopf units, which is the 
largest low-lying El strength ever observed in nuclei, and the B(E1) distribution could 
only be reproduced by a three-body model with a pronounced two-neutron correlation. 

The structure and excitations of these light systems are best described in a shell- 
model approach, and calculations have been reported [79J which include model spaces of 
2fruj, or even 3hw configurations. Many more transition amplitudes are included than in 
ordinary mean-field, e.g. RPA calculations, and the coherence between these amplitudes 
is crucial in enhancing the dipole strength at low energy. In order to obtain realistic 
results for the B(E1) values, extended Woods-Saxon single-particle wave functions have 
to be used in calculations, adjusted in such a way that separation energies reproduce 
the experimental values. This procedure introduces spurious components that have to 
be carefully removed from the wave functions. 

Let us also mention that the continuum RPA model with Skyrme interactions 
has been used in studies of the isovector and isoscalar dipole response in 34 Ca, 28 0, 
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Figure 8. Continuum RPA IS and IV dipole strength functions in the nucleus 60 Ca. 
The scale of the IS dipole strength is shown on the right-hand side, whereas on the left 
we plot the scale of the IV dipole strength. The thick curves are obtained by averaging 
the calculated RPA strength with Lorentzian functions of 1 MeV width. 

60 Ca and 22 C [80] • These nuclei could be considered as benchmark systems for the 
dipole response at the nucleon drip-line. It has been shown that the low-lying neutron 
excitations are characterized by the mixing of isovector and isoscalar modes, and that the 
threshold strength is predominantly of isoscalar nature. For the nucleus 60 Ca the dipole 
strength functions calculated with the SkM* force are shown in Fig. [81 The isoscalar 
strength found at high excitation energies above the IV GDR can be associated with 
the dipole compressional mode. Such extremely neutron-rich systems are, of course, not 
yet accessible in experiments, but the calculation nicely illustrates the features of the 
dipole response at the limit of neutron binding. 

4-2. Low-Energy Dipole Excitations in Oxygen Isotopes 

In the oxygen isotopic chain the drip-line nucleus 24 is located only eight neutrons 
away from the stable isotope 16 and, therefore, the evolution of the low-energy dipole 
strength can in principle be traced up to the neutron drip-line. The possible occurrence 
of low-lying collective dipole strength in neutron-rich oxygen isotopes has attracted 
considerable experimental and theoretical interest in recent years. 

While in the stable nucleus 16 earlier experiments [ETJ [82J, [83] did not find 
any appreciable low-lying dipole strength, more recent electromagnetic excitation 
experiments (at beam energies around 600 MeV/nucleon) have confirmed the expected 
occurrence of low-energy dipole strength in the isotopes 18_22 [EU [85]. Low-energy 
transition strength has been observed in all neutron-rich oxygen isotopes that were 
studied, exhausting up to 12% of the TRK sum rule (170"]) at excitation energies below 15 
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MeV, i.e., below the region of giant resonances. The low-energy El structure in 18 and 
20 O was also studied by intermediate-energy Coulomb excitation at 100 MeV/nucleon, 
and new low-energy dipole states were found in 20 O [EU E7J EE]- The dynamics of 
these low-lying dipole excitations, however, has not been resolved in experiments. In 
particular, it is not clear whether some of these states correspond to a collective soft 
mode, or they all simply result from incoherent single-particle excitations. 

Several modern theoretical approaches have recently been employed in the 
description of the evolution of low-lying dipole strength in oxygen isotopes. In one 
of the first studies, large-scale shell-model calculations were performed for neutron-rich 
oxygen isotopes with up to 3ftw excitations in the Op-lsOd-lpOf model space. Pronounced 
low-lying dipole transition strength below 15 MeV was predicted for 17 0, 18 0, 20 O and 
22 0, exhausting rs 10% of the TRK sum-rule [89j[90]. The calculated low-lying strength 
in 17 and 18 O was found to be consistent with the experimental photoreaction cross 
sections. The continuum QRPA on top of the coordinate-space HFB [301 [91], with 
a Woods-Saxon single-particle potential, a Skyrme force as the residual interaction 
in the ph channel, and a density-dependent 5-force in the pairing channel, has also 
been employed to analyse the low-energy modes in oxygen isotopes. Pairing is not 
very strong in these isotopes, in particular the pairing gap is considerably below the 
empirical 12/v^A estimate in 22 > 24 0. However, since the residual pairing interaction in 
QRPA generates dynamical correlation effects on the response function through pair 
density fluctuations, and therefore provides a contribution to the low-lying multipole 
strength, it is important to include a consistent treatment of pairing correlations within 
the HFB+QRPA framework. Moreover, the energy weighted sum rules are fulfilled 
only if the pairing interaction is consistently included both in the solution for the 
stationary ground state, and in the dynamical linear response [301 E7J- In Ref. [37] 
the self-consistent RHB+RQRPA has been applied in the study of multipole excitations 
of neutron-rich oxygen isotopes and, in particular, in the analysis of the evolution of the 
low-lying isovector dipole strength. 

The overall picture emerging from all these calculations is that the onset of dipole 
strength in the low-energy region is caused by nonresonant independent single-particle 
excitations of the last bound neutrons. This is similar to the case of light nuclei discussed 
in the previous subsection. The difference, however, is that for the oxygen isotopes the 
neutron separation energies are larger, i.e. 3.61 MeV for 24 0, and thus the low-lying 
strength is much less pronounced than for the threshold effect in light systems. 

In order to illustrate the evolution of low-lying dipole strength along the chain of 
neutron-rich oxygen isotopes, we show the results of the self-consistent RHB+RQRPA 
calculation [37] , based on the density- dependent effective interaction DD-ME2 [28J plus 
the Gogny D1S force in the pairing channel. The strength distributions associated with 
the dipole operator Eq. (J66J) are displayed in Fig. [9j With the increase of the number 
of neutrons one finds a pronounced fragmentation of the dipole strength, and low-lying 
strength appears below 15 MeV. The centroid energy of the low-lying El states is 
lowered with the increase of neutron excess, whereas the total strength is enhanced. 
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Figure 9. Evolution of the isovector dipole strength distribution in oxygen isotopes, 
calculated within the RHB+RQRPA model using the DD-ME2 effective interaction. 



The electric dipole strength distributions in 18 0, 20 O, and 22 have also been 
analyzed in calculations which go beyond the mean-field level by including the coupling 
of single-quasiparticle states to vibrational modes [92]. By employing the QRPA-PC 
model with up to four-quasiparticle configurations (two uncorrelated quasiparticles plus 
a collective phonon), it has been shown that the calculated total photoabsorption cross 
section below 15 MeV is in very good agreement with experiment. While the simple 
QRPA analyses predict values which are systematically below the data, the coupling 
with phonons increases the cross section in the low-energy region. Because of the 
repulsion between the simple two-quasiparticle states and the complex configurations 
that include a phonon, the former are shifted to lower energy and this increases the total 
QRPA strength in the low-energy region. The QRPA-PC photoabsorption cross sections 
are shown in Fig. [TUJ We note that the calculation predicts the spreading widths, both 
for the low-energy dipole strength and for the giant dipole resonance. 

The QRPA-PC analysis of Ref. [92J is self-consistent, in the sense that the 
only input is the Skyrme force which determines the ground state, and no further 
adjustment of parameters is necessary in the calculation of the response function. A 
more phenomenological model which also emphasizes the role of phonon coupling is the 
quasiparticle representation of the phonon damping model |93j . In Ref. [50] the time- 
dependent density-matrix (TDDM) model, which is an extension of the time-dependent 
Hartree-Fock theory beyond the mean-field level, was used to calculate the isovector 
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Figure 10. Total photoabsorbtion cross section for the isotopes 18 > 20 > 22 O, calculated 
using the full QRPA plus phonon coupling |92| . 



dipole strength functions of the even-A isotopes 18_24 0. By adjusting the strength of 
the residual interaction, the observed isotopic dependence of low-lying dipole strength 
was reproduced. 

In Tabled] the predictions of several models for the low- lying El strength in 18 ~ 24 
are summarized and compared with data [85l 194"] . The sum of the energy- weighted 
El transition strength below 15 MeV is given in units of the classical TRK sum rule. 
Even though all models agree on the overall effect of the neutron excess on the El 
transition strength, significant differences can be noted in isotopes close to the drip- 
line. In particular, the inclusion of particle-vibration coupling brings the results in 
closer agreement with experiment. 
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Table 1. Sum of the energy-weighted dipole strength for 18 24 up to 15 MeV 
excitation energy, in units of the TRK sum rule. 
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Figure 11. The strength function of the IV dipole operator in 22 O (left). The fully 
self-consistent RHB+RQRPA response (solid line) is compared with the RMF+RRPA 
calculation without pairing (dotted line) , and with the RHB+RRPA calculation which 
includes pairing correlations only in the ground state (dashed line). The proton and 
neutron transition densities for the peak at E = 8.54 MeV are shown in the right 
panel. 



The role of dynamical pairing correlations is illustrated in the example of 22 O. 
The RHB+RQRPA isovector dipole transition strength functions are plotted in the 
left panel of Fig. [TT] for three different calculations: a) the RMF+RRPA calculation 
without pairing, b) pairing correlations included in the RHB calculation of the ground 
state, but not in the RQRPA residual interaction (no dynamical pairing), and c) the 
fully self-consistent RHB+RQRPA calculation. The residual pairing interaction in the 
RQRPA generates pronounced dynamical correlation effects on the responses through 
pair density fluctuations. Moreover, the energy-weighted sum rules are only satisfied 
if the pairing interaction is consistently included both in the static RHB and in the 
dynamical linear response. Pairing is, of course, particularly important for the low-lying 
strength. The inclusion of pairing correlations in the full RHB+RQRPA calculation 
enhances the low-energy dipole strength near the threshold [30l [37] . 

For the main peak in the low-energy region 8.54 MeV), in the right panel of 
Fig. [11] we display the proton and neutron transition densities. In contrast to the 
well known radial dependence of the IVGDR transition densities (proton and neutron 
densities oscillate with opposite phases, the amplitude of the isovector transition density 
is much larger than that of the isoscalar component), the proton and neutron transition 
densities for the main low-energy peak are in phase in the nuclear interior, there 
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is no contribution from the protons in the surface region, the isoscalar transition 
density dominates over the isovector one in the interior, and the strong neutron 
transition density displays a long tail in the radial coordinate. However, a detailed 
analysis of R(Q)RPA amplitudes associated with the low-lying states in oxygen isotopes 
indicated that they originate mainly from the single-nucleon transitions from the 
loosely bound neutron orbits [371 195] ■ Similar results have been obtained with Skyrme 
(Q)RPA calculations. For instance, in the study of 28 performed with the self- 
consistent Hartree-Fock-Skyrme plus RPA, it was shown that the strength around 
the threshold originates essentially from uncorrelated excitations of neutrons with 
small binding energies [80J. Even though in open-shell nuclei the number of partially 
occupied configurations increases because of the smearing of the Fermi surface, both 
non-relativistic and relativistic QRPA calculations do not predict the occurrence of 
pronounced collectivity for the low-lying dipole states in neutron-rich oxygen isotopes. 

4-3. Pygmy Dipole Resonances in Heavier Neutron-Rich Nuclei 

Medium-heavy and heavy neutron-rich isotopes are characterized by the appearance of 
a neutron skin, i.e. a layer of excess neutrons on the nuclear surface (961 EH Ell 99J. 
When approaching the neutron drip-line, in particular, the large proton - neutron 
asymmetry leads to a pronounced difference between the corresponding Fermi energies, 
and neutron orbitals just above the Fermi surface can become unbound. The radial 
wave functions of very weakly-bound or unbound neutron states are extended far beyond 
the nuclear surface and this results in the formation of diffuse surface neutron density 
distributions: skin and halo structures. Experimental evidence for the formation of 
neutron skin is available from antiproton absorption [100], heavy-ion reaction cross 
sections |101j . and from studies of isovector dipole and spin-dipole resonances [102] . 
Estimates about the size of the neutron skin can be deduced from experimental radii of 
charge distributions [1031 1104] and mirror displacement energies [105J. 

The question whether the excess neutrons in the skin can be excited to perform 
collective oscillations against the rest of the nucleus, or they only contribute to the non- 
collective threshold strength, has attracted considerable interest in recent years. In the 
former case one expects that, because the outer neutron orbitals are weakly bound, the 
resulting dipole mode will be rather soft, i.e. its excitation energy will be far below the 
giant resonance region. From the theoretical point of view, such a mode also provides 
a unique test of the isospin-dependent components of effective nuclear interactions, 
which are particularly pronounced in nuclei with a large proton - neutron asymmetry. 
Besides being intrinsically interesting as an exotic mode of excitation, the occurrence 
of low-lying dipole strength plays an important role in predictions of neutron capture 
rates in the r-process nucleosynthesis, and consequently in the calculated elemental 
abundance distribution. Namely, although its transition strength is small compared 
to the total dipole strength, the low-lying collective dipole state located close to the 
neutron threshold can significantly enhance the radiative neutron capture cross section 
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on neutron-rich nuclei, as shown in recent large-scale QRPA calculations [1061 1107] . This 
issue will be discussed in more details in Sec. 14.41 

The possible occurrence of a soft dipole mode, or Pygmy Dipole Resonance (PDR) 
in neutron-rich nuclei, has been analyzd using a variety of theoretical approaches. 
Early studies of the PDR were based on rather simple hydrodynamical models which 
involve classical oscillations of the nucleon fluids. These include the three-fluid 
(protons, neutrons in the same orbitals as protons, and excess neutrons) hydrodynamical 
model [108J, the Steinwedel- Jensen |109] and Goldhaber- Teller models [110j . The low- 
lying mode has been qualitatively described as a collective oscillation of the neutron 
enriched surface layer against the core nucleons. It was also suggested, however, that 
the PDR could arise in nuclear systems with only moderate neutron excess, for instance 
in Ca isotopes [111] . 

More recently, microscopic calculations based on Skyrme effective interactions have 
been employed in studies of the isovector dipole response in neutron-rich nuclei: the 
Hartree-Fock+RPA (EH [EE3J Elj, the continuum RPA approaches [801 ESI EEl fTTTl 
11181 1119j . and the self-consistent HFB+QRPA framework formulated in the canonical 
basis [321 1120j . By employing the continuum RPA, rather large escape widths for direct 
neutron decay from low-energy dipole states were estimated, implying a pronounced 
coupling to the continuum [121j . The results of these studies can be summarized as 
follows: (a) the dipole strength distributions in neutron-rich nuclei are more fragmented 
than in stable nuclei; (b) the centroids are calculated at significantly lower energies; (c) 
the ratio of neutron to proton particle-hole amplitudes of low-lying dipole states is much 
higher than in stable nuclei and, accordingly, the isoscalar (IS) transition densities do 
not vanish and isoscalar probes can excite these states. 

The mixing of isoscalar and isovector states in the low-lying dipole response has 
been analyzed in several studies [95l 11121 H22[ 1120] . More information about the 
isospin structure of the PDR can be obtained from a comparison between the RPA 
dipole strength distribution and the unperturbed response. For 132 Sn this is illustrated 
in Fig. [T21 where we display the discrete dipole spectra for the unperturbed Dirac- 
Hartree response and the relativistic RPA response. When the residual interaction is 
turned on, most of the unperturbed strength is pushed towards higher energies, as one 
expects for isovector states. The pygmy states, however, are shifted below the Dirac- 
Hartree response. Since the residual interaction is attractive in the isoscalar channel, 
it appears that the structure of the PDR is predominantly isoscalar. Experimentally 
the isospin structure of the low-lying El states could be, at least in principle, probed 
by a complementary study of (a, a') and {p,p') scattering |123j . Assuming a simplified 
picture, only isoscalar modes should be excited when the scattered a particle is detected 
under extreme forward angles. 

RPA calculations with zero-range Skyrme forces have shown that, when the 
densities of the core nucleons and the excess neutrons are well separated, more like 
in halo nuclei, oscillations of these densities give rise to pronounced low-energy dipole 
strength. On the other hand, if the two densities overlap, which is the case in neutron- 
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Figure 12. The discrete RRPA dipole strength distribution in 132 Sn, in comparison 
with the unperturbed Dirac-Hartree response. 



skin nuclei, then the coupling between the low-lying excitations and the GDR depletes 
the strength of the former [118j . This result has also been confirmed in studies 
which have used the self-consistent Hartree-Fock+RPA with the finite-range Gogny 
interaction [1241 [125] . 

In open-shell neutron-rich nuclei, where pairing correlations play an important role 
also for low-lying excitations, a fully self-consistent QRPA approach is essential. Two 
such frameworks have been developed recently: the relativistic QRPA formulated in the 
canonical basis of the Relativistic Hartree-Bogoliubov model [37J, and the HFB+QRPA 
based on Skyrme energy-density functionals [321 11261 1127] . These models consistently 
employ for the QRPA residual interactions, both in the ph and pp channels, the same 
effective interactions which determine the nuclear ground state. In this way a direct 
relation is established between the unique ground-state properties of exotic nuclei and 
low-lying collective excitations. In addition, the fully self-consistent formulation of 
QRPA is particularly important for excitations in 1~ channel, because it ensures the 
separation of spurious center-of-mass motion without introducing additional adjustable 
parameters. The relativistic (Q)RPA has been employed in several studies of low-lying 
dipole strength in neutron-rich nuclei. These include the analysis of transition densities 
and velocity fields associated with the PDR in 208 Pb [1281 [129] . the evolution of the PDR 
in exotic isotopes far from the valley of /3-stability [95] , the study of the effects of pairing 
correlations on the low-lying El strength in exotic nuclei, the isotopic dependence of 
PDR excitation energies and transition strength distributions [37], and the relationship 
between the PDR excitation energies and one-neutron separation threshold [130] . An 
alternative approach to the RQRPA, based on the response function formalism and the 
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BCS approximation for the description of pairing correlations, has recently been used 
in studies of low-lying El modes in 26,28 Ne [38] and Ni isotopes [131J. 

The fully self-consistent non-relativist ic HFB+QRPA, based on Skyrme density 
functionals and density-independent delta pairing, has recently been employed in an 
extensive analysis of strength functions and transition densities in the J n = + , 1~, 
and 2 + channels for the even Ca, Ni and Sn isotopes from the proton to the neutron 
drip-lines |120j . It has been shown that the low-energy strength increases with neutron 
number in all multipoles. However, in all channels the correlation between strength and 
collectivity is found to be much weaker than in stable nuclei. 

The spreading width of the low-lying dipole transition strength in neutron-rich 
nuclei has been evaluated with the phonon damping model |132j . the consistent 
Skyrme Hartree-Fock + QRPA with phonon coupling [92J, the quasiparticle phonon 
model (QPM) [133[ I134[ [135J, and the Extended Theory of Finite Fermion Systems 
(ETFFS) [4^ | 1136] . These models do not agree on the effects of the coupling to complex 
configurations on low-lying El strength, e.g. for 48 Ca the ETFFS reduces the pygmy 
strength by 31% with respect to the QRPA result [136J, whereas for neutron-rich Sn 
isotopes the QPM predicts a low-lying dipole strength enhanced by a factor 1.5-2 with 
respect to the QRPA value |135j . The differences result from different approaches 
to the coupling with complex phonon configurations. A consistent description of the 
fine structure of low-lying dipole strength would, of course, necessitate a consistent 
implementation of a nuclear effective interaction in the calculation of ground state 
properties, in the (Q)RPA residual interaction, and in the interaction terms which 
describe the coupling to complex configurations. 

On the experimental side, extensive studies of low-lying electric dipole excitations 
have been performed in recent years. Low-lying El states were observed in neutron- 
capture 7-ray spectra |137[ 1138] , and in resonant scattering of real photons [139] . The 
latter method, although mainly restricted to nuclei with moderate proton - neutron 
asymmetry, provides detailed information about the fine structure of dipole transition 
spectra below the neutron threshold. In particular, pronounced low-energy El strength 
has been observed in 44 > 48 Ca [1361 (HQ], 56 Fe and 58 Ni [HE], 88 Sr [H2], 112 Sn [H3], 
U6,i24g n £44^ N=82 isotones [56l CGES EE], and 204,206,207,208 pb [T33I [1471 fT48l fl49] . 
Recent advances in studies of low-lying dipole modes by photon scattering have been 
reviewed in Ref. |150j . Radioactive nuclear beams provide new opportunities for 
studies of low-lying dipole excitations in heavier nuclei with large proton - neutron 
asymmetry [151j . In a recent experiment of the Coulomb dissociation of secondary 
Sn beams produced by in-flight fission of a primary 238 U beam, the dipole strength 
distribution above the one-neutron separation energy was measured in the unstable 
130 Sn and the doubly-magic 132 Sn |152j . In addition to the giant dipole resonance 
(GDR), evidence was reported for a PDR structure at excitation energy around 10 
MeV both in 130 Sn and 132 Sn, exhausting a few percent of the El energy-weighted sum 
rule. 

Similar to the results obtained with nonrelativistic models, the relativistic QRPA 
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Figure 13. The RRPA dipole strength distribution in 132 Sn, calculated with the 
DD-ME2 effective interaction. In the insertions we plot the ground-state proton and 
neutron density profiles, and the proton and neutron transition densities for the peaks 
at 7.8 MeV and 15.3 MeV excitation energy. 



dipole response of neutron-rich nuclei is characterized by the fragmentation of the 
strength distribution and its spreading into the low-energy region. Fully consistent 
R(Q)RPA calculations have shown that with the increase of the number of neutrons 
along an isotopic chain, a relatively strong El peak appears below 10 MeV. The 
dynamics of this peak is very different from that of the isovector giant dipole resonance 
(IV GDR) [37J, EH]. This is illustrated in Fig. [13] with the example of 132 Sn, where we 
plot the RRPA strength distribution which corresponds to the isovector dipole operator, 
and is calculated with the DD-ME2 effective interaction [28J. In the inserted panels we 
display the neutron and proton ground-state density distributions, and the neutron and 
proton transition densities for the low-lying state at 7.8 MeV, and for the IV GDR 
at 15.3 MeV. For the main peak at 15.3 MeV the transition densities display a radial 
dependence which is characteristic for the isovector dipole mode (IV GDR): the proton 
and neutron densities oscillate with opposite phases. The dynamics of the state at 
7.8 MeV is completely different: the proton and neutron transition densities are in 
phase in the bulk of the nucleus, whereas only neutron excitations contribute to the 
transition density in the surface region. Thus the low-lying pygmy state does not belong 
to statistical El excitations sitting on the tail of the IV GDR, but rather represents a 
new mode - the PDR: the neutron skin oscillates against the core. The neutron skin, 
i.e. the difference between the neutron and proton density distributions in the ground 
state (shown in the right panel in Fig. [T3|) basically determines the properties of the 
PDR |134j . Therefore, for a quantitative description of PDR dynamics it is essential 
to use effective interactions that reproduce available data on the neutron skin. This is 
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the case, for instance, of the relativistic density- dependent interactions DD-ME1 [27J 
and DD-ME2 (28], which have been specifically designed to reproduce the differences 
between the rms-radii of neutron and proton density distributions. 

In light nuclei the low-energy dipole strength predominantly originates from non- 
resonant independent single particle excitations of the loosely bound neutrons. However, 
the structure of the low-lying strength changes with mass. As has been shown in the 
RRPA analysis of Ref. [95J, in heavier nuclei some of the low- lying dipole states display a 
more distributed structure of the RRPA amplitudes. Among several peaks characterized 
by single particle transitions, a single collective dipole state is identified below 10 MeV 
and its RRPA amplitude presents a coherent superposition of many neutron particle- 
hole configurations. For instance, in the case of of 132 Sn (see the dipole strength 
distribution in Fig. [TBI the following neutron ph transitions principally contribute to the 
RRPA amplitude of the state at 7.8 MeV: 3s 1/2 -> 3p 3/2 (51%), 2d 3/2 -> 3p 3/2 (19%), 
24/2 -> 3p 1/2 (11%), 3s 1/2 -> 3p 1/2 (7%), lh n/2 -> lz 13/2 (4%), lg 7/2 -> lh 9/2 (0.9%), 
24/2 -> 3p 3 /2 (0.4%), 2d 5/2 -> 2/ r/2 (0.3%), 2d 3/2 - 4p 1/2 (0.2%), l# 7/2 -> 2/ 5/2 (0.1%), 
etc. On the other hand, the total contribution from all proton transitions to the state at 
7.8 MeV is small: ~ 3%, thus the ratio of neutron to proton contribution is much higher 
than the value N/Z, typical for the IV GDR state. Such a rich structure of the RRPA 
amplitude is in contrast to the situation found in light neutron-rich nuclei, where the 
low-lying dipole peaks below 10 MeV are usually dominated by just one or two neutron 
ph transitions. The level of collectivity can be further enhanced in open-shell nuclei, 
where because of pairing correlations many additional neutron states become partially 
occupied and, therefore, many more 2qp transitions contribute to the RRPA amplitude. 

A similar analysis of neutron particle-hole components of strong low-energy 1~ 
excited states in 132 Sn has also been performed in Ref. |120j . for the self-consistent 
HFB+RPA calculation with the Skyrme SkM* interaction. The distribution of the 
largest neutron ph components and the degree of collectivity for the most pronounced 
low-energy states is comparable to the results of the relativistic RPA, but the two 
models differ in the integrated energy-weighted strength in the low-energy region. While 
relativistic RPA calculations typically predict « 5% of the classical TRK sum rule in 
the energy region below 10 MeV, only about 1% is obtained in the calculation with 
the Skyrme SkM* interaction. In Ref. [120j it has been suggested that this difference is 
related to the larger neutron skin typically calculated with relativistic mean-field models, 
which may be responsible for the more pronounced pygmy resonances. However, this 
does not seem to be the case for the DD-ME2 relativistic interaction (see Fig. [T9|) . which 
does not overestimate the empirical values for the neutron skin in Sn isotopes and, at 
the same time, reproduces the experimental results for the integrated energy-weighted 
dipole strength in the low-energy region. 

The RPA and QRPA analyses of the dynamics of low-lying El strength distributions 
have mostly been performed on the mean-field level, i.e. without taking into account the 
spreading effects which arise from the coupling of single-nucleon states to the collective 
low-lying excitations (phonons). The principal effect of the particle- vibration coupling 
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is an increase of the nucleon effective mass at the Fermi surface, and this is reflected 
in an increase of the density of single-nucleon states close to the Fermi energy. It has 
been argued that the inclusion of particle-vibration coupling in (Q)RPA calculations, 
i.e. extending the (Q)RPA model space to include selected two-quasiparticle <g> phonon 
states, would not only improve the agreement between the calculated and empirical 
widths of the GDR structures, but it could also have a pronounced effect on the low-lying 
El strength. For instance, the coupling to low-lying phonons could fragment the PDR 
structure over a wide region of excitation energies. As a result of this fragmentation only 
an enhancement of the El strength would be observed in the low-energy region, rather 
than a prominent PDR peak. The importance of particle-vibration coupling effects 
for the multipole response of neutron-rich nuclei has particularly been emphasized in 
studies that have used the QRPA plus phonon coupling model based on the Hartree- 
Fock (Q)RPA with Skyrme effective forces pi [153] . In Ref. [153] the QRPA plus 
phonon coupling model was applied in the analysis of dipole excitations in 208 Pb, 120 Sn 
and 132 Sn. In contrast to the results obtained in the relativistic (Q)RPA framework, 
the QRPA plus phonon coupling model predicts low-lying El strength of non- collective 
nature in all three nuclei. 

In Fig. [TH we display the photoabsorbtion cross section for 132 Sn, calculated 
with the fully consistent RPA and RPA-PC models using the Skyrme force SIII. The 
corresponding RPA-PC transition densities for the GDR state at 13.5 MeV and for 
the most pronounced low-energy peak at 9.7 MeV are shown in Fig. [15j Even though 
the transition densities, both for the GDR and for the low- lying peak at 9.7 MeV, 
are similar to those calculated with the relativistic RPA (see Fig. [T3|) the analysis of 
the structure of RPA (RPA-PC) amplitudes shows that none of the peaks below 10 
MeV contain contributions of more than two or three different neutron particle-hole 
(ph) configurations. Predominantly these peaks correspond to just a single-neutron 
transition, and each of them exhausts less than 0.5% of the energy-weighted sum rule. 

Low-lying El excitations in neutron-rich Sn isotopes have also been studied in 
the Quasiparticle Phonon Model [134] . in a model space that included up to three- 
phonon configurations built from a basis of QRPA states, and with separable multipole- 
multipole residual interactions. The single-nucleon spectra were calculated for a Woods- 
Saxon potential with adjustable parameters. Empirical couplings were used for the 
QPM residual interactions. In the QPM spectra for 120_132 Sn the low-energy dipole 
strength was found concentrated in a narrow energy interval such that the PDR could 
be identified. It was shown that, despite significant multi-phonon contributions to the 
mean-energy and transition strength, the PDR states basically retain their one-phonon 
character. 

Because of its relatively large neutron excess, the stable nucleus 208 Pb has also been 
investigated for a possible occurrence of pygmy dipole resonant states. Experimental 
evidence has been reported in elastic photon |154j and photoneutron scattering [155j . 
and in electron scattering [156] : pronounced El strength has been observed in the 
energy region between 9 and 11 MeV, several MeV below the IV GDR in 208 Pb. On 
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Figure 14. Photoabsorption cross section for 132 Sn, calculated with the RPA and 
RPA-PC models. The effective interaction is Skyrme SHI. 



the theoretical side, one of the first microscopic analysis was performed in the Hartree- 
Fock plus RPA model based on the Skyrme interaction SGII |157j . Two pronounced 
peaks were calculated at 8.7 MeV and 9.5 MeV, which appeared as likely candidates for 
the PDR. In a recent self-consistent relativistic RPA study based on the NL3 effective 
interaction, a pronounced low-energy dipole peak was calculated at 7.29 MeV [128] . The 
structure of the RRPA amplitude, the corresponding transition densities and velocity 
fields indicate that this state can be interpreted as a collective PDR mode. The RRPA 
prediction for the PDR state has been confirmed in a subsequent (7,7') experiment, 
which disclosed a resonance-like structure centered at 7.37 MeV, approximately at the 
neutron emission threshold [133J. 

In Fig. [16] we display the isovector dipole strength distribution in 208 Pb, evaluated 
with the self-consistent RRPA model employing the DD-ME2 effective interaction [28J. 
The calculated energy of the main peak at 13.4 MeV is in excellent agreement with 
the experimental value of the excitation energy of IV GDR: 13.3 ± 0.1 MeV |158j . The 
pronounced low-energy peak at 7.2 MeV is close to the experimental centroid of the low- 
lying dipole strength at 7.37 MeV [133] . In the panels on the right in Fig. [16] we plot the 
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Figure 15. RPA-PC transition densities for 132 Sn. In the upper panels the proton 
and neutron (right panel), and isoscalar and isovector (left panel) transitions densities 
for the state at 9.7 MeV are shown. For comparison in the lower panels we plot the 
corresponding transition densities for the IV GDR at 13.5 MeV. 



RRPA transition densities for the low-energy state at 7.2 MeV, and the IV GDR state at 
13.4 MeV, respectively. Obviously, the dynamics of the low-lying mode is very different 
from that of the isovector giant resonance: the proton and neutron transition densities 
are in phase in the nuclear interior and there is large contribution from the neutrons in 
the surface region. For the IV GDR state the total isovector transition density is much 
stronger than the isoscalar component. On the other hand, for the state at 7.2 MeV 
the isoscalar transition density dominates over the isovector one in the interior, and the 
large neutron component in the surface region contributes to the formation of the node 
in the isoscalar transition density [951 IH2[ 1122] . 

Low-lying El excitations have also been observed below 10 MeV in (7, 7') scattering 
on the N=82 isotones: 138 Ba, 140 Ce, 142 Nd, and 144 Sm [HB]. [T23"] . The subsequent analysis 
of the RHB+RQRPA transition densities for the calculated low-lying states in these 
nuclei has shown that a collective PDR mode indeed develops. However, the calculated 
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Figure 16. The RRPA dipole strength distribution in 208 Pb, calculated with the 
DD-ME2 effective interaction. The neutron, proton, isoscalar(IS), and isovector(IV) 
transition densities, for the pygmy state at 7.2 MeV and the IV GDR state at 13.4 
MeV are plotted in the panels on the right. 



PDR are located ~ 1-2 MeV above the experimental centroids [37]. The evolution 
of the isovector dipole strength distribution in N=82 isotones, evaluated in the fully 
self-consistent RHB+RQRPA with the DD-ME2 effective interaction, is illustrated in 
Fig. [T71 The dotted vertical line separates the low-energy region below 10 MeV from 
the region of giant resonances. In contrast to the IV GDR, which weakly decreases 
in excitation energy with the increase of the proton number, i.e. with mass number, 
the centroids of the low-lying structure increase in energy, whereas the total low-energy 
strength decreases when the proton - neutron asymmetry is reduced. One notices that 
the low-lying states are far more sensitive to the variations of the proton number, than 
the IV GDR structure. These observations are consistent with the interpretation of the 
low-energy peaks in terms of the PDR, because the reduction of the asymmetry between 
the neutron and proton density distributions in the ground state should generally result 
in higher PDR excitation energies and in the suppression of its strength |130j . The 
RHB+RQRPA B(E1) strength in the low-energy region below 10 MeV decreases with 
mass number along the N=82 isotone chain (lower panel in Fig. [T71) . but the calculated 
values are systematically above the data [56JI123J. On the other hand, the quasiparticle 
phonon model (QPM) predicts a constant summed El strength in all the measured N=82 
nuclei |159j . The reason is that, in contrast to the fully self-consistent RHB+RQRPA 
approach, the QPM calculations employ the same single-particle spectrum for all the 
N=82 isotones |160j . and therefore cannot describe the details in nuclear structure which 
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Figure 17. The isovector dipole strength distributions in N=82 isotones calculated 
in the RHB+RQRPA model, with the DD-ME2 effective interaction (upper panel). 
The integrated B(E1) transition strength below 10 MeV, for the NL3 and DD-ME2 
interactions, is compared with the data from (7,7') scattering [551 1123] flower panel). 

The theoretical analysis of low-lying excitations and the comparison with data, 
can be used as a sensitive tool to constrain the isovector channel of effective nuclear 
interactions. In the lower panel of Fig. [T7|we compare the theoretical B(E1) strength 
in the low-energy region below 10 MeV with data. The theoretical values have 
been calculated in the consistent RHB+RQRPA model with the very popular non- 
linear meson-exchange effective interaction NL3 [21], and with the new meson-exchange 
interaction DD-ME2 [28] , which explicitly includes a medium dependence of the meson- 
nucleon couplings. Obviously the NL3 interaction, which is known to overestimate the 
size of the neutron skin not only in exotic neutron-rich nuclei but also in 208 Pb, predicts 
too much low-lying B(E1) strength. On the other hand, an interaction like DD-ME2 
which has been adjusted to the empirical differences between the radii of neutron and 
proton density distributions, significantly improves the agreement of the calculated low- 
energy dipole strength with data. The remaining difference might be caused by the 
coupling with more complex phonon configurations [135 ], not taken into account in the 
RHB+RQRPA models, or in the missing El strength in (7, 7') scattering which may be 
quite considerable when dealing with end-point energies close to the neutron separation 
threshold |144j . 
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4-4- Isotopic Dependence of Pygmy Dipole Resonances 

The evolution of low-lying dipole transition strength along an isotopic chain provides 
useful information about the underlying dynamics of soft modes in exotic nuclei. In 
particular, an important question is the location of PDR with respect to the neutron 
separation threshold [130J. This is important not only for the possible detection of 
PDR in experiments, but also for modeling the r-process nucleosynthesis [ 1 TJ . Most 
of the recent photon scattering experiments provide data on the dipole strength only 
below the neutron separation energy, i.e. only a portion of the overall low-lying El 
strength is observed [56J [1361 IHH ED E21 El ESI E3 EH EH El] . Of course for 
a more complete understanding of the structure of the low-lying dipole strength and its 
relation to the PDR mode, data on transition strength above the neutron threshold are 
necessary |130j . 



12 



10 



1 


i 1 i 1 i 1 i 

RHB+RQRPA 




1 1 

Sn 




4 


(DD-ME2) 












































OO E L0 w" Centr0id 










A - A E Low _ p eak ener sy 

<X>S n -RHB 








X A - 




♦♦S n -EXP 












■ S -EXTR. 

™ n 




" o 

1 


o - 


, 


1,1,1 


1 





> 



o 
-1 



112 



116 



120 



124 128 
A 



132 



136 



140 



Figure 18. The calculated PDR peak and centroid energies, and the one-neutron 
separation energies for the sequence of Sn isotopes, as functions of the mass 
number. The DD-ME2 effective interaction has been used in the RHB+RQRPA 
calculations. The RHB results for the neutron separation energies are compared with 
the experimental and extrapolated values [161j . 

This is illustrated in Fig. [18j where we display the RHB+RQRPA results for the 
peak and centroid energies of the PDR in a series of Sn isotopes. The RQRPA predicts 
a monotonic decrease of the PDR with mass number, and only a small kink in the 
peak excitation energies is calculated at the N = 82 shell closure. In the same plot 
we have also included the calculated one-neutron separation energies, in comparison 
with the data and the extrapolated value |161j . The self-consistent RHB calculation, 
with the DD-ME2 mean-field effective interaction in the ph channel and the D1S Gogny 
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force in the pairing channel, reproduces in detail the one-neutron separation energies in 
Sn nuclei. We notice that the separation energies decrease faster than the calculated 
PDR excitation energies. At the doubly closed-shell nucleus 132 Sn a sharp reduction of 
the one-neutron separation energy is observed and reproduced by the RHB calculation, 
whereas the shell closure produces only a much weaker effect on the PDR peak energies. 
The increased fragmentation of the low-lying strength in heavier Sn isotopes results in 
larger differences between the PDR peak and centroid energies. The important result 
here is that for A < 122 the PDR excitation energies are below the corresponding 
one-neutron separation energies, whereas for A > 122 the pygmy resonance is located 
above the neutron emission threshold. This means, of course, that in the latter case the 
observation of the PDR in (7,7') experiments will be strongly hindered [1301 H44] . 




Figure 19. The RHB+RQRPA energy-weighted dipole strength, integrated up to the 
energy cut-off E c =10, 11, and 12 MeV, respectively, and plotted in percents of the 
TRK sum rule. The experimental results are from Refs. [1431 H44] I152j . 

The presently missing data on dipole strength above the neutron threshold 
could be obtained in the near future by using tagged photons at S-DALINAC [123j . 
In addition, photon scattering with high intensity beams at energies below and 
above the neutron separation threshold are planned at the superconducting electron 
accelerator ELBE [1621 1163j . The first studies at ELBE include photon-scattering on 
92,98,100 Mo [164J. It is interesting to note that the data show an enhancement of the 
dipole transition strength around 9 MeV: in 92 Mo the pygmy strength is located below 
the neutron separation energy, whereas in 100 Mo it shifts above the neutron threshold. 

In Fig. [19] we display the isotopic dependence ( 112 Sn- 140 Sn) of the energy weighted 
dipole strength in the low-energy region, integrated up to the cut-off energy E c =10, 
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11, and 12 MeV, respectively, and plotted in percents of the classical TRK sum rule. 
Model calculations are performed in the RHB+RQRPA with the DD-ME2 plus Gogny 
D1S interactions, and the results are compared with the available data from photon 
scattering [143[ 1144] , and Coulomb dissociation of secondary Sn beams from in-flight 
fission [152J. The calculated low-lying El strength is in excellent agreement with the 
recent experimental data for 112 Sn |143j and 130 > 132 Sn [152J, whereas it overestimates 
the (7,7') data for 116,124 Sn [144] . When considering the evolution of low-lying dipole 
strength along an isotopic chain, in a first approximation one could expect that the 
relative strength of the PDR increases monotonically with the number of neutrons, at 
least within a major shell. In the case of Sn isotopes the RHB+RQRPA calculations 
predict, however, that the PDR peak is most pronounced around 124 Sn (depending on the 
cut-off, see Fig. [19]) [37] . A combination of shell effects and reduced pairing correlations, 
leads to a reduction of the strength of the PDR in heavier Sn nuclei below N = 82. 
The local minimum in the low-lying El strength is calculated for 132 Sn, whereas in the 
neighboring isotopes the transition strength increases because of enhanced collectivity, 
i.e. the increase in the number of two-quasiparticle pairs contributing to the RQRPA 
configuration space. We also notice the pronounced difference in the pygmy strength 
between nuclei close to the valley of /^-stability and exotic nuclei: while below the N=82 
shell closure the integrated transition strength is at most ~ 4% of the TRK sum rule 
value (for E c =10 MeV), beyond 132 Sn the PDR strength exhibits a strong enhancement. 




Figure 20. The RHB+RQRPA calculated PDR peak energies for Sn, Te, and Xe 

isotopes. 



It is, of course, interesting to explore other isotopic chains of spherical nuclei where 
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Figure 21. Same as Fig. [T9l but for the Sn, Te, and Xe isotopes. The cut-off energy 
for the low-lying El strength is set at £' c =10 MeV. 

one expects the occurrence of the PDR in the El excitation spectrum. In Figs. [201 and [211 
we plot the calculated PDR peak energies and the integrated low-lying El strength for 
Sn, Te, and Xe isotopes. The cut-off energy for the low-lying strength is arbitrarily 
set at E c =10 MeV. The calculated low- lying dipole strength appears rather sensitive to 
small variations in the number of nucleons. The PDR excitation energies are lowest in 
the isotopic chain with the smallest number of protons, i.e. in Sn nuclei. This behavior 
reflects the nature of the PDR: a larger neutron excess should result in lower PDR 
excitation energy. In the region beyond A=132, the slope of the PDR peak energies 
becomes steeper than for stable nuclei, because the neutrons in outer orbitals are more 
loosely bound and thus the restoring force in the oscillation of the skin against the core 
becomes weaker. For Te and Xe isotopes we plot the PDR peak energies starting from 
118 Te and 120 Xe, respectively. In the lighter systems the PDR could not be uniquely 
identified. The energy-weighted dipole transition strength in the region below 10 MeV 
(Fig. I2T1) is strongest for the Sn chain, and somewhat weaker for Te and Xe. This is, of 
course, to be expected because the PDR strength must be proportional to the neutron 
excess. In all the three chains the local minima in the integrated transition strength are 
obtained at N=82, and the PDR strength rapidly increases beyond the neutron shell 
closure. 

Finally, in Fig. El the RHB+RQRPA results for the PDR in Ni and Pb isotopes are 
shown. The RHB neutron separation energies, calculated with the DD-ME2 plus Gogny 
D1S interactions, are compared with the experimental values |161j . In the sequence of Ni 
nuclei the crossing between the theoretical curve of one-neutron separation energies and 



Exotic modes of excitation in atomic nuclei far from stability 



54 




Figure 22. Same as in Fig. [TH but for the Ni and Pb isotopic chains. The large 
square denotes the experimental position of the PDR in 208 Pb [133] . 



the PDR excitation energies is calculated already at A = 66. In heavier Ni isotopes the 
excitation energy of the PDR is predicted high above the neutron emission threshold. 
One should notice that for the lighter Ni isotopes the agreement between the calculated 
and experimental neutron separation energies is not as good as for the Sn nuclei and, 
therefore, the actual point of crossing between the PDR and the one-neutron separation 
energy could occur for A < 66. The Ni nuclei are not very rigid and, for a more 
quantitative description, one would have to go beyond the simple mean-field plus QRPA 
calculation and include correlation effects. For the Pb isotopes the crossing point is 
calculated at A = 208, in excellent agreement with the data on low-lying El excitations 
in 208 Pb p3j. 

Motivated by the experimental results on the PDR in 130 Sn and 132 Sn [152], the 
relativistic RPA has recently been applied in the study of the isotopic dependence of the 
PDR in tin |165j . focused on the following questions: (a) is there a correlation between 
the development of a neutron skin and the emergence of low-energy dipole strength? 
and (b) can the data be used to discriminate among effective interactions that predict 
different values for the neutron skin in heavy nuclei? The results of the RPA analysis 
are not conclusive because, although a strong linear correlation between the neutron 
skin and the fraction of the energy-weighted sum rule at low energy was observed, an 
anti-correlation actually developed beyond 120 Sn, and it was attributed to the filling 
of the neutron l/tn/2 orbital. It should be pointed out, however, that the analysis 
was performed on the RPA level, without considering the effect of pairing correlations. 
Comparing different effective interactions, it was found that the centroid energy of the 
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PDR is not sensitive to the density dependence of the symmetry energy. The fraction 
of the energy-weighted sum rule exhausted by the PDR, on the other hand, increases 
sharply with increasing neutron skin. Although the experimental error bars are large, 
the data seem to disfavor effective interactions with stiff symmetry energy, i.e. those 
which predict excessively large neutron skins. 

4-5. The Proton Electric Pygmy Dipole Resonance 

Because the proton drip-line is much closer to the line of /^-stability than the neutron 
drip-line, bound nuclei with an excess of protons over neutrons can be only found in the 
region of light Z < 20 and medium mass 20 < Z < 50 elements. For Z > 50, nuclei 
in the region of the proton drip-line are neutron-deficient rather than proton-rich. In 
contrast to the evolution of the neutron skin in neutron-rich systems, because of the 
presence of the Coulomb barrier, nuclei close to the proton drip-line generally do not 
exhibit a pronounced proton skin, except for very light elements. Since in light nuclei 
the multipole response is generally less collective, all these effects seem to preclude the 
formation of the pygmy dipole states in nuclei close to the proton drip-line. Nevertheless, 
a recent analysis based on the RHB+RQRPA approach has shown that proton pygmy 
dipole states can develop in light and medium mass proton-rich nuclei [166J. 

In Fig. [23] we plot the ground state density profiles for the Ar isotopes and for the 
N=20 isotones, respectively, evaluated in the RHB model using the DD-ME2 effective 
interaction, and the Gogny D1S force in the pairing channel. In the both examples 
we observe the formation of the proton skin on the surface of nuclei which have a 
higher ratio of protons over neutrons. The proton skin is, of course, not so pronounced 
as neutron skin in neutron-rich nuclei, because of the Coulomb barrier which tends 
to localize the protons in the nuclear interior. Evidence for a possible formation of 
the proton skin in neutron-deficient or proton-rich nuclei has been reported in recent 
experimental studies |167j . and is supported by model predictions [168[ 1169J. 

Only few studies of dipole excitations in proton-rich nuclei have been reported 
so far. The isovector dipole response of the proton drip-line nucleus 34 Ca has been 
analyzed with the continuum RPA based on Skyrme interactions, and a multiple peak 
structure has been predicted between the low-energy isoscalar dipole response and the 
IV GDR [80]. In the large-scale shell-model calculations for 13 [170] , pronounced El 
strength has been found in the low-energy region below 3 MeV, and related to the 
coherence in the transition amplitudes between the loosely-bound valence nucleons, and 
also between the core and the valence nucleons. The relativistic RPA calculation of the 
dipole response in Ar isotopes [171] , has shown a concentration of strength in the proton- 
rich nuclei 30 Ar and 32 Ar, which has been attributed to excitations from weakly-bound 
single-particle states into the continuum. A recent study based on the self-consistent 
RHB+RQRPA framework has shown that a new collective mode - the Proton Pygmy 
Dipole Resonance (PPDR) could arise in medium-heavy nuclei close to the proton drip- 
line [166] . In Fig. [2H the RQRPA dipole strength distributions in the N=20 isotones 
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Figure 23. Neutron and proton ground-state density profiles for Ar isotopes (left 
panel), and N=20 isotones (right panel), evaluated in the RHB model with the DD- 
ME2 effective interaction. The arrows denote the increasing (decreasing) proton 
(neutron) density distributions along the isotopic, and isotonic chains. 

40 Ca, 42 Ti, 44 Cr, and 46 Fe are shown, calculated with the fully consistent RHB+RQRPA 
with the DD-ME2 plus Gogny D1S effective interactions. The strength distributions are 
dominated by the IV GDR at ~ 20 MeV excitation energy. With the increase of the 
number of protons, low-lying dipole strength appears in the region below the GDR 
and, for 44 Cr and 46 Fe, a pronounced low-energy peak is found at ~ 10 MeV excitation 
energy. In the lower panel of Fig. [24] we plot the proton and neutron transition densities 
for the peaks at 9.98 MeV in 44 Cr and 9.33 MeV in 46 Fe, and compare them with the 
transition densities of the GDR state at 18.82 MeV in 46 Fe. Obviously the dynamics of 
the two low-energy peaks is very different from that of the isovector GDR: the proton 
and neutron transition densities are in phase in the nuclear interior and there is almost 
no contribution from the neutrons in the surface region. As in the case of the PDR 
in neutron-rich nuclei, obviously the low-lying state does not belong to statistical El 
excitations sitting on the tail of the GDR, but could indeed represent a fundamental 
mode of excitation: the proton electric pygmy dipole resonance (PPDR). 

In Fig. [25] we analyse the RQRPA structure of the dipole response in 46 Fe. This 
nucleus is located at the proton drip-line, and recently evidence for ground-state two- 
proton radioactivity was reported in the decay of 45 Fe |172[ 1173] . In the four panels we 
plot the QRPA amplitudes of proton and neutron 2qp configurations 

£,2qp = \^2qp\ ~~ \^2qp\ (77) 

for the three low-lying states at 8.91, 9.33, and 10.11 MeV, as well as for the strongest 
state in the GDR region at 18.82 MeV. For each of the four dipole states, in addition 



Exotic modes of excitation in atomic nuclei far from stability 



57 




Figure 24. The RHB+RQRPA isovector dipole strength distributions in the N=20 
isotones (upper panel). The proton and neutron transition densities for the low-lying 
states in 44 Cr and 46 Fe, and for the IV GDR state in 46 Fe are shown in the lower panel. 
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Figure 25. The distributions of the RQRPA amplitudes in 46 Fe, plotted as functions 
of the unperturbed energy of the respective 2<7p-configurations, for the main peaks in 
the low-energy region and for the IV GDR state. 
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to the excitation energy we have also included the corresponding B(E1) value. The 
amplitudes are shown in a logarithmic plot as functions of the unperturbed energy of 
the respective 2gp-configurations. Only amplitudes which contribute more than 0.01% 
are shown, and we also differentiate between proton and neutron configurations. We note 
that, rather than a single proton 2qp excitation, the low-lying states are characterized 
by a superposition of a number of 2qp configurations. Obviously the pygmy states 
display a degree of collectivity that can be directly compared with the QRPA structure 
of the GDR state at 18.82 MeV. In addition, proton 2qp configurations account for 
~99% of the QRPA amplitude of the pygmy states, whereas the ratio of the proton to 
neutron contribution to the GDR state is ~ 2. For the GDR states the 2qp configurations 
predominantly correspond to excitations from the sci-shell to the /p-shell. The structure 
of the pygmy states, on the other hand, is dominated by transitions from the I/7/2 proton 
state at -0.21 MeV, and from the 2p 3 / 2 proton state at 3.63 MeV (this state is only bound 
because of the Coulomb barrier). The energy weighted sum of the strength below 11 
MeV excitation energy corresponds to 2.7% of the TRK sum rule. 

Another example where a pronounced proton PDR can occur are the proton-rich 
isotopes of Ar. In the left panel of Fig. [26] we display the RHB+RQRPA electric dipole 
strength distribution in 32 Ar. In addition to the rather fragmented GDR structure at 
pa 20 MeV, prominent proton PDR peaks are calculated at 8.14, 8.79, 9.22, and 9.46 
MeV. These peaks form the pygmy structure and exhaust 5.7 % of the TRK sum rule. 
The RQRPA amplitudes of the low-lying states present superpositions of many proton 
2qp configurations, with the neutron contributions at the level of 1%. The dominant 
configurations correspond to transitions from the proton states Ids/2 (-1.94 MeV) and 
2si/2 (-3.98 MeV). In the right panel of Fig. [26] we display the mass dependence of 
the centroid energy of the pygmy peaks and the corresponding values of the integrated 
B(E1) strength below 10 MeV excitation energy. In contrast to the case of medium- 
heavy and heavy neutron-rich isotopes, in which both the PDR and GDR are lowered 
in energy with the increase of the neutron number, in proton-rich isotopes the mass 
dependence of the PDR excitation energy and B(E1) strength is opposite to that of the 
GDR. The proton PDR decreases in energy with the development of the proton excess. 
This mass dependence is intuitively expected because the proton PDR is dominated by 
transitions from weakly-bound proton orbitals. As the proton drip-line is approached, 
either by increasing the number of protons or by decreasing the number of neutrons, 
due to the weaker binding of higher proton orbitals one expects more inert oscillations, 
i.e. lower excitation energies. The number of 2qp configurations which include weakly- 
bound proton orbitals increases towards the drip-line, resulting in an enhancement of 
the low-lying B(E1) strength. 

For heavier nuclei the proton drip-line is located in the region of neutron-deficient, 
rather than proton-rich nuclei, and therefore one does not expect to find low-lying dipole 
strength in medium-heavy and heavy nuclei close to the proton drip-line. 

The effect of the coupling to the continuum on the low-lying dipole strength in nuclei 
close the proton drip-line has also been analyzed in the non-relativistic continuum RPA 
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Figure 26. The RHB+RQRPA isovector dipole strength distribution in 32 Ar (left 
panel). The mass dependence of the PPDR centroid energy for Ar isotopes, and the 
corresponding values of the integrated B(E1) strength below E c =10, 11, and 12 MeV 
are shown in the right panel. 



(CRPA) framework, using Skyrme interactions [80l I174j . In the doubly magic nucleus 
48 Ni both the CRPA and the RRPA predict the occurrence of the proton PDR, but the 
CRPA analysis has shown that, as a result of the coupling to the continuum, the PPDR 
is characterized by a rather large escape width |174j . 



4-6. Di-neutron Correlations near the Drip-Line 

The existence of another exotic mode has been suggested in medium-heavy nuclei close 
to the neutron drip-line: a soft dipole excitation that corresponds to the vibration of 
a di-neutron in the nuclear exterior against the remaining A-2 subsystem |1751 I176j . 
While in nuclei with a pronounced neutron excess the pygmy dipole resonance could 
appear, di-neutron vs core vibrations may occur in very exotic nuclei near the neutron 
drip-line. The latter mode is strongly influenced by neutron pairing correlations, and is 
characterized by a large transition density for pair motion of neutrons. In the case of 
light halo-nuclei pairing correlations between the loosely-bound neutrons in the halo lead 
to a strong enhancement of the soft dipole excitations. Experimental signatures of di- 
neutron correlation in the soft dipole mode have been found in n Li [721173] . However, as 
has recently been shown in Ref . [175] , some features of di-neutron correlations may also 
be present in the ground states of exotic medium-mass nuclei, and therefore influence 
their excitations. 
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In the HFB description of the nuclear ground state (cf. Sec. 12. ip . the spatial 
correlations between a pair of neutrons can be probed by the neutron two-body 
correlation density [175j . 

Pcorr (ra,r'a') = \ K (ra,r'a')\ 2 - \p(ra, rV)| 2 , (78) 

defined by the off-diagonal terms of the pairing tensor (Eq. (151)), and density matrix 
(Eq. (J2J)). The notation a has been defined in Eq. (149p . In the case of medium- 
heavy nuclei close to the neutron drip-line, it has been shown that the two-body 
correlation density clearly reflects the presence of spatial di-neutron correlations in the 
pair-correlated ground state, and these are especially pronounced on the surface and in 
the region of the neutron-skin. The di-neutron correlation originates from a coherent 
superposition of quasiparticle orbits with large orbital angular momenta, which are 
embedded in the nucleon continuum. 

Here we illustrate the role of neutron-pair correlations with the example of the 
dipole response in 158 Sn [176J. The HFB calculation of the ground state has been 
performed with the Skyrme effective interaction SLy4, and the density-dependent delta- 
interaction of Eq. (j!5p in the pairing channel. The overall strength parameter of the 
pairing force has been adjusted to reproduce the 1 S scattering length a = —18 fm in free 
space. By varying the value of the parameter 77, which multiplies the density-dependent 
term (see Eq. (TT5|) ). one controls the effective pairing interaction. A smaller value of 
77, i.e. a weaker density dependence, results in stronger pairing in the interior of the 
nucleus. For the particular case of Sn isotopes, 77 = 0.71 has been adjusted so that 
the average pairing gap becomes comparable to the experimental value A = 1.1 — 1.4 
MeV. The choice 77 = corresponds to a density-independent pairing interaction with 
an extremely large and unrealistic pairing gap A « 15 MeV. 

The response function has been calculated with the continuum QRPA (cf. Sec. 12.3ft . 
which employs the quasiparticle Green's function with exact outgoing boundary 
condition for neutrons |176j . For the ph residual QRPA interaction, the Landau-Migdal 
approximation to the Skyrme functional has been used, whereas the pairing correlations 
have been consistently described by the same density-dependent delta- interaction (fl5|) . 
both in the HFB and QRPA. An additional renormalization factor must be introduced 
in the ph channel in order to remove the spurious center-of-mass contributions. 

In Fig. [27] the HFB+QRPA dipole strength function for 158 Sn is shown for different 
choices of the parameter of the pairing interaction: rj = 0,0.3,0.5 and 0.71. For the 
realistic value of rj = 0.71, in addition to the IVGDR in the high-energy region between 
10 and 15 MeV, the strength distribution exhibits a pronounced low-energy structure 
below 5 MeV. In the case of extremely strong pairing (77=0), the transition strength 
basically contains only a single broad resonance centered at 11 MeV. 

The structure of the characteristic low-energy peak at 4.4 MeV, obtained in the 
realistic calculation with 77 = 0.71, is explored in more details in Fig. [281 where the 
corresponding transition densities are shown [176J. A large contribution of the neutron 
particle-pair transition density is found in the region beyond the nuclear surface. It 
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originates from dynamical pairing correlations among neutrons moving in the external 
region, i.e., the QRPA correlation is determined by the pairing interaction. The 
dynamics of this soft mode can be interpreted as the vibration of di-neutrons against the 
core. Because of the influence of neutron pairing correlations, this mode has a dominant 
particle-particle character. The two-quasiparticle configurations, including orbitals in 
the continuum with orbital angular momenta up to I ~ 10, contribute coherently to 
the large particle-pair transition density. Let us also note that in the continuum QRPA 
calculations of Ref. |176j . the quadrupole core vs dineutron mode has also been predicted 
in the low-energy region, but only for a very strong, and therefore unrealistic, pairing 
interaction. 

5. Isoscalar Modes 

5.1. Isoscalar Dipole Compressional and Toroidal Modes 

Experimental excitation energies of compressional (monopole and dipole) vibrational 
modes in atomic nuclei can in principle be used to deduce the value of the nuclear 
matter compression modulus K nm |177j . This quantity is related to the curvature of the 
nuclear matter equation of state at the saturation point, and controls basic properties 
of atomic nuclei, the structure of neutron stars, the dynamics of heavy-ion collisions 
and of supernovae explosions. The nuclear matter compressibility cannot be measured 
directly, but rather deduced from a comparison of experimental excitation energies of 
isoscalar giant monopole resonances (ISGMR), with the corresponding values predicted 
by microscopic nuclear effective interactions characterized by different values of K nm . 
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Figure 28. The transition densities for the dipole peak at E = 4.4 MeV in 
158 Sn, calculated with the value r\ = 0.71 for the parameter of the density-dependent 
pairing interaction [176j . The curves correspond to the neutron particle- hole transition 
density (solid), neutron particle-pair transition density (diamonds), neutron hole-pair 
transition density (crosses), and the proton particle-hole transition density (dashed). 

Inelastic a-scattering experiments have been used in high precision studies of the 
systematics of ISGMR in nuclei with A > 90. Nuclear structure models provide a 
consistent description of the main moments of strength distributions and the mass 
dependence of excitation energies, and thus relate the ISGMR to the nuclear matter 
compression modulus. There is much less experimental information, and only few 
microscopic theoretical analyses of the structure of compressional modes in lighter nuclei 
with A < 90. While in heavy nuclei the shape of the ISGMR strength distribution is 
typically symmetric, for A < 90 the ISGMR display asymmetric shapes with a slower 
slope on the high energy side of the peak, and with a further decrease of the mass number 
the ISGMR strength distributions become strongly fragmented. The determination of 
K nm is based on microscopic calculations of ISGMR excitation energies. Interactions 
that differ in their prediction of the nuclear matter compressibility, but otherwise 
reproduce experimental data on ground-state properties reasonably well, are used to 
calculate ISGMR in the random phase approximation or the time-dependent framework. 
A fully self-consistent calculation of both ground-state properties and ISGMR excitation 
energies restricts the range of possible values for K nm . It has been pointed out, however, 
that, since K nm determines bulk properties of nuclei and, on the other hand, the 
ISGMR excitation energies depend also on the surface compressibility, measurements 
and microscopic calculations of ISGMR in heavy nuclei should, in principle, provide a 
more reliable estimate of the nuclear matter compressibility |177[ 1178] . 

Recent theoretical studies of nuclear compressional modes have employed the 
fluid dynamics approach [179|, the Hartree-Fock plus random phase approximation 
(RPA) pH dHD HH21 [183], the RPA based on separable Hamiltonians [184j, linear 
response within a stochastic one-body transport theory [185J, the relativistic transport 
approach [186j . and the self-consistent relativistic RPA |187l 11881 11891 1190] . Several 
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analyses have emphasized the importance of a fully self-consistent description of ISGMR, 
and confirmed that the low value of K nm = 210 — 220 MeV, previously obtained 
with Skyrme functionals, is an artefact of the inconsistent implementation of effective 
interactions [1831 1191] . The excitation energies of the ISGMR in heavy nuclei are thus 
best described with Skyrme and Gogny effective interactions with K nm w 235 MeV. In 
Ref. [182] it has been shown that it is also possible to construct Skyrme forces that fit 
nuclear ground state properties and reproduce ISGMR energies, but with K nm « 255 
MeV. In Ref. [183] a new set of Skyrme forces was constructed that spans a wider range 
of values of K nm and the symmetry energy at saturation density a^. RPA calculations 
with these forces have shown that the ISGMR data can be reproduced either with 
forces having a softer density-dependent term (the exponent a = 1/6 in Eq. (|14p ) and 
K nm ~ 230 — 240 MeV, or with forces having a stiffer density-dependent term (ct — 1/3) 
and K nm ~ 250 — 260 MeV. Other forces, in particular those characterized by larger 
values of K nm , are associated with unrealistic values of the effective mass, and do not 
reproduce ground-state properties. On the other hand, it appears that in the relativistic 
framework the interval of allowed values for K nm is more restricted. A recent relativistic 
RPA analysis based on modern effective Lagrangians with explicit density dependence 
of the meson-nucleon vertex functions, has shown that only effective interactions with 
K nm = 250 — 270 MeV reproduce the experimental excitation energies of ISGMR in 
medium-heavy and heavy nuclei, and that K nm « 250 MeV represents the lower limit 
for the nuclear matter compression modulus of relativistic mean-field interactions [190J. 

The isoscalar giant dipole resonance (ISGDR) is a second order effect, built on 
3hw, or higher configurations. It corresponds to a compression wave traveling back and 
forth through the nucleus along a definite direction. Recent data on the compressional 
ISGDR in 90 Zr, 116 Sn, 144 Sm, and 208 Pb [T92l [T93l 191 H25] can also be used to 
constrain the range of allowed values of K nm [196J. The problem, however, is that 
the isoscalar El strength distributions display a characteristic bimodal structure with 
two broad components: one in the low-energy region close to the isovector giant dipole 
resonance (IVGDR) (« 2hu), and the other at higher energy close to the electric 
octupole resonance (~ 3/kj). Theoretical analyses have shown that only the high-energy 
component represents compressional vibrations [34 ] 1197] , whereas the broad structure 
in the low-energy region could correspond to vortical nuclear flow associated with the 
toroidal dipole moment [1981 11991 1200] . However, as has also been pointed out in the 
recent study of the interplay between compressional and vortical nuclear currents [199J, 
a strong mixing between compressional and vorticity vibrations in the isoscalar El states 
can be expected up to the highest excitation energies in the region m 3hu. Nevertheless, 
models which use effective interactions with K nm adjusted to ISGMR excitation energies 
in heavy nuclei, also reproduce the overall structure of the high-energy portion of ISGDR 
data [m[mTOl27)Tl[202] . 

Accurate data on compressional modes are becoming available also for lighter nuclei, 
e.g. 56 Fe, 58 Ni, 60 Ni [2031 1204] . Inelastic a-scattering data on the isoscalar monopole 
and dipole strength distributions have been analyzed in the relativistic quasiparticle 
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Figure 29. The RHB+RQRPA isoscalar dipole transition strength in 56 Fe, 58 Ni, and 
60 Ni calculated with DD-ME2 effective interaction. The arrows denote the positions 
of the experimental centroid energies of the low- and high-energy components 203 1 . 



random-phase approximation (RQRPA) with the DD-ME2 effective nuclear interaction 
in the particle-hole channel and the finite-range Gogny force in the particle-particle 
channel |205j . In Fig. [29] we display the strength functions in 56 Fe and 58 > 60 Ni, for the 
isoscalar dipole operator: 

A 

' " 3 - \ {r\r, j V„,(Q,). (79) 



i=l ^ 



where (r 2 ) denotes the ground-state expectation value, and the inclusion of the second 
term in the operator ensures that the strength distribution does not contain spurious 
components that correspond to the center-of-mass motion. In all three nuclei the 
strength is strongly fragmented and distrubuted over a wide range of excitation energy 
between 10 MeV and 40 MeV, in agreement with the experimental results of Ref. [203J. 
In the experiment between 56% and 72% of the isoscalar El strength has been located 
in these nuclei below 40 MeV excitation energy, and some missing strength probably 
lies at higher energies. Similarly to the results obtained for heavier nuclei [34] 11971 1200J. 
the El strength is basically concentrated in two broad structures: one in the region 10 
MeV < E x < 20 MeV, and the high-energy component above 25 MeV and extending 
above 40 MeV excitation energy. Only the high-energy portion of the calculated 
El strength is sensitive to the nuclear matter compression modulus of the effective 
interaction. The thick arrows denote the locations of the experimental centroid energies 
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(mi/m ) in the low- and high-energy regions of the isoscalar El strength in 56 Fe, 58 Ni, 
and 60 Ni [203] . We notice a good qualitative agreement between the calculated and 
experimental centroids in the high-energy region, especially taking into account that the 
isoscalar strength above E x = 40 MeV has not been observed in the experiment. In the 
low-energy region, however, the theoretical centroid energies are systematically below 
the experimental values by «l-4 MeV. This effect is in agreement with previous RRPA 
calculations in heavier nuclei |200] . and supports the picture of pronounced mixing 
between compressional and vorticity vibrations in the intermediate region of excitation 
energies. 

The role of toroidal multipole form factors and moments in the physics of 
electromagnetic and weak interactions has been extensively discussed in Refs. [206J 
and [207] . They appear in multipole expansions for systems containing convection and 
induction currents. In particular, the multipole expansion of a four-current distribution 
gives rise to three families of multipole moments: charge moments, magnetic moments 
and electric transverse moments. The latter are related to the toroidal multipole 
moments and result from the expansion of the transverse electric part of the current. 
The toroidal dipole moment, in particular, describes a system of poloidal currents on 
a torus. Since the charge density is zero for this configuration, and all the turns of 
the torus have magnetic moments lying in the symmetry plane, both the charge and 
magnetic dipole moments of this configuration are equal to zero. The simplest model is 
an ordinary solenoid bent into a torus. 

Vortex waves in nuclei were analyzed in a hydrodynamic model [208] . By relaxing 
the assumption of irrotational motion, in this pioneering study solenoidal toroidal 
vibrations were predicted, which correspond to the toroidal giant dipole resonance at 
excitation energy E x w (50 - 70) /v4 1/3 MeV. It was suggested that the vortex excitation 
modes should appear in electron backscattering. The isoscalar 1~ toroidal dipole states 
were studied in the framework of the time-dependent Hartree-Fock theory by analyzing 
the dynamics of the moments of the Wigner transform of the density matrix [209J, 
and excitations with dipole toroidal structure were also found in semi-classical studies 
based on nuclear fluid dynamics [1981 1199] - The first fully microscopic analysis of 
toroidal dipole resonances (TGDR) was performed in the framework of the relativistic 
RPA [129] [200J. Compressional and toroidal dipole modes were also studied with the 
Quasiparticle Phonon Model, using separable residual interactions with the Nilsson or 
Woods-Saxon mean-field potentials [184] . Continuum RPA calculations with Skyrme 
interactions in Ni isotopes have shown that vortex waves could also occur for excitations 
with higher multi-polarities 2 + , 3~ and 4 + [210J. 

In Fig. [30] we display the RRPA dipole strength distributions in 208 Pb for the 
isoscalar dipole operator (ISGDR) in the upper panel, and for the isoscalar toroidal 
dipole operator (TGDR) 
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Figure 30. Isoscalar dipole (upper panel) and toroidal dipole strength distributions 
(lower panel) for 208 Pb, calculated with the RRPA based on three density-dependent 
interactions with K nm —230, 250, and 270 MeV. The (a, a') experimental data for the 
centroids of the low-energy and high-energy portions of the isoscalar dipole distribution 
are denoted by dark [21 lj and light [194|, 1195] arrows. 

in the lower panel, denotes a vector spherical harmonic, and a are the Dirac 

a-matrices. As in the case of the dipole operator Eq. (179]) . the second term ensures 
that the TGDR strength distributions do not contain spurious center-of-mass motion 
components. To illustrate the correlation between the nuclear matter compressibility 
and the isoscalar dipole response, the strength functions are calculated with three 
different relativistic effective interactions with K nm =230, 250, and 270 MeV, and the 
volume asymmetry at saturation a 4 =32 MeV |190j . The experimental centroids of the 
low-energy and high-energy portions of the dipole strength distributions, extracted from 
small angle a-scattering spectra, are denoted by arrows |195l 1211] . 

The positions of the calculated peaks in the low-energy region (below 20 MeV) 
depend only weakly on the incompressibility, whereas the structures in the high- 
energy region are much more sensitive to the choice of the compression modulus of 
the interaction. Both dipole strength distributions display two broad structures: one at 
low energies between 8 and 20 MeV, and one in the high-energy region 20 — 30 MeV. 
Obviously, a strong coupling between the two isoscalar 1~ modes can be expected. This 
coupling becomes even more evident if one rewrites the expression in square brackets of 
the toroidal operator Eq. ([80]) as |208] 

V x (fx V)(r 3 - - < r 2 > r) Y lfl , (81) 

and compares it with the isoscalar dipole operator of the compression mode (I79"|) . The 
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relative position of the two resonance structures will, therefore, depend on the interaction 
between the toroidal and compression modes. In Table [2] we compare the RQRPA 
centroid energies of the low-lying portion of the response to the toroidal operator, with 
the corresponding experimental values for 116 Sn, 144 Sm, and 208 Pb |1951I211] . We notice 
that the theoretical centroids, calculated with the DD-ME2 effective interaction, are 
systematically located ~ 1 — 2 MeV below the experimental values. The dynamics of 
the solenoidal toroidal vibrations is illustrated in Fig. [3T1 where we plot the velocity fields 
for the three most pronounced peaks of TGDR response function in 208 Pb (calculated 
with DD-ME2). A vector of unit length is assigned to the largest velocity. All the 
other velocity vectors are normalized accordingly. Since the collective flow is axially 
symmetric, we plot the velocity field in cylindrical coordinates. The z-axis corresponds 
to the symmetry axis of a torus. The lowest peak at 7.2 MeV is dominated by 
vortex collective motion. The velocity fields in the {z,r±) plane correspond to poloidal 
currents on a torus with vanishing inner radius. The poloidal currents determine the 
dynamical toroidal moment. The high-energy peak at 26.1 MeV displays the dynamics 
of dipole compression mode. The "squeezing" compression mode is identified by the 
flow lines which concentrate in the two "poles" on the symmetry axis. The velocity 
field corresponds to a density distribution which is being compressed in the upper half- 
plane, and expands in the lower half-plane. The centers of compression and expansion 
are located on the symmetry axis, at approximately half the distance between the 
center and the surface of the nucleus. Finally, the intermediate peak at 10.0 MeV 
displays the coupling between the toroidal and compression dipole modes. A very 
similar behavior of the velocity distributions as function of excitation energy is also 
observed for 116 Sn and 144 Sm. A direct experimental evidence for the TGDR mode 
remains a challenge for future studies. In principle the vortex type of motion could 
be identified in the measurement of transverse electron scattering form factors. An 
exploratory study with the quasiparticle phonon model (QPM) has shown that the 
cross sections in electron back-scattering could differentiate between the toroidal and 
neutron-skin dipole modes [212j . The respective QPM electron scattering form factors 
at 180° are shown in Fig. [32J, for transitions dominated by the toroidal and neutron-skin 
oscillations. Information about the nature of the low-lying dipole excitations could be 
obtained in the range of incident energies 40 — 90 MeV, even though the predicted values 
for the cross sections are low |212j . 





mi /mo (MeV) 


E x (MeV) [195J 


E x (MeV) [211J 


116 Sn 


13.3 


15.6±0.5 


14.38±0.25 


144 Sm 


12.7 


14.2±0.2 


14.00±0.30 


208p b 


11.2 


13.0±0.1 


13.26±0.30 



Table 2. The RQRPA centroid energies of TGDR calculated with the DD-ME2 
effective interaction in the region below 20 MeV excitation energy, compared with 
the corresponding experimental centroid energies from Refs. |195j and [211j . 
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Figure 31. Velocity distributions for the three most pronounced peaks in the TGDR 
response of 208 Pb (calculated with the DD-ME2 effective interaction). The velocity 
fields correspond to the peaks at (a) 7.2 MeV ,(b) 10.0 MeV and (c) 26.1 MeV. 




E q (MeV) 

Figure 32. Electron scattering form factors of dipole transitions at 180° for 
208 Pb, based on calculations with the quasiparticle phonon model |212j . The two 
curves correspond to predominantly toroidal (solid) and neutron-skin (dashed) density 
oscillations, respectively. 



Exotic modes of excitation in atomic nuclei far from stability 



69 





E exp (MeV) 


B(E2) exp (e 2 -fm 4 ) 


M n /M p 
N/Z 


18Q 


2.0 


45±2 


0.9±0.2 


20Q 


1.7 


28±2 


2.2±0.5 


22Q 


3.2 


21±8 


1.4±0.5 


24Q 


> 3.8 







Table 3. Experimental values of excitation energies, electromagnetic transition 
probabilities, and ratios of the transition matrix elements, for the first 2 + states in 
neutron-rich oxygen isotopes. 

5.2. Low-Lying Quadrupole States in Unstable Nuclei 

In most even-even nuclei the first excited state is a J 7r =2 + . For magic nuclei the electric 
transition from the ground state to the first 2 + state reflects directly the structure of 
the shell gap. The systematics of the 2f energies is very important in studies of shell 
evolution in nuclei far from stability. In particular, low-lying quadrupole states are 
closely related to the phenomenon of nuclear superfluidity: the 2f states are built from 
nucleonic configurations located close to the Fermi surface, which is precisely the energy 
region in which the pairing interaction is most effective [H H] . 

In oxygen isotopes the evolution of shell structure can be explored from the stable 
isotope to the neutron drip-line in just a few mass units (cf. also Sec. 14.21) . A number of 
experimental [2IBI EH EISl EUl E3SI and theoretical [291 Effl] f22CT | I22T1 12221 1223] 
studies of neutron-rich oxygen nuclei have been reported. Recent theoretical analyses 
have predicted the appearance of new magic numbers in 22 (iV=14) [223] . and 24 
(N=16). Shell-model calculations [224 J also show a strong gap ~ 4.3 MeV between the 
1g?5/2 and 2su 2 subshells, and thus 22 appears to be a magic nucleus. Accordingly, 
both QRPA [216] and shell-model [2231 1225] calculations predict a decrease of the 
B(E2] 2± — > 0J) from 20 O to 22 O. This has been confirmed by recent experimental 
results [21511218(1219] . Table [3] summarizes the experimental values of excitation energies 
and B(E2) values for the 2+ state in neutron-rich isotopes. 

The energy of the first 2 + state in 22 has been measured at 3199(8) keV [217] . 
compared to 1670 keV in 20 O, and its small B(E2) value of 21(8) e 2 -fm 4 [215] indicates 
a strengthening of the iV=14 shell gap. Even though the 2± state of 24 has not been 
observed directly, it has been shown that its energy must lie above 3.8 MeV, and this 
points to N=1Q as a new shell closure [218] . 28 O, which is a doubly magic nucleus in 
the standard shell model, was found to be neutron unbound [213J. 

Both the Gogny functionals [125] 1226] . and Skyrme interactions with density 
dependent pairing (cf. Eq. (jT5l) ) [29j [301 12161 1227] . have been employed in recent 
QRPA studies of the structure of 2+ states in exotic nuclei. In Fig. [33] we display the 
quadrupole response function in 20 O, calculated with the Skyrme interactions SLy4 [16J, 
SGII [228] and SIII [15]. In addition to the strong 2+ state at w 3 MeV, the pronounced 
structure above 20 MeV corresponds to the isoscalar giant quadrupole resonance (IS 
GQR). Obviously the details of the calculated quadrupole strength function depend on 
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Figure 33. The QRPA isoscalar quadrupole strength function for the 20 O nucleus, 
calculated with the SLy4, SGII and SIII Skyrme forces. 



the choice of the interaction. The choice of the pairing interaction, in particular, plays 
an important role in the calculation of the 2f states. This is illustrated in Fig. [34] 
which compares the isoscalar quadrupole strength functions of 18 O, calculated with two 
different pairing interactions Eq. (ITS"]) . In addition to the density-dependent (surface- 
type pairing) interaction of Eq. (fl5l) . a density- independent interaction (volume- type 
pairing) i.e. r]=0 in Eq. (1151) . has been used in the pairing channel. The strength 
parameter is adjusted to the empirical pairing gap in 18 0. We note that without any 
density dependence in the pairing channel, the calculated excitation energy of the 2+ 
state is in better agreement with the experimental value of ~ 2 MeV (see Table [3]). This 
shows that an analysis of low-lying quadrupole states in neutron-rich nuclei can be used 
to determine the structure and medium dependence of effective pairing interactions [229J. 

The calculated transition densities of the 2f state in neutron-rich oxygen isotopes 
are shown in Fig. [351 While the peak in the proton transition density does not change 
its position with the increase of the number of neutrons, the radial dependence of the 
neutron transition density clearly reflects the formation of the neutron skin, especially 
in 22 and 24 0. The decrease in magnitude of the neutron transition density in 24 can 
be related to the predicted N=16 magic neutron number, which appears because of the 
2si/ 2 subshell closure, i.e. the 2si/ 2 state is more bound in 24 than in 22 0. 

The measured B(E2; 2^ — > 0^) values are nicely reproduced by QRPA calculation, 
except for the problematic 18 [29]. In this nucleus a large discrepancy between the 
empirical and theoretical B(E2) values has been found in several shell-model [225[ 1230] 
and QRPA calculations [2161 1226] . This could be explained by the presence of deformed 
states in the experimental low- lying spectrum of 18 0. It has been suggested that the low- 
lying states in 16 > 17 > 18 contain sizeable admixtures of highly deformed states [231[ [232] . 
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Figure 34. Isoscalar quadrupole strength function of 18 calculated with two types 
of pairing interaction. See text for description. 



By taking into account the mixing between spherical and deformed states, it should 
be possible to simultaneously reproduce the excitation energies and B(E2) values of 
the low-lying states. For heavier oxygen isotopes the deformed states are predicted at 
higher energies, and thus the mixing is weaker. This explains why the spherical QRPA 
results for the B(E2) values in 20 ' 22 O are in much better agreement with experiment. 
The calculation predicts a decrease in the B(E2) value in 24 0, and this is due to the 
effect of the 2si/ 2 subshell closure. 

In nuclei at the neutron drip-line the Fermi level of neutron states is found close 
to the continuum and, therefore, one expects that continuum effects play an important 
role also in the calculation of low-lying states. This is illustrated in Fig. [36] where we 
plot the low-energy portion of the quadrupole response function of 22 O. The two curves 
correspond to calculations with box boundary conditions, and with the exact treatment 
of continuum states (cf. Sec. 2). The latter predicts the 2 + state at a slightly higher 
energy, and with a weaker transition strength. In Ref. [30] a continuum QRPA based 
on the HFB framework in coordinate space has been formulated, and the quadrupole 
response of the drip-line nucleus 24 has been described using density-dependent zero- 
range forces in the particle-hole and particle-particle channels. It has been shown that 
the low-lying isoscalar quadrupole state is embedded in the neutron continuum, and its 
excitation energy and strength are very sensitive to the density dependence of pairing 
correlations. 

In order to disentangle the proton and neutron contributions to the 2 + excitations, 
the following reduced matrix elements will be useful: 

z 

M p = (2+11^^(^)110), 
i=i 
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Figure 35. Neutron and proton transition densities for the first 2 + state in 18 < 20 < 22 ,24q^ 
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Figure 36. Isoscalar quadrupolc strength function of 22 0, calculated with the 
continuum QRPA (solid), and with a box discretization procedure (dashed). 



N 

M n = (2 + \\Y,r*Y 2 (r t )\\0), (82) 
i=i 

so that for an electromagnetic probe the B(E2) is ~ M 2 , whereas both M n and M v 
contribute in the case of hadron scattering. For instance, experimental evidence for the 
magicity of the N=14 neutron number in 22 O cannot be conclusive without separating 
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the proton and neutron contributions to the 2f state. In Ref. [215] the B(E2) value 
for 2\ was determined by inelastic scattering of 22 O from 197 Au at an energy of 50 
MeV/nucleon. The extracted value, however, depends on theoretical predictions for the 
ratio between M n and M p , because both Coulomb and nuclear interactions contribute to 
the reaction. M n and M p can be separated by means of two experiments which employ 
different probes. An electromagnetic probe is used to measure the B(E2) value directly, 
whereas the second measurement is usually a (p,p') scattering experiment at around 50 
MeV/nucleon. This combination allows to determine both M n and M p , and therefore 
to probe more directly possible shell closures in exotic nuclei. Angular distributions 
for elastic and inelastic proton scattering to the 2\ state of 22 have been measured 
using a secondary radioactive beam |219j . Proton and neutron contributions have been 
disentangled by a comparison of the (p,p') results with a heavy ion scattering experiment 
dominated by electromagnetic excitation, and evidence for a strong N=14 shell closure 
has been found. 

In order to compare the QRPA predictions with proton scattering data, microscopic 
optical potentials can be generated from the HFB and QRPA densities using two 
different methods: the folding model [233J, or the optical model potential (OMP) 
parameterization using the Jeukenne, Lejeune and Mahaux (JLM) interaction [234J. 
The folding model analysis uses the CDM3Y6 interaction folded with the HFB densities 
to generate the isoscalar and isovector parts of the OMP. The spin-orbit potential and 
the transition potentials are determined by folding the QRPA transition densities with 
the nucleon-nucleon interaction. The imaginary part of the OMP is generated with 
the Koning and Delaroche [235j phenomenological parameterization. Cross sections are 
calculated using DWBA with the ECIS97 [236] code. 

The elastic angular distribution is nicely reproduced, even at large angles [21911227] . 
Since the B(E2) can be described by the proton transition density, the neutron transition 
density is renormalized to reproduce the data. This procedure assumes that the QRPA 
provides a reliable description of the shape of the transition density for collective states, 
and provides an empirical value of the M n /M p ratio for the 2f state, deduced from the 
combination of the electromagnetic and the (p,p') measurements |216j . The resulting 
M n /M p values are included in Table We note that the value of M n /M p divided 
by N/Z, is considerably smaller for the first 2 + state of 22 0, compared to 20 O. In 
22 the contributions of protons and neutrons are comparable, because the measured 
(M n /M p ) ratio is close to N/Z. This is different from 20 O, where the much higher 
value of M n /M p divided by N/Z shows that neutrons predominantly contribute to the 
quadrupole excitation. Combined with the high energy of the 2\ state in 22 O, these 
results confirm the N=14 shell closure in neutron- rich oxygen nuclei. The dependence 
on the potential used to describe the (p,p') angular distributions can be checked by using 
the complex optical and transition potentials obtained by inserting the calculated ground 
state and transition densities into the JLM density-dependent optical potential [234J. 
Renormalizing the neutron transition density to reproduce the inelastic data, the same 
value of the ratio M n /M p is obtained as with the folding potential. In this way two 
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reliable optical potentials are used to test the HFB+QRPA matter and transition 
densities. 

A similar study of the N=16 sub-shell closure in 24 could be performed in a 
24 0(p,p') experiment, but this will have to wait for the next generation of radioactive 
beam facilities. The generality of the method used to microscopically determine the ratio 
M n /M p from two complementary sets of data, will allow to extend our understanding of 
neutron-shell closure to regions of heavier nuclei far from stability. An ingenious proton- 
scattering setup with a liquid hydrogen target and 7-detectors opens the possibility to 
determine M n /M p in nuclei very far from stability, e.g. in 28 Ne [237J. 
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Figure 37. The QRPA excitation energies (left panels) and transition strength 
(right panels) for the 2\ and states in the sequence of even-even Sn isotopes 
(squares) J238], compared with data (circles). 

The Skyrme-HFB plus QRPA framework has been applied to the study of low-lying 
quadrupole states in N = 20 isotones, including the neutron-rich nuclei 32 Mg and 30 Ne 
[127] . The calculation reproduces both the excitation energies of the first 2 + states and 
the B(E2; Oj" — ► 2+) values, and it has been concluded that pairing effects account 
largely for the anomalous B(E2) value and the very low excitation energy in 32 Mg. 

QRPA calculations for the low-lying 2 + and 3~ states have also been performed 
for heavier neutron-rich nuclei. For example, the QRPA with Skyrme forces in the ph 
channel, and a density-dependent pairing force in the pp channels, has been employed 
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in the study of the evolution of the 2f and 3]~ states along the Sn isotopic chain |238j . 
The results are shown in Fig. [371 i n comparison with experimental data [2 39] , including 
recent Coulomb excitation measurements of the unstable 126 > 128 > 130 Sn isotopes [240] . We 
notice that the calculation nicely reproduces the excitation energies of the 2f state, 
but overestimates the experimental energies of the 3~[ state by more than 50%. On 
the other hand, the measured B(E2) values in open shell nuclei are systematically 
underestimated by large factors with the QRPA model, whereas this discrepancy is less 
pronounced for the B(E3) values. The reason for such a large difference between the 
experimental and calculated B{E2) values, is that the QRPA does not take into account 
the anharmonicities of the low-lying quadrupole vibrations, i.e. the 2f is described as 
a pure 1-phonon state, whereas in this mass region the 2f states contain pronounced 
admixtures of 2-phonon states which cannot be described in a simple QRPA framework. 

In the recent Skyrme-HFB plus QRPA analysis of the J n = + , 1~, and 2 + 
multipoles for the even Ca, Ni and Sn isotopes from the proton to the neutron drip-lines 
[120J, it has been shown that the strength functions in the 2 + channels are qualitatively 
different from those in the lower-multipole channels. The low-energy strength grows with 
the neutron number, but unlike in the + channel, both neutrons and protons contribute 
to the transitions to these states, even near the drip line, because the isoscalar (IS) peaks 
are much larger than the isovector (IV) peaks. An interesting result is also that at the 
neutron drip-line the IS and IV strength functions in Ni and Sn have distinct low-energy 
peaks, whereas in Ca the peaks coincided. In all nuclei near the neutron drip-line, the 
states in the low-energy peaks are mostly above the neutron-emission threshold, yet the 
neutron tails cut off at much smaller radii than do those in the + and 1~ channels. 
Also the transition densities to the 2 + states are different: they have no real nodes and a 
proton component that is of the same order as the neutron component. The strong low- 
energy states near the neutron drip-line have transition densities that resemble those 
of surface vibrations are often quite collective. The detailed analysis of Ref. [120J has 
emphasized the complicated relationship among collectivity, strength, and transition 
density in neutron-rich nuclei. 

5. 3. Giant Quadrupole Resonance and Higher Multipoles 

The giant quadrupole resonance (GQR) corresponds to a highly collective oscillation 
of the neutron and proton density distributions between prolate and oblate ellipsoidal 
shapes. In the isoscalar mode the proton and neutron densities oscillate in phase, with 
an empirical excitation energy: Egqr ~ 64 MeV-A -1 / 3 and, in heavy nuclei, this mode 
typically exhausts 50 to 100 % of the energy weighted sum rule (EWSR) [Bj. The 
isovector GQR is found at much higher excitation energies and, even in heavy nuclei, it 
is much more fragmented. 

In neutron-rich nuclei far from stability one expects that the neutron skin has a 
pronounced effect on the high-energy quadrupole vibrational mode, however, no data 
on the GQR in unstable nuclei are available at present. There have also been only few 
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E GQR (MeV) 


Tfwhm (MeV) 


34 Ca 


17.5 


1.6 


40 Ca 


16.2 


0.6 


48 Ca 


16.7 


0.3 


60 Ca 


14.9 


1.3 



Table 4. Continuum-RPA results for the GQR centroids and widths of Ca 
isotopes |118j 

theoretical studies of giant quadrupole resonances in exotic nuclei. In Refs. [29] EH EZ] 
the QRPA and RQRPA calculations have been performed for the isoscalar GQR in 
oxygen isotopes. Since in light systems this mode can be highly fragmented, for a study 
of the effect of neutron excess on the GQR one should consider heavier nuclei, in which 
the GQR displays a single peak at energies between 10 and 15 MeV. Such studies have 
been carried out with the continuum- RPA using Skyrme forces [ 118 ], and recently RPA 
calculations with the Gogny interaction have also been performed [125j . 

The isospin dependence of the excitation energy and width of the GQR in exotic 
nuclei is illustrated in Table H] for the example of 34 > 40 > 48 > 60 Ca [118j . In addition to the 
systematic lowering of the GQR excitation energies with the increase of the number 
of neutrons ( 48 Ca is an exception, because of the neutron shell closure at N = 28), 
the calculation predicts an enhancement of the low-energy quadrupole strength both 
in the proton-rich 34 Ca, and neutron-rich 60 Ca. Continuum-RPA calculations predict 
the escape width of a resonant state, and in Table H] we notice a pronounced increase 
of the escape width in the weakly-bound nuclei 34 Ca and 60 Ca. The calculated widths, 
however, do not contain the spreading contribution, i.e. the width that results from the 
coupling of the GQR to more complex states like, for instance, two-phonon admixtures. 
In this sense the widths in Table H] are not realistic, and only illustrate the effect of the 
coupling to continuum states. No systematic calculation of the GQR spreading width 
in exotic nuclei has been reported so far, the only exception is the study of the GQR in 
28 in Ref. [41J, which has shown that the spreading width is enhanced with respect to 
the isotopes close to the stability line. 

Giant quadrupole resonances in exotic nuclei have also been calculated with the 
RPA based on the Gogny interaction [125J. For the doubly-magic 78 Ni, 132 Sn and 100 Sn 
the Gogny- RPA results predict GQR excitation energies that are 1 — 1.5 MeV above the 
empirical relation Ecqr ~ 64 MeV-A -1 / 3 . In Ref. [125] RPA-Gogny calculations have 
also been performed for the isoscalar and isovector octupole strength distributions, and 
in Fig. [38] we display the isoscalar octupole states in 78 Ni, 132 Sn, 100 Sn and 208 Pb. The 
strength is clearly separated into two regions: the 3^~ state dominates the low-energy 
strength, whereas the strong peaks above 20 MeV excitation energy correspond to the 
High Energy Octupole Resonance (HEOR). The properties of the HEOR in nuclei far 
from stability ( 78 Ni, 132 Sn, 100 Sn) are not significantly different from those in stable 
nuclei, e.g. 208 Pb, probably because the HEOR corresponds to 3ftw excitations and 
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these high-energy configurations may not be very sensitive to changes in the number of 
neutrons/protons with respect to stable nuclei. 
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Figure 38. Fraction of the EWSR exhausted by the isoscalar octupole states in 78 Ni, 
132 Sn, 100 Sn and 208 Pb (I55I . 



5.4- Pairing Vibrations in Drip-Line Nuclei 

Two-neutron transfer reactions such as (t, p) or (p, t) have been used for many years in 
studies of nuclear pairing correlations |241] . The corresponding pair-transfer modes are 
usually described in terms of pairing vibrations or pairing rotations |242[ I243] . Pairing 
vibrations are L = modes induced by the addition or removal of a pair of neutrons, 
and can be associated with the fluctuation of the pairing field A (see Sec. I2.ll) around its 
equilibrium value. Around magic nuclei, such as 208 Pb, the J 7r =0 + spectrum generated 
by the pair- addition (cjct) and pair-subtraction (cjCj) operators in Eq. (I37I) . is harmonic 
and corresponds to the so-called pairing vibration mode, which can be viewed as a 
vibration in an abstract "gauge" space, instead of the ordinary three-dimensional space 
in which shape vibrations take place. 

In a microscopic approach the collective two neutrons transition can be described 
in the QRPA framework: the excitation operator of Eq. (1371) includes particle-hole, 
particle-particle and hole-hole excitations. In this case the non-conservation of the 
particle number, which is implicit in the quasiparticle formalism, can be used as a tool 
to study particle-violating transitions. The pairing vibrational state reads 

\A + 2,n)=(j2 4 B) <M - E \A > ■ (83) 

\k F <k<l k<l<k F J 

High-energy collective pairing modes - Giant Pairing Vibrations (GPV) have been 
predicted and studied theoretically |241j 12431 I244[ I245j I246j . but have never been 
observed despite a number of (p, t) experiments performed in the 70's and 80's |241U243j . 



Exotic modes of excitation in atomic nuclei far from stability 



78 



The experimental setup for such studies must achieve a proper balance between the low 
energy of the incident proton beam (below 50 MeV) necessary for the excitation of the 
L = mode, and the energy which is required to populate the energy region of the GPV. 
These experiments are therefore rare [247] , and the discovery of deep hole-states of non- 
collective character made the detection of the GPV even harder [248L [249J. Beams of 
exotic nuclei could provide the solution. Incident beams of a few MeV/nucleon of weakly- 
bound projectiles (providing high Q-values) could be used to populate the GPV [241] 
using reactions such as ( 6 He,a). However, the intensity of radioactive beams is typically 
several orders of magnitude lower than that of stable beams, and the background 
resulting from break-up reactions could be large with weakly-bound projectiles. There 
is a renewed interest in GPV and improved experimental investigations are currently 
being planned, both with stable and exotic beams. 

Pairing vibrations in exotic nuclei could provide valuable structure information [241] 
124411250] , in particular on pairing correlations in systems far from stability and the effects 
of the coupling to the continuum. Pairing vibrations generally depend on the strength of 
the pairing interaction between the two transferred neutrons: transfer cross sections are 
enhanced when the two neutrons form a strongly bound pair. The theoretical analysis 
of two-nucleon transfer modes in nuclei far from stability is complicated by the effect 
of continuum coupling. The right tool to study pairing vibrations in exotic nuclei is 
the continuum-QRPA, because it provides a consistent microscopic treatment of both 
pairing and continuum effects. The strength function which describes the two-particle 
transfer from the ground state of a nucleus with A nucleons to the excited states of a 
nucleus with A+2 nucleons reads: 



where G 22 denotes the {pp,pp) component of the Green's function, and F 12 is the 
perturbing external field (see Sec. 12.31) . 

The strength function for the neutron-pair transfer to 22 O is shown in Fig. [39] [250J. 
The interaction in the ph channel is Skyrme SLy4, and a zero-range density-dependent 
pairing interaction is used in the pp channel. The solid curve denotes the unperturbed 
strength, which displays peaks characteristic for the addition of two neutrons in specific 
configurations. The first peak results from the filling of the 2sy2 orbital, and the 
peak at 10.8 MeV corresponds to the two-neutron (lc^) 2 quasiparticle configuration. 
The effect of the residual interaction on the pair transfer mode is seen in the QRPA 
strength function (dashed curve), and it demonstrates the collective nature of pairing 
vibrations [243J. Namely, the residual interaction shifts the position of the two- 
quasiparticle resonant state, located at 10.8 MeV in the unperturbed response, to lower 
energy and increases its strength. The strong peak at zero-energy corresponds to the pair 
transfer to the spurious Goldstone mode associated with particle-number fluctuations. 
The peak at ~ 8.6 MeV represents a neutron-pair transfer predominantly to the 
states, whereas the broad resonant structure around 16 MeV is built mainly on the 
neutron resonant state 1/7/2- This broad two-quasiparticle resonance is characteristic 




(84) 
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for a giant pairing vibration [244[ I245[ 246J . 
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Figure 39. The strength function for the two-neutron transfer on 22 0. The solid 
(dashed) curve corresponds to the unperturbed (QRPA) response. 



To compare the results with the measured cross section, one must calculate the 
nuclear structure form factor, which describes the wave function of the two transferred 
neutrons, as well as the reaction part which includes the optical potential. For closed- 
shell nuclei the two-particle transfer modes are described by the particle-particle (pp) 
RPA [25H [252], whereas the QRPA [IM3l [211] is used for open-shell nuclei. Most of the 
cross section calculations employ the Distorted Wave Born Approximation (DWBA). 

In nuclear response theory the transition from the ground state to the excited state 
\u) in the same nucleus is described by the transition density 

5p" (rcr) = <0|c f (rex) c (rex) \v) (85) 

where c 1 (rcr) is the particle creation operator in coordinate space. The corresponding 
quantity used in the description of pair transfer processes is the pair transition density 
defined by 

5k u (rcr) = (0|c (rcr) c (rcr) \v) (86) 

where the operator c' (rcr) = — 2erc' (r — cr) creates a particle in the time-reversed state. 
The pair transition density determines the transition from the ground state of a nucleus 
with A nucleons to the state \v) in the nucleus with A+2 nucleons. This quantity is 
calculated in the QRPA. 

The DWBA calculation of the cross section for the two-neutron transfer requires 
the form factor which describes the correlation between the two neutrons and the initial 
nucleus |253j . Calculations based on the continuum- QRPA include effects of both 
pairing correlations and continuum coupling. The calculation of the form factor and 
cross section will be illustrated for the transfer reaction 22 0(t,p). The DWBA cross 
section is calculated for the 22 0+t Becchetti and Greenlees optical potential [254] in the 
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entrance channel, and the 22 0+p Becchetti and Greenlees |255] in the exit channel. The 
calculation is performed in the zero-range approximation, in which the two-neutrons and 
the residual fragment are located at the same point in space, and therefore the range 
function is determined by a simple constant D |253j . The zero-range approximation 
provides a satisfactory description of the shape of the angular distribution [2431 12531 1256J. 
but its magnitude is generally underestimated. The form factor for the pair transfer 
is obtained by folding the pair transition density 5k v (Eq. (186]) ) with the interaction 
acting between the transferred pair and the residual fragment [241 J. In the zero-range 
approximation one uses the 5-force for this interaction, and therefore the form factor 
coincides with the pair transition density ( 186]) [253J. In order to illustrate the effect 
of the continuum on the form factor, the (t,p) angular distribution for the mode at 
Rj 8.6 MeV (see Fig. [39]) is calculated using both box boundary conditions and with the 
exact treatment of the continuum. The resulting two-neutron transfer cross sections are 
shown in Fig. HQ] and we notice a pronounced effect of the continuum in the diffraction 
minima. 

The continuum-QRPA could be used in the analysis of the forthcoming data on 
the GPV, obtained with beams of exotic nuclei at around 5 MeV/nucleon. Preliminary 
calculations [244j have shown that ( 6 He,a) reactions can excite the GPV with a cross 
section of the order of few millibarns, whereas this mode is not excited in ( 14 C, 12 C) 
reactions. 
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Figure 40. Two-neutron transfer cross section on 22 0, for the state at w 8.6 MeV. The 
solid (dashed) curve corresponds to the QRPA result obtained with the box (exact) 
boundary conditions. 
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6. Charge-exchange Resonances 

6.1. Proton- Neutron QRPA 

Starting from the ground-state of the (N,Z) nuclear system, charge-exchange excitations 
induced by isospin lowering r_ and raising r + operators correspond to transitions to final 
states in the neighboring (N=pl,Z±l) nuclei, respectively. These transitions can occur 
either spontaneously in /3-decay, or they can be induced in charge-exchange reactions, 
e.g. (p, n) or ( 3 He,t) in the r_ case, when the final state lies outside the /9-decay energy 
window. Charge-exchange resonances in stable nuclei have been the subject of numerous 
experimental and theoretical studies (see, for instance, Ref. [6]). However, a systematic 
experimental information on these resonances is still missing, mainly because of the 
absence of selective probes. For instance, there is very little data on the charge-exchange 
dipole resonance and, despite considerable experimental effort, the charge-exchange 
monopole mode has not yet been observed. This resonance is particularly interesting, 
because its excitation energies could provide useful information on the isospin mixing 
in nuclear ground-states, the symmetry term of the nuclear equation of state, and the 
isovector terms in effective nucleon-nucleon interactions. Charge-exchange transitions 
are important for weak-interaction processes in nuclei, in particular charged current 
processes in nuclear astrophysics and neutrino physics. These include the /3-decay of 
nuclei that lie on the r-process path of stellar nucleosynthesis, and neutrino-nucleus 
scattering. 

The only exception are the Isobaric Analog Resonance (IAR) and the Gamow- 
Teller Resonance (GTR), for which rather detailed experimental information has become 
available in several mass-regions. They are induced by the following transition operators 

A 

F IAR = J2r-(i), (87) 
i=i 

and 

A 

F GTR =J2 a ( i> > T -W ' ( 88 ) 

i=l 

respectively. Since these are L = modes, they can be selectively excited by using zero- 
degree charge-exchange reactions, like (p,n) or ( 3 He,t). When the incident projectile 
energy is increased, the excitation of the GTR is favored over the IAR. A comprehensive 
review of the properties of the IAR and GTR in stable nuclei can be found in Ref. [257] . 

Charge-exchange resonances in exotic nuclei are virtually unexplored. The 
experimental investigation of these modes in nuclei far from stability will become 
possible with the new generation of radioactive-beam facilities. So far, experimental 
studies have been restricted to the states which are accessible in /3-decay. In neutron-rich 
nuclei, for example, the neutron Fermi level is located at much higher energy than the 
proton valence states, and the /3~-decay becomes progressively much faster. In proton- 
rich nuclei the same happens for the /3 + -decay. Closer to the nucleon drip-lines more and 
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more single-nucleon states enter the /3-decay energy window, and the experiment probes 
most of the charge-exchange strength. The importance of this phenomenon for the 
understanding of the isospin properties of nuclei has been first pointed out in Ref. [258J, 
where the HF-RPA has been employed in the study of charge-exchange transitions in 
weakly-bound systems. 

The extension of the Skyrme-RPA model to the charge-exchange channel [259J, 
has been used in calculations of the response to different operators [260] . The discrete 
RPA has also been employed in studies of charge-exchange modes. In addition, the 
formalism of the coupling of p — h configurations to more complex states of 2p — 2h 
character, has been developed in Refs. [5H, 12611 [262J. In the case of open-shell nuclei, 
most of the charge-exchange QRPA calculations have used simple separable forces as 
the ph and pp residual interactions. The same functional form has been used in both 
channels, with two different strength parameters g p h and g pp , as in the pioneering work 
of Ref. [263J where the formalism has been developed for the first time. More recently, 
the PN-QRPA based on Skyrme forces in the ph channel, with a simplified residual 
interaction of a separable form, and with the BCS treatment of the pp channel, has also 
been extended to the description of deformed nuclei. /3-decay rates in different isotopic 
chains have been studied |264l 12651 1266] 1267] . in order to determine to what extent 
the intrinsic deformation influences /5-decay properties, and whether the decay spectra 
differentiate between spherical, prolate and oblate shapes. Medium-heavy nuclei like 
Kr, Sr have been studied, as well as heavier systems e.g. Hg, Pb and Po [268J. It has 
been shown that prolate and oblate deformations lead to significantly different Gamow- 
Teller spectra, but in many cases these signatures of the intrinsic deformation are not 
too sensitive to the choice of the Skyrme force used in the calculation. 

The path of the r-process nucleosynthesis is governed by delicate balance between 
neutron-capture reactions and /3-decay rates. This has motivated the calculation of the 
GT /3-decay of the so-called "waiting-point" nuclei in Ref. [31], which was also the first 
attempt to perform a self-consistent PN-QRPA analysis based on the HFB treatment 
of the nuclear ground state. A Skyrme interaction and a zero-range pairing interaction 
were used in the mean-field and pairing channels, respectively. The associated QRPA 
equations were solved in the canonical basis, in which they read as in Eq. ( l55l) . with the 
matrices A and B given by: 



Here p, p f , and n, n' denote proton and neutron canonical states, respectively, V ph is 
the proton-neutron ph residual interaction, and V pp is the corresponding pp interaction. 



At 




bL 



(89) 



Since the canonical basis does not diagonalize the mean-field Hamiltonian h, nor the 
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Figure 41. IAR excitation energy in the sequence of Sn isotopes. The results from 
the fully self-consistent Skyrme QRPA calculation are compared with the experimental 
excitation energies [271] , 



pairing field A, the non-diagonal matrix elements H^ n , and H^, appear in the matrix 
A: 

Hll, = (u K u K , - v K v K i)h KK , - (u K v K > + v K u K >)A KK , . (90) 

The PN-QRPA calculation of Ref. [31] has shown that the excited 1 + states and /3- 
decay half-lives are very sensitive to the T = component of the residual pp interaction. 
Since the T = pairing is not manifest in the ground states of nuclei with N different 
from Z, it has to be introduced in the PN-QRPA independently from the T = 1 
channel, and this of course breaks the self-consistency of the calculation. In particular, a 
finite-range interaction with adjustable parameters has been employed, and the overall 
strength was tuned to reproduce some selected experimental /3-decay half-lives. The 
same approach has been used in Ref. p69j to analyze properties of GT resonances 
predicted by several Skyrme parameterizations, in correlation with the corresponding 
values of the Landau parameters in infinite nuclear matter. 

Very recently a fully self-consistent non-relativistic charge-exchange QRPA has 
been developed and tested in the description of IAR |270j . In this case only the 
T = 1 component of the residual pairing interaction contributes and, because of isospin 
invariance, its strength must be identical to that used for the calculation of HF-BCS 
ground-state. In Fig.HHwe show the excitation energies of IAR of Sn isotopes, calculated 
with the PN-QRPA based on the Skyrme interaction SLy4 and compared with the 
experimental values. We notice a very good agreement with data, and this shows that 
the model can be systematically employed in calculation of charge-exchange processes 
in stable and exotic open-shell nuclei. 



Exotic modes of excitation in atomic nuclei far from stability 



84 



6.2. Proton- Neutron Relativistic QRPA 

In the relativistic mean-field framework the spin-isospin dependent interaction terms 
are generated by the p- and 7r-meson exchange. Because of parity conservation, the one- 
pion direct contribution vanishes in the mean-field calculation of a nuclear ground state. 
Its inclusion is important, however, in calculations of excitations that involve spin and 
isospin degrees of freedom. The particle-hole residual interaction in the PN-RQRPA is 
derived from the Lagrangian density 

U - 

£tt+ p = -gp^YPy.^ -ipjsYd^Til) . (91) 

In Ref. |272j the proton-neutron relativistic quasiparticle RPA (PN-RQRPA) has 
been formulated in the canonical single-nucleon basis of the time-dependent RHB model, 
and applied in studies of charge-exchange excitations in open shell nuclei. The spin- 
isospin dependent residual two-body interaction in the particle-hole channel reads 

nf 2 (SiVi)(S 2 V2)D w (n,r 2 ) , (92) 



m t 

where D p ^) denotes the meson propagator 

1 e - m p(7r)|ri-r 2 | 



and 



iV)(ri,r 2 ) = ^ [ri _ r2 | > (93) 



(94) 




Transitions between the + ground state of a spherical even-even parent nucleus, and 
the state with angular momentum and parity J n of the odd-odd daughter nucleus are 
considered. With respect to the RHB calculation of the ground state of the even- 
even nucleus, the charge-exchange channel includes the additional one-pion exchange 
contribution. For the pseudovector pion-nucleon coupling the standard values: m w = 
138.0 MeV, and f% NN /4:iT = 0.08, are used. The derivative type of the pion- 
nucleon coupling necessitates the inclusion of the zero-range Landau-Migdal term, which 
accounts for the contact part of the nucleon-nucleon interaction 

VW(ri,r 2 ) = g ( — J fif 2 Ei • E 2 $(n - r 2 ) , (95) 

with the parameter g' 0.6 usually adjusted to reproduce data on excitation energies of 
Gamow- Teller resonances |272[ I273I I274] . In the non-relativistic limit the corresponding 
two-body interaction reduces to the familiar form G' a\Ti ■ o^^- We note that a ph 
residual interaction based on p- and 7r-meson exchange has also been used in a number 
of non-relativistic RPA studies of charge-exchange excitations [2751 1276l I277j . 

In the particle-particle channel of the PN-RQRPA equations both the T = 1 and 
T = pairing interactions contribute. In Ref. |272j the finite-range Gogny force with 
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the parameter set D1S [T3] has been used in the T = 1 channel, both in the RHB 
calculation of the ground state of the even-even system and as the PN-RQRPA pp 
residual interaction. For the T = proton-neutron pairing force a similar interaction 
has been employed: a short-range repulsive Gaussian combined with a weaker long-range 
attractive Gaussian 



where Il5 =1)T= o projects onto states with S = 1 and T = 0. This interaction was also 
used in the non-relativistic QRPA calculation of /3-decay rates for spherical neutron- 
rich r-process waiting-point nuclei [31] . The ranges of the two Gaussians /ii=1.2fm and 
yU2=0.7fm are taken from the Gogny interaction Eq. ffl3l) . and the choice of the relative 
strengths: g% = 1 and #2 = — 2 makes the force repulsive at small distances. The 
overall strength parameter Vq can be adjusted, for instance, to experimental /3-decay 



The two-quasiparticle PN-RQRPA configuration space includes states with both 
nucleons in the discrete bound levels, states with one nucleon in the bound levels and 
one nucleon in the continuum, and also states with both nucleons in the continuum. 
In addition to configurations built from two-quasiparticle states of positive energy, 
the RQRPA configuration space contains pair-configurations formed from the fully 
or partially occupied states of positive energy and the empty negative-energy states 
from the Dirac sea. The inclusion of configurations built from occupied positive-energy 
states and empty negative-energy states is essential for the consistency of the relativistic 
(proton-neutron) QRPA (current conservation, decoupling of spurious states, sum rules). 
The PN-RQRPA model is fully consistent: the same interactions, both in the ph and 
pp channels, are used in the RHB equation that determines the canonical quasiparticle 
basis, and in the PN-RQRPA equation. In both channels the same strength parameters 
of the interactions are used in the RHB and RQRPA calculations. 

6.3. Gamow- Teller Resonances 

Collective spin and isospin excitations in atomic nuclei have been the subject of many 
experimental and theoretical studies (for an extensive review see Ref. [257J). Nucleons 
with spin up and spin down can oscillate either in phase (spin scalar S — mode) or 
out of phase (spin vector 5=1 mode). The spin vector, or spin-flip excitations can 
be of isoscalar (S = 1, T = 0) or isovector (S = 1, T — 1) nature. These collective 
modes provide direct information on the spin and spin-isospin dependence of the effective 
nuclear interaction. 

The Gamow- Teller resonance represents a fundamental charge-exchange mode and 
corresponds to a collective spin-isospin state J n = 1 + formed when the excess neutrons 
coherently change the direction of their spin and isospin without changing their orbital 
motion. This collective mode was predicted already in 1963 [279] . but it was only 
in 1975 that the first experimental indications of the GT resonance were observed in 



2 




(96) 



half-lives [311 1278]. 
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(p, n) charge-exchange reactions at intermediate energies |28Uj . More recently, advanced 
experiments with ( 3 He,t) reactions have achieved very good energy (AE « 50 keV) 
and angular (AG ~ 0.3°) resolution [281J, allowing not only the detection of the fine 
structure of GT transitions [282] . but also studies of GT strength in exotic proton-rich 
nuclei [283J. The detailed knowledge of GT strength distributions in regions away from 
the valley of /3-stability is essential for an understanding of nuclear processes relevant 
for nucleosynthesis. In particular, the low-lying GT strength is directly related to (3- 
decay rates, as well as to the electron-capture process leading to the stellar collapse and 
supernovae explosion. At present, however, very little is known about charge-exchange 
excitations in exotic nuclei, and the modeling of nuclear weak-interaction processes in 
stars must rely on theoretical predictions of the GT strength distributions. 

Recent theoretical studies of Gamow- Teller excitations and /3-decay rates have been 
based on: (i) the shell-model approach, (ii) the non-relativistic proton-neutron (Q)RPA, 
and (iii) the relativistic proton-neutron (Q)RPA. Data on charge-exchange excitations 
in light and medium-mass nuclei are very successfully reproduced by large-scale shell- 
model calculations [284] . However, as the number of valence nucleons increases, the 
dimension of shell-model configuration space becomes far too large for any practical 
calculation. Present shell-model calculations are thus restricted up to the region of 
p/-shell nuclei with A = 45 — 65 [285] . Both the GT~ and GT + response, as well as 
/3-decay rates in medium-mass nuclei have been successfully described by shell-model 
Monte Carlo calculations (SMMC) [2861 1287] . Studies of charge-exchange excitations 
have only recently been reported in the relativistic mean-field plus RPA framework. 
These include the relativistic RPA analysis of isobaric analog states and Gamow- Teller 
resonances (GTR) in the doubly closed-shell nuclei 48 Ca, 90 Zr, and 208 Pb [2731 1288] . 
performed in a restricted configuration space that did not include negative-energy states 
from the Dirac sea. Although these configurations do not affect the excitation energies 
of the charge-exchange modes, they have a pronounced effect on the Gamow- Teller sum 
rule [2721 1289] . GT resonances in doubly closed shell nuclei have also been studied with 
the relativistic RPA in the response function formalism [274J. 

The PN-QRPA can be employed in calculations of charge-exchange excitations in 
mass regions that are presently beyond the reach of the most advanced shell-model codes, 
both in stable and exotic nuclei [290] . Most QRPA studies of the GTR have been based 
on Skyrme effective interactions with BCS-type pairing [51 [2581 EE [29H [2921 [293] . 
It has been shown that the choice of the spin-isospin terms of the Skyrme energy 
functional affects the calculated strength distribution and excitation energy of the 
GT resonance, i.e., the properties of the GTR are not entirely determined by the 
Landau-Migdal residual interaction [269] . In addition, the inclusion of particle-particle 
correlations in the QRPA residual interaction is important in calculations of the GT 
transition strength [291 ESH1 [296], /?-decay rates [HH [2971 EHE] and double /5-decay 
amplitudes [299, 300J. The inclusion of proton-neutron pairing changes significantly the 
rates of the neutrinoless double /3-decay, allowing for larger values of the expectation 
value of light neutrino masses [301] . Both the two- neutrino and neutrinoless double 
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/3-decay matrix elements are suppressed by the particle-particle interaction [295] . The 
importance of including proton-neutron pairing in calculations of GT excitations has also 
been illustrated in relativistic QRPA studies [272] . Pairing correlations have recently 
also been included in the proton-neutron continuum QRPA based on a phenomenological 
mean-field potential and the isovector part of the Landau-Migdal ph interaction [302J. 
In addition to pairing, deformation plays an important role in the description of GT 1 * 1 
transitions in many nuclei [2661 1303j . The sensitivity of GT strength distributions and 
double /3-decay matrix elements to the deformed mean-field has recently been analyzed 
with the QRPA |304] , and it has been shown that nuclear deformation could result in a 
suppression of two-neutrino double /3-decay rates. 

The GT resonance represents a coherent superposition of high-lying J n = 1 + 
proton-particle - neutron-hole configurations. In the relativistic formalism the GT 
operator reads 

A 

Q% = Sr ± , (97) 
i=i 

where S is defined in Eq. ( 1941) . In Fig. [42] we display the GT~ strength distributions in 
the magic nuclei 48 Ca, 90 Zr, and 208 Pb, calculated in the PN-RRPA with the DD-ME2 
effective interaction and, for comparison, data on the excitation energies of the GTR are 
also included in the figure. In addition to the high-energy GT resonance - a collective 
superposition of direct spin-flip (J — i+| — > j = l — transitions, the response functions 
display a concentration of strength in the low-energy tail. The low-lying GT excitations 
correspond to core-polarization (J — I ± | — > j — I ± |), and back spin-flip ( j = I — | 
— ► j = I + |) neutron- hole - proton-particle transitions. The strength parameter of 
the zero-range Landau-Migdal force Eq. f[9"5l) has been adjusted to reproduce the GTR 
excitation energy in 208 Pb (g' = 0.52), but we notice a very good agreement with data 
also for 48 Ca and 90 Zr. The adjusted value of g' in general depends on the choice of 
the effective interaction. By employing a set of RMF effective interactions with density- 
dependent meson-nucleon couplings [35], in Ref. [272] it has been shown that there is 
a linear correlation between the value of the nuclear asymmetry energy at saturation 
d4, and the value of g' adjusted to reproduce the GTR excitation energies: effective 
interactions with higher values of 04 require higher values of g' . 

As an example of Gamow- Teller resonances in open shell nuclei, in Fig. [43] we 
plot the PN-RQRPA strength distributions for the 112 > 118 > 124 > 130 Sn target nuclei, in 
comparison with the data for the centroids of the high-energy direct spin-flip strength, 
obtained in the Sn( 3 He,t)Sb charge-exchange reactions [271] . Direct spin-flip transitions 
dominate the high-energy region above 10 MeV. The low-energy tail of the strength 
distribution corresponds to core-polarization, and back spin-flip transitions. The solid 
curves have been calculated without including the T = proton-neutron pairing in 
the RQRPA residual interaction. The resulting high-energy GT strength in 118 Sn and 
124 Sn is divided into two main components, because of the splitting between different ph 
configurations. GTR configuration splitting (an appearance of two or more collective 
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Figure 42. Gamow- Teller strength distributions for 48 Ca, 90 Zr, and 208 Pb. PN- 
RRPA results are shown in comparison with experimental data (arrows) for the GTR 
excitation energies in 48 Ca [505] . 90 Zr [3061 HE! , and 208 Pb [TM 13051 1300] . 
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bumps with comparable intensities in the GTR strength function) was investigated in 
Ref. [310] in the framework of the shell optical model. For Sn nuclei, in particular, this 
effect was predicted to occur as the valence neutron start to occupy the level with the 
highest j in the shell: h n / 2 . The configuration splitting was attributed to the fact that 
the unperturbed energies of the (lgy 2 )(lgg^ 2 )~ 1 and (l/ig^ 2 ) (l^ii/ 2 ) — 1 configurations 
are almost degenerate. The residual interaction removes this degeneracy and, as a 
result, the main GT component separates in two distinct peaks. The ground-state 
pairing correlations have a strong influence on the occupation of the lh^^ level, and 
therefore the energy spacing between the two peaks will depend on T = 1 pairing. For 
118 Sn the calculated energy splitting of the GTR was 2.6 MeV |310j . Subsequently, the 
fragmentation and splitting of the GTR in Sn nuclei was experimentally investigated in 
Ref. |271j . The theoretically predicted configuration splitting of the main GT component 
could not be observed, however, because the total widths of the resonances ^5 — 6 MeV 
exceed the predicted splitting. The splitting of the main GT component shown in Fig. H3] 
(solid line) of ~ 3 MeV is in agreement with the result of Ref. [310] . 

The other curves shown in Fig. H3] (dotted, dashed and dot-dashed) have been 
calculated for different values of the strength-parameter of the T = proton-neutron 
pairing interaction in Eq. (196]) : Vn = 200, 250, and 300 MeV. In Ref. [31] the overall 
strength parameter V of the interaction was adjusted to the measured half-lives of 
neutron-rich nuclei in regions where the r-process nucleosynthesis path comes closest to 
the valley of stability: V = 230 MeV near N = 50, and V = 170 MeV in the N = 82 
region. In our illustrative calculation of Sn nuclei the inclusion of the T = pairing has 
a strong influence on the low-energy tail of the GT distribution in 112 < 118 ' 124 Sn, and the 
configuration splitting between the two high-energy peaks in 118 Sn and 124 Sn disappears. 
This happens because the T = pairing interaction does not affect configurations 
based on the (lg%/ 2 ) or bital (fully occupied), whereas it lowers configurations based on 
(1/1^/2) and (2g?5/ 2 ) (partially occupied). This calculation therefore demonstrates that 
the T = proton-neutron pairing strongly reduces the predicted configuration splitting 
of the main high-energy GT component. In addition to the main GTR which decreases 
in energy with increasing mass number, part of the strength associated with direct spin- 
flip transitions is concentrated at ~ 10 MeV in 112 Sn and 118 Sn. In Ref. |272j it has 
been shown that the centroid of the GT strength composed of direct spin-flip transitions 
practically does not depend on the strength of the T = proton-neutron pp residual 
interaction. 

The structure and evolution of the low-energy tail of the GT distribution determines 
the /3-decay rates of very neutron-rich nuclei, and thus sets the time scale of the r-process 
nucleosynthesis. Since the vast majority of nuclides which lie on the path of the r-process 
are out of experimental reach, nuclear structure models must be developed that can 
provide predictions of weak-interaction rates of thousands of nuclei with large neutron 
to proton asymmetry. Two microscopic approaches can be employed in large-scale 
calculations of /3-decay rates: the interacting shell model and the QRPA. The advantage 
of using the shell model is the ability to take into account the detailed structure of the 
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/3-strength function |311j . whereas the QRPA approach is based on global effective 
interactions and provides a systematic description of /3-decay properties of arbitrarily 
heavy nuclei along the r-process path. In a recent review of modern QRPA calculations 
of /3-decay rates for astrophysical applications |312j . the importance of performing 
calculations based on self-consistent mean-field models has been emphasized, rather 
than on empirical mean-field potentials, e.g. the Woods-Saxon potential. In a self- 
consistent framework both the nuclear ground states, i.e. the masses which determine 
the possible r-process path, and the corresponding /3-decay properties are calculated 
from the same energy density functional or effective nuclear interaction. This approach 
ensures the consistency of the nuclear structure input for astrophysical modeling, and 
allows reliable extrapolations of the nuclear spin-isospin response to regions of very 
neutron-rich nuclei. 

6.4- Effect of the Dirac Sea on the Gamow-Teller Sum Rule 

In nuclei all over the periodic table the GT strength distribution, when measured in the 
excitation energy region where the most pronounced GT peaks occur, is quenched by 
more than 20% with respect to the model independent Ikeda sum rule |279] 

(S^-S^)=3(N-Z), (98) 

where Sp± denotes the total sum of Gamow-Teller strength for the (3^ transition. In 
the early (p,n) studies, for instance, only « 60% of the sum rule was observed [313j . 
Two physically different mechanisms had been suggested as a possible explanation of 
the quenching of the total GTR strength: (i) nuclear configuration mixing - the high- 
lying 2p — 2h states mix with the lp — lh GT states and shift the GT strength to 
high-energy region far beyond the resonance |314[ I315[ 1316] ; (ii) excitation of a nucleon 
into the high-energy A-isobar [317[ 1318] , with the A-isobar - nucleon-hole configurations 
(A — h) coupling to the GT mode and removing part of the strength from the low-lying 
excitation spectrum |319L 1320] . In more recent (p, n) scattering experiments data on GT 
strength in the high energy region up to 50 MeV became available |307l 1321] . It has been 
shown that the measured GT strength exhausts 88% and 84% of Ikeda sum rule in 90 Nb 
and 27 Si, respectively. Most of the GT strength missing in early experimental studies 
is, therefore, recovered in the energy region where the multiconfiguration spreading 
mechanism is effective, and only a small fraction of the GT quenching may have its 
origin in A — h transitions lying high above the ordinary ph excitations. 

When the Gamow-Teller strength is calculated in the relativistic RPA framework, 
the total GT strength in the nucleon sector is reduced by ~ 12% in nuclear matter, and 
by « 6% in finite nuclei, as compared to the Ikeda sum rule [2721 12731 12891 13221 13231 1324] . 
This reduction has been attributed to the effect of Dirac sea negative-energy states, 
i.e. the missing part of the sum rule is taken by particle-hole excitations formed from 
ground-state configurations of occupied states in the Fermi sea and empty negative- 
energy states in the Dirac sea. The effect is illustrated in Fig. HH where we display the 
running sum of GTR strength relative to the Ikeda sum rule for 90 Zr, 132 Sn, and 208 Pb, 
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Figure 44. The running sum of the GTR strength for 90 Zr, 132 Sn, and 208 Pb. 
The dashed lines corresponds to PN-RRPA calculations with only positive-energy ph 
configurations. For the calculation denoted by the solid lines the RRPA space includes 
configurations formed from occupied states in the Fermi sea and empty negative-energy 
states in the Dirac sea. The total sum of the GT strength is compared to the model 
independent fkeda sum rule (dotted lines). 



evaluated in the PN-RRPA with the DD-ME2 effective interaction. The horizontal 
dotted lines denote the value 3(N — Z) of the Ikeda sum rule. The solid and dashed 
lines correspond to the values of the GTR sum calculated from — oo to the excitation 
energy denoted on the abscissa. The big jump in the calculated GTR sum occurs, of 
course, when the main GTR peak is included. The PN-RRPA calculation represented 
by the dashed lines includes only positive energy ph configurations. Even by extending 
the sum up to 80 MeV, the total sum is reduced ~ 7-8% with respect to the Ikeda sum 
rule. The Ikeda sum rule is completely exhausted by the calculated GT strength only 
when the relativistic RPA/QRPA space includes ph excitations formed from ground- 
state configurations of the fully or partially occupied states of positive energy, and the 
empty negative-energy states from the Dirac sea (solid lines in Fig. 

In Fig. H5] we plot the discrete GT spectrum for 132 Sn. Two regions of excitation 
energies are shown. The panel on the right contains the positive energy irp — vh strength, 
with a pronounced Gamow- Teller resonance peak. The panel on the left displays the 
negative energy spectrum built from na — vh transitions (a denotes a negative energy 
state). Even though these transitions are much weaker than the GTR (notice that the 
vertical scales are different for the two panels), there are many of them and their overall 
sum represents the strength missing in the positive energy sector. 
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Figure 45. The PN-RQRPA strength distribution of discrete GT _ states for 132 Sn. 
The panel on the right-hand side contains the positive energy np — vh strength. The 
panel on the left displays the negative energy spectrum built from transitions to the 
empty states in the Dirac sea. 



6. 5. Spin-Isospin Resonances and the Neutron Skin of Nuclei 

In medium-heavy and heavy neutron-rich nuclei the ground-states are often 
characterized by an extended neutron density distribution, and in some cases evidence 
has been reported for the formation of a neutron skin on the surface of the nucleus. 
The determination of neutron density distributions provides not only basic nuclear 
structure information, but it also places important additional constraints on effective 
interactions used in nuclear models. Extremely accurate data on charge densities, and 
therefore on proton distributions in nuclei, have been obtained from elastic scattering 
of electrons |325j . Data of comparable precision on neutron density distributions are, 
however, not yet available. It is much more difficult to measure the distribution of 
neutrons, though experimental information on the differences between radii of the 
neutron and proton density distributions has been reported |326l 13271 1328] . Various 
experimental methods have been used, or suggested, for the determination of the neutron 
density in nuclei |329] , but no existing measurement of neutron densities or radii has an 
accuracy better than a few percent. 

One of the modern approaches that provides information about the neutron skin in 
nuclei is based on studies of giant resonances. In particular, in Ref. (326] excitations of 
the giant dipole resonance (GDR) were analyzed, and the spin-dipole resonance (SDR) 
was studied in |327j . The GDR cross section strongly depends on the difference between 
neutron and proton density distributions |330] 1331] I332j . In Ref. [327] it has been 



demonstrated that there is a predictable correlation between the SDR cross section and 
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Figure 46. RHB plus proton-neutron RQRPA results for the energy spacings 
between the Gamow- Teller resonances and the respective isobaric analog resonances 
for the sequence of even-even 112_124 Sn target nuclei. The experimental data are from 
Ref. [271] . 

the difference between the rms radii of the neutron and proton density distributions. 
By normalizing the results to 120 Sn, data on neutron-skin thickness along the stable Sn 
isotopic chain were obtained, in good agreement with theoretical predictions. 

Recently a new method has been suggested for determining the difference between 
the radii of the neutron and proton density distributions along an isotopic chain, based 
on measurement of the excitation energies of the Gamow- Teller resonances relative to the 
isobaric analog states [333]. The Gamow- Teller resonance (GTR) J T = 1 + represents 
a collective spin-isospin oscillation with the excess neutrons coherently changing the 
direction of their spins and isospins without changing their orbital motion. The 
simplest charge-exchange excitation mode, however, does not require the spin-flip and 
corresponds to the well known isobaric analog state (IAS) J n = + . The spin-isospin 
characteristics of the GTR and the IAS are related through the Wigner supermultiplet 
scheme. The Wigner SU(4) symmetry implies the degeneracy of the GTR and IAS, 
and furthermore the resonances would completely exhaust the corresponding sum rules. 
The Wigner SU(4) symmetry is, however, broken by the spin-orbit term of the effective 
nuclear potential. The energy difference between the GTR and the IAS decreases with 
increasing asymmetry (N — Z)/A. It is implicit, therefore, that the energy difference 
between the GTR and the IAS reflects the magnitude of the effective spin-orbit potential. 
A number of relativistic mean-field calculations have shown that the magnitude of the 
spin-orbit potential is considerably reduced in neutron- rich nuclei [334] . and this is 
reflected in the larger spatial extension of the neutron density, which becomes very 
diffuse on the surface. The neutron-skin increases correspondingly. The energy spacings 
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Figure 47. The PN- RQRPA and experimental [271j differences between the excitation 
energies of the GTR and IAR, as a function of the calculated differences between the 
rms radii of the neutron and proton density distributions of even-even Sn isotopes 
(upper panel). The calculated differences r n — r p are compared with experimental 
data |327| (lower panel). 



between neutron spin-orbit partner states decrease, and this reduction is quantitatively 
in accordance with the gradual weakening of the spin-orbit term of the effective potential. 

There is a direct connection between the increase of the neutron-skin thickness in 
neutron-rich nuclei, and the decrease of the energy difference between the GTR and 
the IAS. In Fig. H6] we display the calculated differences between the centroids of the 
direct spin-flip GT strength and the respective isobaric analog states for the sequence 
of even-even Sn target nuclei. For A = 112 - 124 the results of RHB plus PN- RQRPA 
calculation (DD-ME2 density-dependent effective interaction, Gogny T = 1 pairing, 
T = pairing interaction Eq. (|96|) with Vq = 250 MeV, the Landau-Migdal parameter 
g' = 0.52 adjusted to reproduce the excitation energy of the GT resonance in 208 Pb), are 
compared with the experimental values obtained in a systematic study of the ( 3 He,t) 
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charge-exchange reaction over the entire range of stable Sn isotopes [271] . The calculated 
energy spacings are in very good agreement with the data, although for the lighter Sn 
isotopes it appears that the calculated values differ somewhat from the experimental 
trend. However, the theoretical energy spacings might depend on the details of the 
effective interaction and, in fact the data shown in Fig. HS1 provide valuable information 
that can be used to constrain the spin-isospin channel of the effective interaction. 

In Fig. H7] the calculated and experimental energy spacings between the GTR and 
IAS are plotted as a function of the calculated differences between the rms radii of 
the neutron and proton density distributions of even-even Sn isotopes (upper panel). 
We note the remarkable uniform dependence of the energy difference between the GTR 
and IAS on the size of the neutron-skin. This means that, in principle, the value 
of r n — r p can be directly determined from the theoretical curve for a given value of 
Eqtr — Eias- In the lower panel the calculated differences between neutron and proton 
rms radii are compared with available experimental data |327j . The agreement between 
the theoretical and experimental values suggests that the neutron-skin thickness can be 
determined from the measurement of the excitation energies of the GTR relative to IAS. 
This method is, of course, not completely model independent, but it does not require 
additional assumptions. Since the neutron-skin thickness is determined in an indirect 
way from the measurement of the GTR and IAS excitation energies in a sequence of 
isotopes, in practical applications at least one point on the theoretical curve should be 
checked against independent data on r n — r p . 

7. Exotic Nuclear Modes in Astrophysical Processes 

Collective modes play an important role in astrophysical processes that involve both 
stable and exotic nuclei far from stability. In particular, spin-isospin excitations such 
as the Gamow- Teller resonance are essential for weak-interaction processes, e.g. (3- 
decay of exotic nuclei on the path of r-process nucleosynthesis, electron capture on 
neutron-rich nuclei at temperatures and densities characteristic for stellar core collapse, 
neutrino-induced reactions on heavy neutron-rich nuclei in the post-collapse supernova 
environment, including the process of neutrino nucleosynthesis. A recent review of 
nuclear weak- interaction processes in stars can be found in Ref. }311j . 

In this section we will review the role of non-charge-exchange nuclear collective 
modes in astrophysical processes. Particularly interesting is the effect of the low-lying 
dipole transition strength on the r-process nucleosynthesis, and in the propagation of 
ultra high-energy cosmic rays. We will also discuss exotic collective modes in the crust 
of a neutron star. 

7.1. Low-Energy Dipole Strength and the r-process 

Approximately half of the nuclides with A > 60 found in nature are formed in the 
rapid neutron-capture process (r-process) nucleosynthesis. The nuclear input for r- 
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process calculations necessitates the knowledge of the properties of thousands of nuclei 
far from stability, including the characteristics of strong, electromagnetic and weak 
interaction processes [335J. Most of these nuclei are not accessible in experiments 
and, therefore, many nuclear astrophysics calculations crucially depend on accurate 
theoretical predictions for the nuclear masses, bulk properties, nuclear excitations, (n, 7) 
and (7, n) rates, a- and /5-decay half-lives, fission probabilities, electron and neutrino 
capture rates, etc. Early calculations of the nuclear processes relevant for astrophysical 
applications were based only on phenomenological models |336J I337[ 13381 339J. Only 
more recently large-scale microscopic calculations became standard in the prediction 
of nuclear masses [340] . and dipole strength distributions [1061 1107] . The availability 
of large scale microscopic calculations opens the possibility for global predictions of 
the nuclear ingredients for the r-process, based on an universal nuclear energy density 
functional. However, fully microscopic calculations of nuclear observables over the whole 
isotope chart are not yet feasible, especially for excited states. Only the excitation 
spectra of all even-even spherical nuclei can be obtained from fully consistent microscopic 
Skyrme-QRPA calculations [29], whereas for deformed, odd-even, and odd-odd nuclei 
a series of approximations must be designed specifically for astrophysical applications, 
and implemented in the microscopic QRPA. 

Two principal candidates have been suggested for the astrophysical site of the r- 
process. In the first scenario the r-process takes place in explosive stellar events, such 
as the core collapse supernova [335] , in an environment characterized by high neutron 
densities (N n ~ 10 20 cm -3 ), so that successive neutron-captures proceed into regions 
of neutron-rich nuclei far away from the valley of /3-stability. If the temperature or 
the neutron density that characterize the r-process are low enough to break the (n,7) 
- (7,n) equilibrium, the waiting-point approximation is not valid any more, and the r- 
abundance distribution directly depends on the neutron-capture rates on exotic neutron- 
rich nuclei [336J. 

Recently an alternative scenario for the r-process has attracted renewed interest. 
It is related to the decompression of cold neutron star matter, in particular its crust 
(see [335j for a more detailed description). In this scheme the production of heavy 
nuclei follows a completely different path from the core-collapse supernovae scenario. 
The inner crust of a neutron star is composed of nuclear clusters immersed in a neutron 
gas. When a decompression of this crust occurs, nuclear clusters and the neutron gas 
are both ejected, and this leads to a decrease of the matter density. The /3-equilibrium 
is broken and nuclei with Z in the range between 40 and 70 are produced. At very 
low density drip-line nuclei are formed, immersed in a neutron flux of N n ~ 10 35 cm -3 . 
The competition between the neutron-capture and the /3-decay, starting from these 
drip-line nuclei, leads to the production of heavy nuclear systems. Because this type 
of nucleosynthesis starts from drip-line nuclei, where the waiting-point approximation 
cannot be applied, theoretical predictions of various structure phenomena are essential 
ingredients for r-process modeling. 

Several types of nuclear observables are therefore required for the description of 
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r-process abundances. In addition to the /3-decay rates, it is also necessary to describe 
the (n,7) rates. This process can be divided into two steps: the neutron capture and 
the photo-deexcitation |341j . Nuclear masses and neutron-nucleus optical potentials 
enter into the calculation of the neutron-capture rates. The description of the photo- 
deexcitation process necessitates predictions of the El strength functions, as well the 
level densities in daughter nuclei. Neutron capture rates are evaluated in the framework 
of the Hauser-Feshbach statistical model, which is based on the fundamental assumption 
that the capture process occurs with the intermediary formation of a compound nucleus 
in thermodynamic equilibrium. In this approach the Maxwellian-averaged (11,7) rate, 
at temperatures characteristic for the r-process environment, strongly depends on the 
electromagnetic interaction, i.e. on the photo-deexcitation probability. The modeling of 
the r-process abundances requires a reliable extrapolation of the El-strength functions 
towards the neutron-drip line. Figure HE] shows a schematic picture of the (n,7) reaction 
in the statistical model framework. 




T E1 (E) p(E) dE 



(Z,A+1) 



Figure 48. Illustration of the statistical radiative (n.,7) neutron-capture. T„ is 
the neutron-capture coefficient, and T 7 is the photo-transmission coefficient. T 7 is 
determined by the El strength Tei(E), and the level density p(E). 

Low-energy dipole modes are expected to play a more important role in the r- 
process than the higher lying strength, including the giant dipole resonance (GDR). 
Namely, the low-lying modes located close to neutron separation energy S n are directly 
sustained by neutron-capture reactions, as shown in Fig. HH1 The importance of the 
soft dipole modes has been emphasized in Ref. [336] . where it has been shown that the 
presence of a low-lying resonant component of the El strength leads to an increase of 
the radiative neutron-capture rate by factors 10 - 100, for nuclei with S n between 2 MeV 
and 4 MeV. The r-abundance distribution is affected because the presence of low-energy 
dipole resonances accelerates neutron capture and allows the production of heavy nuclei 
around A = 130. 

Large-scale calculations of El-strength functions for astrophysical applications are 
usually performed using phenomenological Lorentzian models |336] . Several refinements 
can be introduced, such as including the energy dependence and/or the temperature 
dependence of the width of the Lorentzian [336[ I337[ I338[ 1339] . The Lorentzian GDR 
approach presents, however, several problems. On one hand, in this framework it is not 
possible to predict the enhancement of the El strength at energies close the neutron 
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separation energy. On the other hand, even if a Lorentzian-type function provides a 
suitable description of the El strength in stable nuclei, the location of its maximum and 
its width for each nucleus remain to be predicted from some systematics or underlying 
model. For astrophysical applications these properties have often been obtained from a 
droplet-type model [342] . This approach is clearly not reliable when dealing with exotic 
nuclei, and this was already demonstrated in Refs. [106[ 1112] . In order to achieve a 
better description of the r-process, one has to improve the nuclear structure modeling. 
Generally speaking, the more microscopic the underlying theory, the more reliable will 
be the extrapolations towards the neutron drip-line. Microscopic calculations of the El- 
strength functions for the whole nuclear chart have recently been reported [106L 1107] . In 
a first step the dipole response was calculated with the QRPA based on the HF-Skyrme 
plus BCS description of nuclear ground states |106j . In neutron- rich nuclei pronounced 
El strength was predicted in the low-energy region below the giant dipole resonance. 

The dipole strength can also be calculated with the QRPA based on HFB ground 
states [29l 1343] . In this microscopic approach photoabsorption cross sections were 
determined |107j . and one was able to judge the ability of the different forces to reproduce 
experimental data. In Refs. [3401 13441 1345[ 1346] nuclear deformation was also taken into 
account. Based on the Skyrme-HFB approach, a number of new effective forces have 
recently been introduced [3401 1345[ 1346] . Their parameters were exclusively adjusted 
to the 2135 experimental masses [347] . with some additional constraints related to 
the stability of neutron matter and the incompressibility of nuclear matter. The new 
effective interactions BSk2-BSk7 are summarized in Ref. [335J. 

Photo- induced reaction cross sections were compiled in Refs. [3481 1349] . and 
represent the most reliable source of data with which HFB+QRPA predictions can 
be compared. These include the GDR parameters (the peak energy, peak cross section, 
and the full width at half maximum) observed in photonuclear reactions measured by 
bremsstrahlung, quasimonoenergetic and tagged photons for 84 nuclei. Among these, 48 
nuclei are spherical and can be used to test the (HF-BCS or HFB)+QRPA predictions. 
Figure 09] shows the calculated GDR centroid energies of these 48 spherical nuclei using 
the HF-BCS+QRPA with the SLy4 parameterization, in comparison with the data. 

For the BSk6 and BSk7 interactions, in particular, the comparison shows that these 
functionals not only reproduce the experimental masses with great accuracy (the rms 
deviation is only 0.676 MeV for the set of 2135 known masses), but also are suited for 
the description of El collective excitations [107J. The agreement with data is of the 
same level of accuracy as in the case of phenomenological models, i.e. the rms deviation 
for the GDR excitation energies is around 500 keV. However, since the microscopic 
HFB+QRPA approach potentially encompasses a wide range of phenomena, including 
the soft dipole mode, it is clearly preferable in large-scale calculations. 

Turning now to astrophysical applications, microscopic predictions need to be 
extended on the whole nuclear chart, including also nuclei with an odd number of protons 
and/or neutrons, and deformed systems. The approximations that are currently used 
in these calculations include: the filling approximation for odd-A and odd-odd nuclei, 
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Figure 49. Comparison between the experimental GDR energies and the HF- 
BCS+QRPA calculations with the SLy4 force [106] . 



a phenomenological damping of the El strength, and a phenomenological treatment 
of deformation. These estimates may be improved in the near future, for instance by 
developing a deformed QRPA model that could be used in large-scale calculations of 
thousands of nuclei. 

The QRPA provides a reliable description of the GDR centroid energy, but it 
is necessary to go beyond this approximation to describe the damping of collective 
motion. The GDR are empirically known to have rather large widths and therefore 
finite lifetimes, which can be described by several models |4"3~t |92 | [350] . In Ref. [106] . for 
instance, the QRPA strength function was folded by an arbitrary Lorentzian adjusted 
to the empirical GDR width. The damping of the El strength can also be described 
by the approximate procedure developed in Ref. [43]. In the large-scale calculation of 
Ref. [107] the QRPA strength was folded by a Lorentzian function representing the self- 
energy operator [43j 1351] . It is also necessary to introduce a temperature-dependent 
correction factor in the expression for the GDR width [3381 13521 1353] . In deformed 
axially symmetric nuclei the GDR splits into two major components as a result of the 
different resonance conditions characterizing the oscillations of protons against neutrons 
along the axis of rotational symmetry, and an arbitrary axis perpendicular to it. In the 
phenomenological approach, the Lorentzian-type formula is generalized to a sum of two 
Lorentzian functions [354J. 

Large-scale QRPA calculations based on the BSk7 Skyrme interaction have recently 
been performed for some 8300 nuclei with 8 < Z < 110, extending between the 
proton and neutron drip-lines [107] . In the region of neutron-deficient nuclei, as well 
as along the valley of /^-stability, the calculated El-strength functions are very similar 
to the empirical Lorentzian approximation. However, in neutron-rich nuclei the QRPA 
predictions start to deviate from the simple Lorentzian shape. In particular, low-lying 
transitions are found at excitation energies well below the GDR, and their strength 
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increases with neutron excess. This effect has been discussed in detail in Sec. 14.41 In 
Fig. [50] we plot the El strength function in Sn isotopes, calculated in the HFB+QRPA 
model with the BSk7 effective interaction. For A > 140 a significant portion of the 
strength is concentrated at low energies: E ~ 5 — 7 MeV. It should be noted that 
phenomenological models cannot predict these low-energy components. For 150 Sn, for 
instance, all phenomenological systematics that are used in calculations of neutron- 
capture cross sections, predict a 7-ray strength peaked around 15 MeV with a full 
width at half maximum of about 4.5 MeV [355] . This is obviously very different from 
the microscopic prediction shown in Fig. [50] More generally, HFB+QRPA calculations 
confirm that the neutron excess affects the spreading of the isovector dipole strength, as 
well as the centroid of the strength function. The energy shift is larger than predicted 
by the usual A~ 1//6 or A^ 1 ^ 3 dependence given by the phenomenological liquid drop 
approximations |342] . The basic features of the QRPA El strength function for nuclei 
with a large neutron excess are qualitatively independent of the choice of the effective 
interaction. 



20 
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Figure 50. El strength functions for the Sn isotopic chain, calculated in the 
HFB+QRPA model with the BSk7 force. Only isotopes between A=115 and A=150, 
with a step of A A = 5 are shown [107j . 

The radiative neutron capture cross sections are calculated with the Hauser- 
Feshbach statistical model [356J, starting from nuclear ground state properties 
determined consistently in the microscopic HFB model with the same BSk7 Skyrme 
force [340] . The calculation includes the improved nuclear level density prescription 
based on the microscopic statistical model, also used to estimate the nuclear 
temperature |357] . The direct-capture contribution, as well as the possible overestimate 
of the statistical predictions for resonance-deficient nuclei, could have a significant effect 
on the radiative neutron capture by exotic nuclei [336] ■ However, we will focus on the 
role of the dipole strength, which is almost exclusively probed by the statistical model. 

The results of the Hauser-Feshbach calculation can be compared to the 
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experimental (11,7) cross-sections. However, it is more convenient to measure the 
inverse (7,11) reaction. Such an experiment has recently been performed on 181 Ta: 
181 Ta(7,n) 180 Ta [358] . In FigioTl we show the comparison between the data and the 
Hauser-Feshbach calculations, obtained using either the QRPA El strength, or the 
phenomenological strength distributions. We notice that at low energies (E < 10 
MeV) the microscopic calculation of the dipole strength produces results which are 
in excellent agreement with data, but significantly different from those obtained with 
the phenomenological approach. The tail of the cross section between 7.5 and 10 MeV 
is attributed to the presence of the low-energy dipole (pygmy) mode in the QRPA 
calculation, which does not appear in the phenomenological models. 




E [MeV] 

Figure 51. Comparison of the experimental photoneutron cross section for 181 Ta, 
with the Hauser-Feshbach predictions obtained with the QRPA, hybrid, and Lorentzian 
models |358j . 

The Maxwellian-averaged neutron-capture rates calculated with the HFB+QRPA 
El strength functions are compared in Fig. [52] with those based on the Hybrid 
phenomenological formula, for all nuclei with 8 < Z < 110. The Hybrid El strength 
differs from the QRPA estimate in the location of the centroid energy, as well as in 
the low-energy tail. Obviously the El strength obtained with QRPA enhances the 
capture rates by a factor up to 10 close to the neutron drip-line. For r-process nuclei 
characterized by neutron separation energies S n < 3 MeV, neutron capture proceeds 
much faster than predicted by the phenomenological Hybrid formula. This is due to the 
shift of the GDR to lower energies as compared to the usually adopted liquid-drop A^ 1 ^ 
rule, as well as to the appearance of dipole modes at low energies. Both effects tend to 
enhance the El strength at energies below the GDR, i.e in the energy window relevant 
for the neutron capture process. For less exotic nuclei this effect is much smaller, and the 
differences are mainly due to the predicted location of the GDR and the strength of the 
low-energy tail. When compared to the HF-BCS+QRPA results [106], the HFB+QRPA 
model [10TJ predicts very similar neutron-capture rates, even close to the neutron drip- 
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line (see the lower panel in Fig. IB^l) . This demonstrates the consistency of the results 
obtained with various microscopic approaches. 




Figure 52. Upper panel: ratio of the Maxwellian-averaged (n.,7) rates (at a 
temperature of 1.5 10 9 K) obtained with the HFB+QRPA El strength to those 
calculated with the Hybrid formula [336] . as a function of the neutron separation 
energy S n for all nuclei with 8 < Z < 110. Lower panel: the ratios of the 
HFB+QRPA neutron-capture rates and those calculated with the HF-BCS+QRPA 
model of Ref. [106] . 



Many further improvements may be useful, but this will require intensive theoretical 
and computational advances. For instance, large-scale microscopic QRPA predictions 
of the El strength in deformed nuclei, the inclusion of the particle-vibration coupling 
effects on the low-energy strength, an improved treatment of odd and odd-odd nuclei, 
etc. The aim is to achieve a fully microscopic description of the r-process based on an 
universal nuclear energy density functional. 



7.2. Nuclei as Ultra-High Energy Cosmic Rays 

Cosmic rays are energetic particles that originate in the Universe, with observed energies 
up to ~ 3 10 20 eV |359[I360| . Astrophysical sites able to accelerate particles to such ultra 
high energies are currently under discussion. Among them, violent processes related to 
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neutron stars are possible candidates, because they generate strong magnetic fields able 
to confine protons of energy 10 20 eV. Figure [53] illustrates various natural accelerators in 
the Universe, and displays their typical magnetic field B with respect to their curvature 
radius p. The product Bp value gives the energy which can be reached by the accelerated 
particles, from the well known relation Bp = p/Q where p is the momentum of the 
particle, and Q its charge. Sites which could possibly accelerate protons to 10 20 eV are 
located in the shaded band. 
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Figure 53. Various accelerators in the Universe, located on the plot of their magnetic 
field versus their size. 

The composition of Ultra High Energy Cosmic Rays (UHECR) is not known. 
According to two principal hypotheses, they are composed either only of protons, 
or a mixture of protons and nuclei, ranging from hydrogen to iron. It has been 
known for almost four decades that UHECR interact with the 2.7 K cosmic microwave 
background radiation (CMB), leading for instance to a spectacular decrease in their 
flux above energies around 10 20 eV - the so-called Greisen, Zatsepin and Kuzmin (GZK) 
suppression |361[ I362] . The following UHECR issues need to be addressed: the initial 
composition, the acceleration mechanism, and the intergalactic propagation. Going 
through these steps, the UHECR may reach an Earth based detection setup, such as 
the Auger detector [363] . The first two issues include many unsolved questions about the 
composition and acceleration processes of UHECR. It is therefore necessary to describe 
very accurately their propagation in order to provide a test of the composition and 
acceleration scenarios by comparing the predictions with data measured on Earth. 

In the rest frame of a nucleus, at typical UHECR energies of 10 19 -10 21 eV, the CMB 
photons are boosted to the energy range between a few hundred keV, up to a few hundred 
MeV. The interaction process between the UHECR and the CMB is dominated by the 
giant dipole resonance (GDR) at photon energies below 30-50 MeV, and to a lesser 
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extent by the quasideuteron emission at intermediate energies (between 50 MeV and 
150 MeV), and the pion photoproduction at energies above 150 MeV [3641 [365] . Figure 
1541 displays the mean free path of 56 Fe nuclei for various photonuclear processes, and 
shows the predominant role of the GDR. Nuclei photodesintegrate by emitting nucleons 
through (7,n), (7,p), (7,2n), ... reactions. It is therefore necessary to accurately describe 
the dipole strength for the nuclei on the photodisintegration path from Fe to protons. 
It should be noted that many nuclei along this path are unstable, for instance 44 Ti. 

In order to describe the changes in the abundance of heavy nuclei as a result of 
the interaction of the UHECR with the CMB, a nuclear reaction network that includes 
all interactions of interest must be used. The chosen set of nuclear species are coupled 
by a system of differential equations corresponding to all the reactions affecting each 
nucleus, i.e. mainly photodisintegrations and /3-decays |366j . All nuclei lighter than the 
seed nuclei and located between the valley of stability and the proton drip-line must be 
included in the network. Under the most natural astrophysical assumptions, UHECR 
are accelerated out of the ambient gas, possibly enriched in Fe close to neutron stars or 
depleted in metals (i.e. nuclei heavier than H) if significant photodisintegration occurs 
during the acceleration stage itself. Therefore, if nuclei are indeed present among the 
UHECR, it is expected that they typically include the most abundant elements found 
in the interstellar medium, i.e. essentially those lighter than Fe. The interaction of 
UHECR with the CMB is thus expected to include all possible nuclei resulting from the 
photodisintegration of the heaviest species and therefore involve all stable and neutron- 
deficient unstable isotopes with A < 56. 




Lo re ntz factor, r 



Figure 54. Mean free path of 56 Fe nuclei against the various processes in the CMB: 
Giant Dipole Resonance (GDR), Quasi-Deuteron process (QD), Baryonic Resonance 
(BR) and Pion Photoproduction (PF) [367] , 

The UHECR photodisintegration was originally investigated by Puget, Stecker and 
Bredekamp (PSB) |364[ 1368] . However, in the PSB model two major approximations 
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were employed to estimate the intergalactic UHECR propagation. The first is related 
to the total photoabsorption cross section which is parameterized as a simple Gaussian 
function |364[ I3~68j. abruptly cut below the theoretical reaction threshold. The second 
is based on the use of a reduced reaction network, including only one nuclide for each 
value of A, to estimate the time evolution of the UHECR composition. More precisely, 
assuming that the /3-decay of the unstable nuclei produced by photodisintegration 
is always faster than the corresponding photoemission rate, a unique nuclear path 
is followed from the initial 56 Fe source to the final protons |364] . neglecting the 
contribution of unstable nuclei. It is necessary to achieve an accurate description 
of the photodisintegration rates using photoreactions that are extensively studied 
in nucleosynthesis, where phenomenological parameterizations of the photoabsorption 
cross sections have been optimized during the last decades, and where large-scale 
microscopic predictions have also become available, as described in section I7TT1 |106l [TUT] . 
New compilations of experimental photoabsorption data also help to determine the 
degree of accuracy with which the present reaction models predict the corresponding 
cross sections. Important progress has also been made in the field of nucleosynthesis by 
solving large reaction networks on the nuclear chart and thus following the time evolution 
of the composition of the material at given astrophysical sites. Similar tools can therefore 
be used in the field of UHECR in order to take into account the contribution of unstable 
nuclei during the photodisintegration path. 

The total photon transmission coefficient characterizing the probability to populate 
by photoabsorption a compound nucleus excited state is obviously one of the key 
ingredients for the evaluation of the photoreaction rates. In the specific astrophysical 
conditions considered for UHECR energies of 10 19-21 eV, this function is dominated 
by the El transition. To estimate the accuracy of the different methods available 
for the evaluation of an ET-strength function, four models have been considered: the 
Lorentzian [369], the generalized Lorentzian [352], the HF-BCS+QRPA |106j . and the 
HFB+QRPA [107J. The former two are phenomenological, and the latter two are 
microscopic. 

The photoreaction cross sections are estimated with the Talys nuclear reaction 
code |235| 370J, which takes into account all types of direct, pre-equilibrium and 
compound mechanism to estimate the total reaction probability, as well as the 
competition between the various open channels. The quasideuteron process is neglected 
because of the limited photon energy range [365J. The predictions are compared with 
available experimental data [371] for nuclei with A<56. It should be noted that, even 
for stable nuclei, the data on such nuclei are scarce. For instance, total photoabsorption 
cross sections around the GDR peak energy are available for only 10 nuclei |371j . An 
extensive study has been performed to compare the predictions with the data [366J. Both 
the microscopic and the Lorentzian approaches correctly describe the data, because the 
nuclei involved lie close to the valley of stability. 

The intergalactic UHECR propagation is calculated considering the interaction with 
the CMB. For illustrative purposes, we will restrict ourselves in a first step to study the 
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propagation of a UHECR source made of 56 Fe only. All stable and neutron-deficient 
unstable nuclei with A<56 are included in the reaction network. Fig. [55] shows the 
evolution of the average mass number < A > as a function of the distance from the 56 Fe 
source, calculated with the four GDR prescriptions. For a given source distance, < A > 
is the average value of the calculated nuclei abundances. The full reaction network 
is solved at each time-step, taking into account all open photoemission channels, i.e 
(7,n), (7,p), (7,a), (7,2n), (7,2p), (7,2a), (7,np), (7,11a), (7,pa). In other words, the 
abundance of each nuclide is derived by taking into account the contribution of all 
production channels, from the source nucleus downwards the table of nuclides, with 
the appropriate weight determined according to the corresponding cross sections. The 
values obtained are thus equilibrium values, representing the composition which would 
result from the propagation of an infinite number of nuclei up to the time considered. 
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Figure 55. Evolution of the average mass number < A > with respect to the distance 
from the 56 Fe source for three different energies. Left: E = 3.6 x 10 20 eV; center: 
E = 10 21 eV; right: E = 2.61 x 10 21 eV [366] . 



To illustrate the role of the soft dipole mode, let us consider a 56 Fe nucleus with an 
energy of 3.6 x 10 20 eV. In this regime only the lowest energy part of the .El-strength 
overlaps with the photon density n(e). The distance of propagation is mainly sensitive to 
the low-energy component of the .El-strength function, and using different prescriptions 
leads to significant differences in the propagation distance. This emphasizes the necessity 
accurately describe the low energy tail of the El strength, even for nuclei close to the 
valley of stability. In contrast, results for 56 Fe at higher energy mainly depend on 
the location of the GDR peak or integrated photoabsorption, and for this reason the 
propagation distance is less sensitive to the photoreaction details. 

The UHECR propagation distance has been estimated using the complete reaction 
network. The initial PSB calculations were based on the reduced PSB path illustrated in 
Fig. [561 In this approximation, as explained above, only one stable isotope is considered 
for an isobaric chain and the corresponding isobars are not affected by competitive 
channels. However, as shown in Fig. [561 about 85 nuclei are involved in the 56 Fe 
photodisintegration at E = 10 21 eV, and a number open channels, including /3-decay, 
can compete (the Lorentz dilation of time allows /3-unstable nuclei with half-lives of the 
order of an hour to survive over a Mpc scale, and thus have a chance to interact with 



Exotic modes of excitation in atomic nuclei far from stability 



107 



a CMB photon). Most of the stable nuclei involved in the photodisintegration process 
have more neutrons than protons. Neutron photo-emissions are therefore favored and 
the corresponding unstable nuclei will /3 + -decay towards the valley of stability. 
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Figure 56. Nuclei involved in the photodisintegration process of 56 Fe at E = 10 21 eV. 
Unstable nuclei are denoted by shaded squares, and the PSB path is indicated 
by arrows. The mass number A of each nucleus is written in the corresponding 
square [366] . 



Significant discrepancies are therefore expected between the recent calculation in 
Ref. |366j . and the original PSB results based on the reduced path and the Gaussian 
parameterization of GDR strengths. In particular, as seen in Fig. [56] for A>45, about 
70% of the nuclei are shortcut by the simplified PSB path. For E = 3.6 x 10 20 eV, a 
significant difference is found between the PSB predictions and the one using an accurate 
El strength description, such as the QRPA or the generalized Lorentzian. This effect 
is due to the prediction of the low-energy part of the dipole strength, which in the PSB 
case does not agree with the data. 

A full propagation calculation based on a Monte-Carlo simulation, has been 
performed in Ref. |367j . A mixed source of protons and nuclei with A<56 has 
been considered, as well as quasideuteron, pion photoproduction, and pair-production 
processes. As shown in Fig. |57] the propagated spectra are in good agreement with 
the cosmic ray data. The effect of the prediction of the dipole strength compared to 
the PSB parameterization is non- negligible, and the sensitivity to the description of the 
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El strength will become even more important, because the data provided by the Auger 
detector will drastically reduce the statistical error bars at high energy. 
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Figure 57. Propagated spectra from a source composed of protons and nuclei (solid 
line), compared with available data [367] . 

In Ref. [367J it has also been shown that a mixed composition of protons and 
nuclei provides a reasonable interpretation of the high energy part of the CR spectrum 
around the so-called ankle area. Better agreement with observed data is obtained with 
composite source models, than with just a uniform proton source. Moreover, the work 
of Ref. [3 T2j has shown that by including nuclei in UHECR a much better agreement 
with complementary experimental observables is obtained. All these studies support 
the hypothesis that UHECR may be composed of nuclei. 

7. 3. Supergiant Resonances in the Inner Crust of Neutron Stars 

Some of the most exotic nuclear excitations could arise in the crust of neutron stars. 
In microscopic calculations the inner crust matter is usually described in the Wigner- 
Seitz (WS) approximation |373l I374] . i.e. the inner crust is modeled by non-interacting 
cells that contain a neutron-rich nucleus immersed in a dilute gas of neutrons and 
relativistic electrons. For baryonic densities in the range from 1.4 10~ 3 po about 0.5po ; 
where po=0.16 fm -3 is the nuclear matter saturation density, nuclear clusters can be 
considered spherical [3741 I375J. At higher densities, the inner crust matter can develop 
various non-spherical phases (e.g. rods, slabs, tubes, bubbles) |373j . The structure of 
the WS cells has been analyzed as a function of their density |374j . which is related 
to the distance of the WS cell from the center of the neutron star. The equilibrium 
condition for the WS cells has been derived, with the values of Z and N that minimize 
the energy of the system. The resulting values of Z are typically located between 10 and 
50, depending on the density of the WS cell. The corresponding number of neutrons 
N is about several hundreds. Therefore the cell with Z=50 and N=1750 is quoted as 
i800g n j n ce n j-jjg density f tjhe neutron gas far from the nuclear cluster is ~ 0.018 
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fm -3 , and the radius of the cell is 27.6 fm. The WS cell can therefore be represented 
as a drip-line nucleus immersed in a neutron gas, and considered as a system between a 
finite nucleus and the uniform neutron matter. Nuclear collective modes are expected to 
develop in such systems, and several interesting questions arise: What is the structure 
of collective states in these extremely neutron-rich systems? What is the relationship 
between the nuclear cluster and the neutron gas? Can the study of very neutron-rich 
nuclei be helpful in understanding these excitations? 

The cooling of low-mass neutron stars is strongly influenced by the superfluid 
properties of inner crust matter [376] . These properties and their effect on the 
specific heat have been analyzed in various theoretical frameworks, e.g. semiclassical 
pairing models |377j I378j . Bogoliubov-type calculations based on a Woods-Saxon mean- 
field |3791 1380] . and the self-consistent Hartree-Fock-Bogoliubov (HFB) approach [3811 
382J. These calculations have shown that pairing correlations can reduce by orders 
of magnitude the specific heat of baryonic matter in the inner crust. It must be 
emphasized, however, that the specific heat was evaluated for a system of non- 
interacting quasiparticles. This quantity can also be strongly affected by collective 
modes determined by the residual interaction between quasiparticles, especially if these 
modes appear at low-excitation energy. This effect has been studied in Ref. [383[ [384|. 
It should be noted that the specific heat of the inner crust is also determined by the 
motion of electrons and, to a lesser extent, by lattice vibrations [3731 13761 13801 [385J . 
These degrees of freedom of the inner crust matter will not be taken into account in the 
present discussion. 

In Ref. [383J RPA calculations were performed for 580 Sn, and a pronounced low- 
lying quadrupole peak was obtained. In order to take into account pairing effects, it is 
convenient to calculate the collective response with the HFB+QRPA model formulated 
in the coordinate representation [29J. This representation is particularly suited to 
describe systems with a large number of quasiparticle states, such as WS cells. In the 
first step of the calculation the HFB equations for the ground-state of the given WS cell 
are solved, considered as an isolated system. The HFB calculations are performed in the 
coordinate representation and the Dirichlet-Neumann boundary conditions are imposed 
at the contour of the cell [374] . These are the only discretization conditions which 
can produce a constant density around the outer boundary of the WS cell. Figure [58] 
displays the HFB results for the particle densities of the 1800 Sn and 982 Ge WS cells p84j. 
The density profiles are very diffuse, because of the presence of the neutron skin on the 
surface of the nuclear cluster (~ 7 fm). At larger radii the densities remain constant 
and correspond to the neutron gas component. 

The response of several cells has been calculated within the HFB+QRPA model, 
e.g. 1500 Zr, 1800 Sn and 950 Gc [384J. In these cells a very collective low-lying state has 
been predicted, located in the energy region between 2 and 4 MeV. This SuperGiant 
Resonance (SGR) typically exhausts around 70 % of the EWSR [384J. The analysis 
of the SGR structure has shown that the contributions of the nuclear cluster and the 
neutron gas to the neutron transition density are comparable. However, the magnitude 
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of the strength (but not the energy position) is mainly determined by the neutron gas, 
because it is given by J dr r L+2 5p{r) for a given multipole L (cf. Eq ( J75il ). and this 
strongly favors effects located at large radii. 

The monopole, dipole, and quadrupole neutron response of the cell 1800 Sn are shown 
in Fig. [59j We display both the unperturbed HFB response built from non-interacting 
quasiparticle states, and the QRPA response function. When the residual interaction 
between quasiparticles is turned on, the unperturbed spectrum, which is distributed 
over a large energy region, is collected in a strong peak located at about 3 MeV. 
All multipolarities exhibit the SGR, and in the case of the quadrupole response the 
SGR peak collects more than 99% of the total quadrupole strength. This mode is 
extremely collective, there are more than one hundred two-quasiparticle configurations 
contributing to the QRPA amplitude, and its reduced transition probability B(E2) 
pd 25 x 10 3 Weisskopf units. This value is two orders of magnitude larger than the 
B(E2)'s found in ordinary nuclei. It should be noted that the extrapolation of the 
energy position based on the Giant Quadrupole Resonance systematics in finite nuclei, 
i.e. 65A -1 / 3 MeV [6], predicts the low-energy peak at about 5 MeV. The additional 
lowering of the peak, which is due to the comparable contributions of the nuclear cluster 
states and the neutron gas, shows that the WS cells cannot be simply considered as giant 
nuclei. 

In order to study collective excitations at higher baryonic densities of the inner 
crust matter, the response of the cell 982 Ge has also been calculated in Ref. [384J. In 
this cell the cluster and the neutron gas are less separated than in the case of 1800 Sn 
(see Fig. [58J) • Because of the large degeneracy of states located close to the Fermi level, 



Exotic modes of excitation in atomic nuclei far from stability 



111 



> 



60 r 
50 - 
40 - 
30 - 
20 - 
10 - 




2 4 6 8 10 12 14 16 18 20 

E* (MeV) 




10 12 14 16 18 20 



x 10 

2000 

,-n 1750 
> 1500 
S 1250 

"*d 1000 
X 750 

*W 500 
W 250 




E (MeV) 



2 4 6 8 10 12 14 16 18 20 

E* (MeV) 
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distributions for the cell 1800 Sn. The solid and dashed curves represent the QRPA 
strength and the HFB unperturbed strength, respectively. 



many particle-hole configurations add coherently to form the SGR. It is also interesting 
to investigate whether the contribution of the free neutron gas to the collective mode 
could be eventually described in a semiclassical picture by using the hydrodynamic 
model. This requires a coherence length that is much smaller than the size of the 
system [386j, and is only fulfilled for the most external WS cells, such as 1800 Sn. For the 
WS cells of higher density, such as 1500 Zr and 982 Ge, the hydrodynamic model cannot 
reproduce the energy position of the SGR [384J. 

The predicted energies of the SGR are of the same order of magnitude as the 
average pairing gap of the neutron superfluid and, consequently these modes could 
have a significant effect on the entropy and the specific heat of baryonic matter in 
the inner crust. A quantitative estimate of these effects will require finite-temperature 
HFB+QRPA calculations [387]. There are, however, open questions about the validity 
of the non-interacting WS approximation in the description of collective modes. Their 
amplitude is strongly sensitive to the surface, where the neutron gas connects two 
contiguous cells. Tests have been performed by modifying the Dirichlet-Neumann 
boundary conditions, but no pronounced deviation of the strength has been found. 
The role of the coupling between two contiguous WS cells through the neutron gas 
should be investigated. The cells which have an impact on the specific heat remain also 
to be determined and, finally, neutron-rich nuclei which will be produced by the next- 
generation radiocative-beam facilities, could provide data that will constrain models and 
functionals which describe WS cells. 

8. Concluding Remarks and Outlook 

Experimental studies with radioactive isotope beams have disclosed a wealth of structure 
phenomena in nuclei far from the line of /3-stability. Among the most interesting results 
are those on the evolution of collective modes of excitation in regions of unstable nuclei, 
and the possible existence of entirely new types of excitation in weakly-bound systems. 
The best studied example is the evolvement of the dipole response in neutron-rich nuclei, 
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and the possible occurrence of the pygmy dipole resonance: the resonant oscillation 
of the weakly-bound neutron skin against the isospin saturated proton-neutron core. 
However, even in this case the available data on the low-energy dipole response do not 
discriminate between various theoretical interpretations. In nuclei near the neutron drip- 
line new types of excitation could be induced by the strong surface pairing properties 
as, for instance, the soft dipole pairing mode. On the proton-rich side, even if the 
protons do not develop a pronounced skin structure, in relatively light nuclei close to 
the proton drip-line the weakly-bound proton orbitals, including states which are bound 
only because of the presence of the Coulomb barrier, could contribute to the evolution 
of the proton pygmy dipole resonance. The existence of exotic isoscalar modes has also 
been predicted: a low-lying dipole mode which could correspond to the toroidal dipole 
resonance, giant pairing vibrations, quadrupole excitations and higher multipoles. The 
increasing interest in exotic collective modes has led to important advances both in the 
experimental methods and theoretical tools that are used in the study and interpretation 
of these phenomena. 

A modern microscopic time-dependent mean-field theory of collective excitations 
in nuclei far from stability necessitates a fully self-consistent implementation of effective 
nuclear interactions: both the equations which determine the ground state of a nucleus, 
and the residual interaction which governs the small-amplitude vibrations around the 
equilibrium, must be derived from the same effective interactions in the particle-hole and 
particle-particle channels. Only fully self-consistent calculations ensure the separation 
of spurious states from the physical excitations, and provide reliable microscopic 
predictions of collective modes directly based on ground-state properties. The presently 
available models which meet these conditions are: the HFB plus QRPA based on 
Skyrme or Gogny effective interactions, and in the relativistic framework the RHB 
plus RQRPA based on effective Lagrangians with density-dependent meson-nucleon 
couplings. Even though some of these state-of-the-art models include the coupling to the 
particle continuum, i.e. they take into account the escape width of vibrational states 
that lie above the particle emission threshold, one really must go beyond the mean- 
field approximation in order to quantitatively describe decay properties of collective 
excitations. At present only few theoretical approaches take into account in a fully 
consistent way the coupling of the simple lp — lh (or two-quasiparticle) states to more 
complex 2p — 2h configurations, which basically determines the spreading width of 
resonances. The way the damping mechanism is modified in nuclei far from stability, 
e.g. neutron-rich nuclei, is therefore largely unexplored. 

Studies of low-energy collective excitations in nuclei with relatively moderate 
neutron excess provide crucial information on the manifestation of exotic modes in 
neutron-rich nuclei, and their subsequent evolution towards systems with large isospin 
asymmetry near the nucleon drip-lines. Evidence for the existence of low-lying dipole 
strength in neutron-rich nuclei, which might indeed correspond to the pygmy dipole 
resonance (PDR), has become available from several experiments based on (7,7') 
resonant scattering, and data on low-lying El strength in exotic nuclei has recently 
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been obtained in studies with radioactive ion beams. Experimentally, however, very 
little is known about the nature of the observed low-lying dipole states, the degree of 
collectivity and the isospin character, and extensive investigations have to be carried 
out. For instance, studies of low-energy transition strength below the particle emission 
threshold in (a, a'7) coincidence experiments, which enable a clear separation of El 
excitations from states of other multipolarities, will allow a more detailed analysis of 
the structure of these states [388J. 

A number of theoretical studies have recently been devoted to properties of the 
PDR, and some models predict an enhancement of collectivity for these states, and 
pronounced mixing of isoscalar and isovector components. An important issue in 
neutron-rich nuclei is also the relationship between the PDR excitation energy and 
the neutron separation threshold. It appears that in nuclei not so far from stability the 
PDR is located below or very close to the neutron separation energy, whereas in weakly- 
bound systems with large neutron excess the PDR energy is high above the neutron 
threshold. Only part of the low-lying El strength is observed in current experiments. 
While the (7, 7') spectra are composed of peaks below the threshold energy, the Coulomb 
excitation of fission fragments gives the transition strength above the neutron separation 
energy. Future experiments therefore need to provide complete low-lying dipole spectra, 
both below and above the neutron threshold energy. A comparison of the complete 
spectra of measured low-lying strength in nuclei far from stability with the predictions 
of self-consistent microscopic theories, will present a very sensitive test for the isovector 
channel of nuclear effective interactions. Namely, it is difficult to adjust the isovector 
terms of effective interactions only to data in stable nuclei. Stringent constraints on the 
microscopic approach to nuclear dynamics and effective nuclear interactions will emerge 
from studies of the structure and stability of exotic nuclei with extreme isospin values. 

Except for the PDR, few data have been reported so far on other possible exotic 
modes of excitation, but dedicated experiments are being planned and designed at 
radioactive-beam laboratories. For instance, in the near future more information could 
become available on modes which arise in nuclei near the drip-lines: the proton PDR, 
and the di-neutron vs core oscillations. An important question which has not been 
addressed so far is the evolution of low-energy modes with temperature in hot nuclei, a 
topic that could be important for astrophysical applications. The dipole toroidal mode 
could be probed by the measurement of transverse electron scattering form factors at 
180°. An entirely unexplored field is the evolution of modes that involve the spin and/or 
isospin degrees of freedom in nuclei far from stability. It must also be emphasized that 
while most theoretical studies of exotic modes of excitation have assumed spherical 
symmetry, the evolution of deformation in unstable neutron-rich nuclei could give rise 
to interesting collective phenomena such as, for instance, the low-energy scissors mode 
of oscillation of the neutron skin [389J . 

Of particular importance is the role that new exotic modes or, more generally 
collective excitations in nuclei far from stability, play in astrophysical processes: 
the description of Gamow- Teller resonances is essential in calculations of /3-decay, 
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electron capture and neutrino-nucleus interaction rates; the low-energy dipole transition 
strength in neutron-rich nuclei has a pronounced effect on the calculated r-process 
abundances and on the propagation of ultra-high energy cosmic rays; on the proton- 
rich side the proton pygmy dipole resonance could contribute to the nucleosynthesis in 
rapid proton capture processes, as well as in the two-proton capture in astrophysical 
conditions characteristic for explosive hydrogen burning in novae and x-ray bursts [390J . 
Theoretical studies of electron capture rates on neutron-rich nuclei, at temperatures and 
densities characteristic for core collapse, have recently shown that these rates can be so 
large that electron capture on nuclei dominates over capture on free protons. Further 
studies of capture rates, particularly in a self-consistent approach extended to finite 
temperatures, are clearly desirable. The effect of giant resonances in extremely neutron- 
rich systems on the cooling time of neutron stars should also be investigated. Neutrino- 
nucleus reactions in the low-energy range 1-100 MeV play an important role in many 
astrophysical processes, including stellar nucleosynthesis, and the study of low-energy 
neutrino reactions on medium- heavy and heavy nuclei is of great current interest. Since 
the description of a neutrino-nucleus reaction becomes increasingly complicated as the 
target mass number increases, accurate self-consistent mean-field approaches must be 
developed and applied in calculations of cross sections for all relevant neutrino-induced 
reactions. More generally, microscopic nuclear structure theory must be integrated into 
various astrophysical models of nucleosynthesis processes, supernova dynamics, and 
neutrino-induced reactions, by providing accurate global predictions for bulk nuclear 
properties and nuclear excitations. 
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